Bursts of gamma ray showers have been observed in coincidence with downward propagating negative leaders in lightning flashes by the Telescope Array Surface Detector (TASD). The TASD is a 700~square kilometer cosmic ray observatory located in southwestern Utah, U.S.A. In data collected between 2014 and 2016, correlated observations showing the structure and temporal development of three shower-producing flashes were obtained with a 3D lightning mapping array, and electric field change measurements were obtained for an additional seven flashes, in both cases co-located with the TASD. National Lightning Detection Network (NLDN) information was also used throughout. The showers arrived in a sequence of 2--5 short-duration ($\le$10~$\mu$s) bursts over time intervals of several hundred microseconds, and originated at an altitude of $\simeq$3--5 kilometers above ground level during the first 1--2 ms of downward negative leader breakdown at the beginning of cloud-to-ground lightning flashes. The shower footprints, associated waveforms and the effect of atmospheric propagation indicate that the showers consist primarily of downward-beamed gamma radiation. This has been supported by GEANT simulation studies, which indicate primary source fluxes of $\simeq$$10^{12}$--$10^{14}$ photons for $16^{\circ}$ half-angle beams. We conclude that the showers are terrestrial gamma-ray flashes (TGFs), similar to those observed by satellites, but that the ground-based observations are more representative of the temporal source activity and are also more sensitive than satellite observations, which detect only the most powerful TGFs.

We report an in-depth Angle Resolved Photoemission Spectroscopy (ARPES) study on $2H$-TaS$_2$, a canonical incommensurate Charge Density Wave (CDW) system. This study demonstrates that just as in related incommensurate CDW systems, $2H$-TaSe$_2$ and $2H$-NbSe$_2$, the energy gap ($\Delta_{\text{cdw}}\,$) of $2H$-TaS$_2$ is localized along the K-centered Fermi surface barrels and is particle-hole asymmetric. The persistence of $\Delta_{\text{cdw}}\,$ even at temperatures higher than the CDW transition temperature $\it{T}_{\text{cdw}}\,$ in $2H$-TaS$_2$, reflects the similar pseudogap (PG) behavior observed previously in $2H$-TaSe$_2$ and $2H$-NbSe$_2$. However, in sharp contrast to $2H$-NbSe$_2$, where $\Delta_{\text{cdw}}\,$ is non-zero only in the vicinity of a few "hot spots" on the inner K-centered Fermi surface barrels, $\Delta_{\text{cdw}}\,$ in $2H$-TaS$_2$ is non-zero along the entirety of both K-centered Fermi surface barrels. Based on a tight-binding model, we attribute this dichotomy in the momentum dependence and the Fermi surface specificity of $\Delta_{\text{cdw}}\,$ between otherwise similar CDW compounds to the different orbital orientations of their electronic states that are involved in CDW pairing. Our results suggest that the orbital selectivity plays a critical role in the description of incommensurate CDW materials.

Photons, electrons, and their interplay are at the heart of photonic devices and modern instruments for ultrafast science [1-10]. Nowadays, electron beams of the highest intensity and brightness are created by photoemission with short laser pulses, and then accelerated and manipulated using GHz radiofrequency electromagnetic fields. The electron beams are utilized to directly map photoinduced dynamics with ultrafast electron scattering techniques, or further engaged for coherent radiation production at up to hard X-ray wavelengths [11-13]. The push towards improved timing precision between the electron beams and pump optical pulses though, has been stalled at the few tens of femtosecond level, due to technical challenges with synchronizing the high power rf fields with optical sources. Here, we demonstrate attosecond electron metrology using laser-generated single-cycle THz radiation, which is intrinsically phase locked to the optical drive pulses, to manipulate multi-MeV relativistic electron beams. Control and single-shot characterization of bright electron beams at this unprecedented level open up many new opportunities for atomic visualization.

Non-stoichiometry effect on the extreme magnetoresistance is systematically investigated for the Weyl semimetal WTe2. Magnetoresistance and Hall resistivity are measured for the as-grown samples with a slight difference in Te vacancies and the annealed samples with increased Te vacancies. The fittings to a two-carrier model show that the magnetoresistance is strongly dependent on the residual resistivity ratio (i.e., the degree of non-stoichiometry), which is eventually understood in terms of electron doping which not only breaks the balance between electron-type and hole-type carrier densities but also reduces the average carrier mobility. Thus, compensation effect and ultrahigh mobility are probably the main driving force of the extreme magnetoresistance in WTe2.

We consider a variety of settings involving chains with one or more defects stemming from the introduction of nodes bearing internal resonators. Motivated by experimental results in woodpile elastic lattices with one or two defects, we consider a variety of different theoretical scenarios. These include multi-defect chains and their ability to transmit, reflect, and especially trap energy; they also include settings with linear vs. nonlinear defects of variable interaction exponent. Moreover, they involve defects which are spatially separated and either statically, or more effectively dynamically, enable the confinement of energy between the separated defects. Wherever possible, comparisons of the experiments with numerical simulations, as well as with theoretical intuition are also offered, to provide a justification for the observed findings.

To date, PMN J2134-0419 (at a redshift z=4.33) is the second most distant quasar known with a milliarcsecond-scale morphology permitting direct estimates of the jet proper motion. Based on two-epoch observations, we constrained its radio jet proper motion using the very long baseline interferometry (VLBI) technique. The observations were conducted with the European VLBI Network (EVN) at 5 GHz on 1999 November 26 and 2015 October 6. We imaged the central 10-pc scale radio jet emission and modeled its brightness distribution. By identifying a jet component at both epochs separated by 15.86 yr, a proper motion of mu=0.035 +- 0.023 mas/yr is found. It corresponds to an apparent superluminal speed of beta_a=4.1 +- 2.7 c . Relativistic beaming at both epochs suggests that the jet viewing angle with respect to the line of sight is smaller than 20 deg, with a minimum bulk Lorentz factor Gamma=4.3. The small value of the proper motion is in good agreement with the expectations from the cosmological interpretation of the redshift and the current cosmological model. Additionally we analyzed archival Very Large Array observations of J2143-0419 and found indication of a bent jet extending to ~30 kpc.

The correlation of co-located hydrogen masers (H-masers) is difficult to measure because their common-mode noise induced by the environment will be cancelled out during the comparison measurement. With the development of fibre-based high-precision time and frequency transfer technique, the correlation of co-located hydrogen masers can be directly measured with the help of remote H-masers. Recently, a fiber-based frequency synchronization network was constructed in the Beijing region by connecting 5 H-masers from 4 institutions. The correlation coefficient of atomic clocks is defined and the correlation between two co-located H-masers is measured using both experimental and simulative methods. The results show that the correlation is not prominent until the averaging time is larger than $\sim10^3$s; then, the coefficient grows rapidly for averaging times ranging from $\sim10^3$s to $\sim10^5$s and decreases beyond $\sim10^5$s up to 5 days.

We present a 45 ks Chandra observation of the quasar ULAS J1342+0928 at z=7.54. We detect 14.0^{+4.8}_{-3.7} counts from the quasar in the observed-frame energy range 0.5-7.0 keV (6-sigma detection), representing the most distant non-transient astronomical source identified in X-rays to date. The present data are sufficient only to infer rough constraints on the spectral parameters. We find an X-ray hardness ratio of HR = -0.51^{+0.26}_{-0.28} between the 0.5-2.0 keV and 2.0-7.0 keV ranges and derive a power-law photon index of Gamma = 1.95^{+0.55}_{-0.53}. Assuming a typical value for high-redshift quasars of Gamma = 1.9, ULAS J1342+0928 has a 2-10 keV rest-frame X-ray luminosity of L_{2-10} = 11.6^{+4.3}_{-3.5} x 10^{44} erg/s. Its X-ray-to-optical power-law slope is alpha_{OX}=-1.67^{+0.16}_{-0.10}, consistent with the general trend indicating that the X-ray emission in the most bolometrically powerful quasars is weaker relative to their optical emission.

The fractional or non-integer valence state in 3d and 4d transition metal oxides always gives rise to novel physics and multifunctional applications, such as charge order, magnetoelectric effect, and double exchange. The ground state becomes even more complex for 5d systems, as a result of the enhanced spin-orbit couplings. In this paper, experimental and theoretical studies are performed on La3Ir3O11, a new iridate with a non-1/2 fractional valence state (+4.33) for Ir ions and a three-dimensional strongly frustrated structure composed of edge-shared Ir2O10 dimers. Density functional theory calculations explain the spin-orbit driven Mott insulating transport, and show that the electronic state is neither Jeff=1/2 state nor J=0 state. Such a peculiar electronic state leads to a novel magnetic ground state, which hosts strong antiferromagnetic fluctuation yet without a long-range magnetic order. These behaviors are eventually demonstrated by our theoretical model beyond the simple Curie-Weiss law. Our work would shed new light on the studies of quantum spin liquid.

We report on a measurement of the cosmic ray energy spectrum by the Telescope Array Low-Energy Extension (TALE) air fluorescence detector. The TALE air fluorescence detector is also sensitive to the Cherenkov light produced by shower particles. Low energy cosmic rays, in the PeV energy range, are detectable by TALE as "Cherenkov Events". Using these events, we measure the energy spectrum from a low energy of $\sim 2$ PeV to an energy greater than 100 PeV. Above 100 PeV TALE can detect cosmic rays using air fluorescence. This allows for the extension of the measurement to energies greater than a few EeV. In this paper, we will describe the detector, explain the technique, and present results from a measurement of the spectrum using $\sim 1000$ hours of observation. The observed spectrum shows a clear steepening near $10^{17.1}$ eV, along with an ankle-like structure at $10^{16.2}$ eV. These features present important constraints on galactic cosmic rays origin and propagation models. The feature at $10^{17.1}$ eV may also mark the end of the galactic cosmic rays flux and the start of the transition to extra-galactic sources.

An intermediate-scale energy spectrum anisotropy has been found in the arrival directions of ultra-high energy cosmic rays of energies above $10^{19.2}$ eV in the northern hemisphere, using 7 years of data from the Telescope Array surface detector. A relative energy distribution test is done comparing events inside oversampled spherical caps of equal exposure, to those outside, using the Poisson likelihood ratio. The center of maximum significance is at $9^h$$16^m$, $45^{\circ}$. and has a deficit of events with energies $10^{19.2}$$\leq$$E$$<$$10^{19.75}$ eV and an excess for $E$$\geq$$10^{19.75}$ eV. The post-trial probability of this energy anisotropy, appearing by chance anywhere on an isotropic sky, is found by Monte Carlo simulation to be $9$$\times$$10^{-5}$ ($3.74$$\sigma_{global}$).

Radio-loud active galactic nuclei (AGNs), hosting powerful relativistic jet outflows, provide an excellent laboratory for studying jet physics. Very long baseline interferometry (VLBI) enables high-resolution imaging on milli-arcsecond (mas) and sub-mas scales, making it a powerful tool to explore the inner jet structure, shedding light on the formation, acceleration and collimation of AGN jets. In this paper, we present Very Long Baseline Array (VLBA) observations of ten radio-loud AGNs at 43 and 86~GHz, which were selected from the {\it Planck} catalogue of compact sources and are among the brightest in published VLBI images at and below 15 GHz. The image noise levels in our observations are typically 0.3 mJy beam$^{-1}$ and 1.5 mJy beam$^{-1}$ at 43 and 86 GHz, respectively. Compared with the VLBI data observed at lower frequencies from the literature, our observations with higher resolution (the highest resolution up to 0.07 mas at 86 GHz and 0.18 mas at 43 GHz) and at higher frequencies detected new jet components at sub-parsec scales, offering valuable data for studies of the physical properties of innermost jets. These include compactness factor of the radio structure (the ratio of core flux density to total flux density), and core brightness temperature ($T_{\rm b}$). In all these sources, the compact core accounts for a significant fraction ($> 60\%$) of the total flux density. Their correlated flux density at the longest baselines is higher than 0.16 Jy. The compactness of these sources make them good phase calibrators of mm-wavelength ground-based and space VLBI.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

Glassiness is ubiquitous and diverse in characteristics in nature. Understanding their differences and classification remains a major scientific challenge. Here, we show that scaling of magnetic memories with time can be used to classify magnetic glassy materials into two distinct classes. The systems studied are high-temperature superconductor-related materials, spin-orbit Mott insulators, frustrated magnets, and dilute magnetic alloys. Our bulk magnetization measurements reveal that most densely populated magnets exhibit similar memory behavior characterized by a relaxation exponent of 1-n ~ 0.6(1). This exponent is different from 1-n ~ 1/3 of dilute magnetic alloys that was ascribed to their hierarchical and fractal energy landscape and is also different from 1-n=1 of the conventional Debye relaxation expected for a spin solid, a state with long range order. Furthermore, our systematic study on dilute magnetic alloys with varying magnetic concentration exhibits crossovers among the two glassy states and spin solid.

The Telescope Array (TA) experiment is located in the western desert of Utah, USA, and observes ultra high energy cosmic rays (UHECRs) in the Northern hemisphere. At the highest energies, $E>10$~EeV, the shape of cosmic ray energy spectrum may carry an imprint of the source density distribution along the line of sight different in different directions of the sky. In this study, we search for such directional variations in the shape of the energy spectrum using events observed with the Telescope Array's surface detector. We divide the TA field of view into two nearly equal-exposure regions: the "on-source" region which we define as $\pm 30^\circ$ of the supergalactic plane containing mostly nearby structures, and the complementary "off-source" region where the sources are further away on average. We compare the UHECR spectra in these regions by fitting them to the broken power law and comparing the resulting parameters. We find that the off-source spectrum has an earlier break at highest energies. The chance probability to obtain such or larger difference in statistically equivalent distributions is estimated as $6.2\pm1.1\times10^{-4}$ ($3.2\sigma$) by a Monte-Carlo simulation. The observed difference in spectra is in a reasonable quantitative agreement with a simplified model that assumes that the UHECR sources trace the galaxy distribution from the 2MRS catalogue, primary particles are protons and the magnetic deflections can be neglected.

We report on the tuning of magnetic interactions in superlattices composed of single and bilayer SrIrO$_3$ inter-spaced with SrTiO$_3$. Magnetic scattering shows predominately $c$-axis antiferromagnetic orientation of the magnetic moments for the bilayer justifying these systems as viable artificial analogues of the bulk Ruddlesden-Popper series iridates. Magnon gaps are observed in both superlattices, with the magnitude of the gap in the bilayer being reduced to nearly half that in its bulk structural analogue, Sr$_3$Ir$_2$O$_7$. We assign this to modifications in the anisotropic exchange driven by bending of the $c$-axis Ir-O-Ir bond and subsequent local environment changes, as detected by x-ray diffraction and modeled using spin wave theory. These findings explain how even subtle structural modulations driven by heterostructuring in iridates are leveraged by spin orbit coupling to drive large changes in the magnetic interactions.

We report $^{75}$As nuclear magnetic resonance (NMR) / nuclear quadrupole resonance (NQR) and transmission electron microscopy (TEM) studies on LaFeAsO$_{1-x}$F$_{x}$. There are two superconducting domes in this material. The first one appears at 0.03 $\leq$ $x$ $\leq$ 0.2 with $T_{\rm c}$$^{max}$ = 27 K, and the second one at 0.25 $\leq$ $x$ $\leq$ 0.75 with $T_{\rm c}$$^{max}$ = 30 K. By NMR and TEM, we demonstrate that a $C4$-to-$C2$ structural phase transition (SPT) takes place above both domes, with the transition temperature $T_{\rm s}$ varying strongly with $x$. In the first dome, the SPT is followed by an antiferromagnetic (AF) transition, but neither AF order nor low-energy spin fluctuations are found in the second dome. In LaFeAsO$_{0.97}$F$_{0.03}$, we find that AF order and superconductivity coexist microscopically via $^{75}$As nuclear spin-lattice relaxation rate (1/$T_1$) measurements. In the coexisting region, 1/$T_1$ decreases at $T_{\rm c}$ but becomes to be proportional to $T$ below 0.6$T_{\rm c}$, indicating gapless excitations. Therefore, in contrast to the early reports, the obtained phase diagram for $x \leq$ 0.2 is quite similar to the doped BaFe$_{2}$As$_{2}$ system. The electrical resistivity in the second dome can be fitted by $\rho = {{\rho }_{0}}+A{{T}^{n}}$ with $n$ = 1 and a maximal coefficient $A$ at around $x_{opt}$ = 0.5$\sim$0.55 where $T_{\rm s}$ extrapolates to zero and $T_{\rm c}$ is the maximal, which suggest the importance of quantum critical fluctuations associated with the SPT. We have constructed a complete phase diagram of LaFeAsO$_{1-x}$F$_{x}$, which provides insight into the relationship between SPT, antiferromagnetism and superconductivity.

In September 2016, the microquasar Cygnus X-3 underwent a giant radio flare, which was monitored for 6 days with the Medicina Radio Astronomical Station and the Sardinia Radio Telescope. Long observations were performed in order to follow the evolution of the flare on a hourly scale, covering six frequency ranges from 1.5 GHz to 25.6 GHz. The radio emission reached a maximum of 13.2 +/- 0.7 Jy at 7.2 GHz and 10 +/- 1 Jy at 18.6 GHz. Rapid flux variations were observed at high radio frequencies at the peak of the flare, together with rapid evolution of the spectral index: alpha steepened from 0.3 to 0.6 within 5 hours. This is the first time that such fast variations are observed, giving support to the evolution from optically thick to optically thin plasmons in expansion moving outward from the core. Based on the Italian network (Noto, Medicina and SRT) and extended to the European antennas (Torun, Yebes, Onsala), VLBI observations were triggered at 22 GHz on five different occasions, four times prior to the giant flare, and once during its decay phase. Flux variations of 2-hour duration were recorded during the first session. They correspond to a mini-flare that occurred close to the core ten days before the onset of the giant flare. From the latest VLBI observation we infer that four days after the flare peak the jet emission was extended over 30 mas.

The number of measurements required to reconstruct the states of quantum systems increases exponentially with the quantum system dimensions, which makes the state reconstruction of high-qubit quantum systems have a great challenge in physical quantum computing experiments. Compressive sensing (CS) has been verified as a effective technique in the reconstruction of quantum state, however, it is still unknown that if CS can reconstruct quantum states given the less data measured by nuclear magnetic resonance (NMR). In this paper, we propose an effective NMR quantum state reconstruction method based on CS. Different from the conventional CS-based quantum state reconstruction, our method uses the actual observation data from NMR experiments rather than the data measured by the Pauli operators. We implement measurements on quantum states in practical NMR computing experiments and reconstruct states of 2,3,4 qubits using fewer number of measurements, respectively. The proposed method is easy to implement and performs more efficiently with the increase of the system dimension size. The performance reveals both efficiency and accuracy, which provides an alternative for the quantum state reconstruction in practical NMR.

A small fraction of Tidal Disruption Events (TDE) produce relativistic jets, evidenced by their non-thermal X-ray spectra and transient radio emission. Here we present milliarcsecond-resolution imaging results on TDE J1644+5734 with the European VLBI Network (EVN). These provide a strong astrometric constraint on the average apparent jet velocity <0.27, that constrains the intrinsic jet velocity for a given viewing angle.

We have collected and analyzed the complete archive of {\itshape XMM-Newton\} (116), {\itshape Chandra\} (151), and {\itshape RXTE\} (952) observations of the Small Magellanic Cloud (SMC), spanning 1997-2014. The resulting observational library provides a comprehensive view of the physical, temporal and statistical properties of the SMC pulsar population across the luminosity range of $L_X= 10^{31.2}$--$10^{38}$~erg~s$^{-1}$. From a sample of 67 pulsars we report $\sim$1654 individual pulsar detections, yielding $\sim$1260 pulse period measurements. Our pipeline generates a suite of products for each pulsar detection: spin period, flux, event list, high time-resolution light-curve, pulse-profile, periodogram, and spectrum. Combining all three satellites, we generated complete histories of the spin periods, pulse amplitudes, pulsed fractions and X-ray luminosities. Some pulsars show variations in pulse period due to the combination of orbital motion and accretion torques. Long-term spin-up/down trends are seen in 12/11 pulsars respectively, pointing to sustained transfer of mass and angular momentum to the neutron star on decadal timescales. Of the sample 30 pulsars have relatively very small spin period derivative and may be close to equilibrium spin. The distributions of pulse-detection and flux as functions of spin-period provide interesting findings: mapping boundaries of accretion-driven X-ray luminosity, and showing that fast pulsars ($P<$10 s) are rarely detected, which yet are more prone to giant outbursts. Accompanying this paper is an initial public release of the library so that it can be used by other researchers. We intend the library to be useful in driving improved models of neutron star magnetospheres and accretion physics.

Using hundreds of XMM-Newton and Chandra archival observations and nearly a thousand RXTE observations, we have generated a comprehensive library of the known pulsars in the Small and Large Magellanic Clouds (SMC, LMC). The pulsars are detected multiple times across the full parameter spaces of X-ray luminosity ($L_X= 10^{31-38}$~erg/s) and spin period ( P$<$1s -- P$>$1000s) and the library enables time-domain studies at a range of energy scales. The high time-resolution and sensitivity of the EPIC cameras are complemented by the angular resolution of Chandra and the regular monitoring of RXTE. Our processing %$\sim$15 year pipeline uses the latest calibration files and software to generate a suite of useful products for each pulsar detection: event lists, high time-resolution light curves, periodograms, spectra, and complete histories of $\dot{P}$, the pulsed fraction, etc., in the broad (0.2-12 keV), soft (0.2-2 keV), and hard (2-12 keV) energy bands. After combining the observations from these telescopes, we found that 28 pulsars show long-term spin up and 25 long-term spin down. We also used the faintest and brightest sources to map out the lower and upper boundaries of accretion-powered X-ray emission: the propeller line and the Eddington line, respectively. We are in the process of comparing the observed pulse profiles to geometric models of X-ray emission in order to constrain the physical parameters of the pulsars. Finally we are preparing a public release of the library so that it can be used by others in the astronomical community.

Black hole X-ray binaries undergo occasional outbursts caused by changing inner accretion flows. Here we report high-angular resolution radio observations of the 2013 outburst of the black hole candidate X-ray binary system J1908+094, using data from the VLBA and EVN. We show that following a hard-to-soft state transition, we detect moving jet knots that appear asymmetric in morphology and brightness, and expand to become laterally resolved as they move away from the core, along an axis aligned approximately $-11$\degree\ east of north. We initially see only the southern component, whose evolution gives rise to a 15-mJy radio flare and generates the observed radio polarization. This fades and becomes resolved out after 4 days, after which a second component appears to the north, moving in the opposite direction. From the timing of the appearance of the knots relative to the X-ray state transition, a 90\degree\ swing of the inferred magnetic field orientation, the asymmetric appearance of the knots, their complex and evolving morphology, and their low speeds, we interpret the knots as working surfaces where the jets impact the surrounding medium. This would imply a substantially denser environment surrounding XTE J1908+094 than has been inferred to exist around the microquasar sources GRS 1915+105 and GRO J1655-40.

We study the dynamics of isospectral granular chains that are highly tunable due to the nonlinear Hertz contact law interaction between the granular particles. The system dynamics can thus be tuned easily from being linear to strongly nonlinear by adjusting the initial compression applied to the chain. In particular, we introduce both discrete and continuous spectral transformation schemes to generate a family of granular chains that are isospectral in their linear limit. Inspired by the principle of supersymmetry in quantum systems, we also introduce a methodology to add or remove certain eigenfrequencies, and we demonstrate numerically that the corresponding physical system can be constructed in the setting of one-dimensional granular crystals. In the linear regime, we highlight the commonalities and differences in the elastic wave transmission characteristics of such isospectral systems, and emphasize that the presented mathematical framework allows one to suitably tailor the wave transmission through a general class of granular chains. Moreover, we show how the dynamic response of these structures deviates from its linear limit as we introduce Hertzian nonlinearity in the chain and how nonlinearity breaks the notion of linear isospectrality.