Source author record

Y. Su

Y. Su appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

70works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

70 published item(s)

preprint2026arXiv

An Enigmatic PeVatron in an Area around HII Region G35.6$-$0.5

Identifying Galactic PeVatrons (PeV particle accelerators) from the ultra-high-energy (UHE, >100 TeV) $γ$-ray sources plays a crucial role in revealing the origin of Galactic cosmic rays. The UHE source 1LHAASO J1857+0203u is suggested to be associated with HESS J1858+020, which may be attributed to the possible PeVatron candidate supernova remnant (SNR) G35.6$-$0.4 or HII region G35.6$-$0.5. We perform detailed analysis on the very-high-energy and UHE $γ$-ray emissions towards this region with data from the Large High Altitude Air Shower Observatory (LHAASO). 1LHAASO J1857+0203u is detected with a significance of 11.6$σ$ above 100 TeV, indicating the presence of a PeVatron. It has an extension of $\sim 0.18^\circ$ with a power-law (PL) spectral index of $\sim$2.5 in 1-25 TeV and a point-like emission with a PL spectral index of $\sim$3.2 above 25 TeV. Using the archival CO and HI data, we identify some molecular and atomic clouds that may be associated with the TeV $γ$-ray emissions. Our modelling indicates that the TeV $γ$-ray emissions are unlikely to arise from the clouds illuminated by the protons that escaped from SNR G35.6$-$0.4. In the scenario that HII region G35.6$-$0.5 could accelerate particles to the UHE band, the observed GeV-TeV $γ$-ray emission could be well explained by a hadronic model with a PL spectral index of $\sim$2.0 and cutoff energy of $\sim$450 TeV. However, an evolved pulsar wind nebula origin cannot be ruled out.

preprint2026arXiv

An Ultrahigh-energy $γ$-ray Bubble Powered by a Super PeVatron

We report the detection of a $γ$-ray bubble spanning at least 100$\rm deg^2$ in ultra high energy (UHE) up to a few PeV in the direction of the star-forming region Cygnus X, implying the presence Super PeVatron(s) accelerating protons to at least 10 PeV. A log-parabola form with the photon index $Γ(E) = (2.71 \pm 0.02) + (0.11 \pm 0.02) \times \log_{10} (E/10 \ {\rm TeV})$ is found fitting the gamma-ray energy spectrum of the bubble well. UHE sources, `hot spots' correlated with very massive molecular clouds, and a quasi-spherical amorphous $γ$-ray emitter with a sharp central brightening are observed in the bubble. In the core of $\sim 0.5^{\circ}$, spatially associating with a region containing massive OB association (Cygnus OB2) and a microquasar (Cygnus X-3), as well as previously reported multi-TeV sources, an enhanced concentration of UHE $γ$-rays are observed with 2 photons at energies above 1 PeV. The general feature of the bubble, the morphology and the energy spectrum, are reasonably reproduced by the assumption of a particle accelerator in the core, continuously injecting protons into the ambient medium.

preprint2026arXiv

Constraining the Cosmic-ray Energy Based on Observations of Nearby Galaxy Clusters by LHAASO

Galaxy clusters act as reservoirs of high-energy cosmic rays (CRs). As CRs propagate through the intracluster medium, they generate diffuse $γ$-rays detectable by arrays such as LHAASO. These $γ$-rays result from proton-proton ($pp$) collisions of very high-energy cosmic rays (VHECRs) or inverse Compton (IC) scattering of positron-electron pairs created by $pγ$ interactions of ultra-high-energy cosmic rays (UHECRs). We analyzed diffuse $γ$-ray emission from the Coma, Perseus, and Virgo clusters using LHAASO data. Diffuse emission was modeled as a disk of radius $R_{500}$ for each cluster while accounting for point sources. No significant diffuse emission was detected, yielding 95\% confidence level (C.L.) upper limits on the $γ$-ray flux: for WCDA (1-25~TeV) and KM2A ($>25$~TeV), less than $(49.4, 13.7, 54.0)$ and $(1.34, 1.14, 0.40) \times 10^{-14}$~ph~cm$^{-2}$~s$^{-1}$ for Coma, Perseus, and Virgo, respectively. The $γ$-ray upper limits can be used to derive model-independent constraints on the integral energy of CRp above 10~TeV (corresponding to the LHAASO observational range $>1$~TeV under the $pp$ scenario) to be less than $(1.96, 0.59, 0.08) \times 10^{61}$~erg. The absence of detectable annuli/ring-like structures, indicative of cluster accretion or merging shocks, imposes further constraints on models in which the UHECRs are accelerated in the merging shocks of galaxy clusters.

preprint2026arXiv

Constraints on heavy decaying dark matter from 570 days of LHAASO observations

The Kilometer Square Array~(KM2A) of the Large High Altitude Air Shower Observatory (LHAASO) aims at surveying the northern gamma-ray sky at energies above 10 TeV with unprecedented sensitivity. Gamma-ray observations have long been one of the most powerful tools for dark matter searches, as e.g., high-energy gamma-rays could be produced by the decays of heavy dark matter particles. In this letter, we present the first dark matter analysis with LHAASO-KM2A, using the first 340~days of data from 1/2-KM2A and 230~days of data from 3/4-KM2A. Several regions of interest are used to search for a signal and account for the residual cosmic-ray background after gamma/hadron separation. We find no excess of dark matter signals, and thus place some of the strongest gamma-ray constraints on the lifetime of heavy dark matter particles with mass between 10^5 and 10^9~GeV. Our results with LHAASO are robust, and have important implications for dark matter interpretations of the diffuse astrophysical high-energy neutrino emission.

preprint2026arXiv

Deep view of Composite SNR CTA1 with LHAASO in $γ$-rays up to 300 TeV

The ultra-high-energy (UHE) gamma-ray source 1LHAASO J0007+7303u is positionally associated with the composite SNR CTA1 that is located at high Galactic Latitude $b\approx 10.5^\circ$. This provides a rare opportunity to spatially resolve the component of the pulsar wind nebula (PWN) and supernova remnant (SNR) at UHE. This paper conducted a dedicated data analysis of 1LHAASO J0007+7303u using the data collected from December 2019 to July 2023. This source is well detected with significances of 21$σ$ and 17$σ$ at 8$-$100 TeV and $>$100 TeV, respectively. The corresponding extensions are determined to be 0.23$^{\circ}\pm$0.03$^{\circ}$ and 0.17$^{\circ}\pm$0.03$^{\circ}$. The emission is proposed to originate from the relativistic electrons and positrons accelerated within the PWN of PSR J0007+7303. The energy spectrum is well described by a power-law with an exponential cutoff function $dN/dE = (42.4\pm4.1)(\frac{E}{20\rm\ TeV})^{-2.31\pm0.11}\exp(-\frac{E}{110\pm25\rm\ TeV})$ $\rm\ TeV^{-1}\ cm^{-2}\ s^{-1}$in the energy range from 8 TeV to 300 TeV, implying a steady-state parent electron spectrum $dN_e/dE_e\propto (\frac{E_e}{100\rm\ TeV})^{-3.13\pm0.16}\exp[(\frac{-E_e}{373\pm70\rm\ TeV})^2]$ at energies above $\approx 50 \rm\ TeV$. The cutoff energy of the electron spectrum is roughly equal to the expected current maximum energy of particles accelerated at the PWN terminal shock. Combining the X-ray and gamma-ray emission, the current space-averaged magnetic field can be limited to $\approx 4.5\rm\ μG$. To satisfy the multi-wavelength spectrum and the $γ$-ray extensions, the transport of relativistic particles within the PWN is likely dominated by the advection process under the free-expansion phase assumption.

preprint2026arXiv

Discovery of a new $γ$-ray source LHAASO J0341+5258 with emission up to 200TeV

We report the discovery of a new unidentified extended $γ$-ray source in the Galactic plane named LHAASO J0341+5258 with a pre-trial significance of 8.2 standard deviations above 25 TeV. The best fit position is R.A.$=55.34^{\circ}\pm0.11^{\circ}$ and Dec$=52.97^{\circ}\pm0.07^{\circ}$. The angular size of LHAASO J0341+5258 is $0.29^\circ \pm 0.06^\circ_{stat} \pm0.02^\circ_{sys}$. The flux above 25 TeV is about $20\%$ of the flux of Crab Nebula. Although a power-law fit of the spectrum from 10 TeV to 200 TeV with the photon index $α=2.98 \pm 0.19_{stat} \pm 0.02_{sys}$ is not excluded, the LHAASO data together with the flux upper limit at 10 GeV set by the Fermi LAT observation, indicate a noticeable steepening of an initially hard power-law spectrum %($α\leq 1.75$) spectrum with a cutoff at $\approx 50$ TeV. We briefly discuss the origin of UHE gamma-rays. The lack of an energetic pulsar and a young SNR inside or in the vicinity of LHAASO J0341+5258 challenge, but do not exclude both the leptonic and hadronic scenarios of gamma-ray production.

preprint2026arXiv

Discovery of the Ultra-high energy gamma-ray source LHAASO J2108+5157

We report the discovery of a UHE gamma-ray source, LHAASO J2108+5157, by analyzing the LHAASO-KM2A data of 308.33 live days. Significant excess of gamma-ray induced showers is observed in both energy bands of 25-100 TeV and $\gt$100 TeV with 9.5 sigma and 8.5 sigma, respectively. This source is not significantly favored as an extensive source with the angular extension smaller than the point-spread function of KM2A. The measured energy spectrum from 20 to 200 TeV can be approximately described by a power-law function with an index of -2.83$\pm$ 0.18stat. A harder spectrum is demanded at lower energies considering the flux upper limit set by Fermi-LAT observations. The position of the gamma-ray emission is correlated with a giant molecular cloud, which favors a hadronic origin. No obvious counterparts have been found, deeper multiwavelength observations will help to shed new light on this intriguing UHE source.

preprint2026arXiv

Energy calibration of LHAASO-KM2A using the cosmic ray Moon shadow

We present a precise measurement of the westward, rigidity-dependent shift of the Moon's shadow using three and a half years of cosmic-ray data collected by the Kilometer Square Array (KM2A) of the Large High Altitude Air Shower Observatory (LHAASO). These measurements enable us to calibrate the detector energy response in the range 20-260 TeV, with results showing excellent agreement with the response derived from Monte Carlo (MC) simulations of the KM2A detector. We also measure a best-fit parameter $ε= 0.015 \pm 0.08$, corresponding to a 95% confidence interval of [-14%, +17%] for the energy-scale estimation. This result establishes the exceptional accuracy of the KM2A-MC in simulating the detector's response within this energy range.

preprint2026arXiv

Energy-Dependent Shifts of Medium-Scale Anisotropies in Very-High-Energy Cosmic Rays Observed by LHAASO-KM2A

Small deviations from isotropy in the arrival directions of Galactic cosmic rays serve as a unique probe of the local magnetic environment. In this Letter, we report observations of medium-scale anisotropies (MSA) at energies above 10 TeV using the LHAASO-KM2A array. Our analysis identifies four regions of excess and four regions of deficit, each spanning angular scales of approximately ten degrees. Crucially, we detect significant energy-dependent shifts in the centroids of two excess regions: Region B and the newly identified Region $\mathrm{\widetilde{D}}$. We also characterize the energy evolution of the fractional relative intensity across both excess and deficit regions. These findings imply that the observed anisotropies are shaped by the specific realization of the local turbulent magnetic field within the cosmic ray scattering length. Such energy-dependent behaviors impose strict constraints on local turbulence models and cosmic ray propagation theories.

preprint2026arXiv

Evidence for particle acceleration approaching PeV energies in the W51 complex

The $γ$-ray emission from the W51 complex is widely acknowledged to be attributed to the interaction between the cosmic rays (CRs) accelerated by the shock of supernova remnant (SNR) W51C and the dense molecular clouds in the adjacent star-forming region, W51B. However, the maximum acceleration capability of W51C for CRs remains elusive. Based on observations conducted with the Large High Altitude Air Shower Observatory (LHAASO), we report a significant detection of $γ$ rays emanating from the W51 complex, with energies from 2 TeV to 200 TeV. The LHAASO measurements, for the first time, extend the $γ$-ray emission from the W51 complex beyond 100 TeV and reveal a significant spectrum bending at tens of TeV. By combining the ``$π^0$-decay bump" featured data from Fermi-LAT, the broadband $γ$-ray spectrum of the W51 region can be well-characterized by a simple pp-collision model. The observed spectral bending feature suggests an exponential cutoff at $\sim400$~TeV or a power-law break at $\sim200$~TeV in the CR proton spectrum, most likely providing the first evidence of SNRs serving as CR accelerators approaching the PeV regime. Additionally, two young star clusters within W51B could also be theoretically viable to produce the most energetic $γ$ rays observed by LHAASO. Our findings strongly support the presence of extreme CR accelerators within the W51 complex and provide new insights into the origin of Galactic CRs.

preprint2026arXiv

Evidence for the rare decay $B^{+} \to \bar Λp μ^{+} μ^{-}$

A search for the rare decay $B^{+} \to \bar Λp μ^{+} μ^{-}$ is performed using proton-proton collision data recorded by the LHCb experiment at a center-of-mass energy of $\sqrt{s}= 13$ TeV, corresponding to an integrated luminosity of 5.4 fb$^-1$. An excess of events is found with respect to the background-only expectation, with a signal significance of 3.5 standard deviations, in the low invariant-mass region of $m(\bar Λp)<2.8$ GeV/$c^2$. The branching fraction is measured to be ${\cal B}_{\rm low}({B^{+}}\to {\bar Λp μ^{+} μ^{-}})=\left(1.70 ^{+0.65}_{-0.56}(\rm stat) \pm 0.17(\rm syst) \pm 0.14(\rm ext)\right) \times 10^{-8}$, where the last uncertainty is due to external inputs on ${\cal B}(B^+\to{J/ψ\bar Λp})\times {\cal B}({J/ψ}\to{μ^+μ^-})$. With no significant signal observed in the high $m(\bar Λp)$ region above 2.8 GeV/$c^2$, an upper limit is set to be ${\cal B}_{\rm high}(B^{+}\to {\bar Λp μ^{+} μ^{-}})<2.8\,(3.7) \times 10^{-9}$ at the $90\%$ ($95\%$) confidence level.

preprint2026arXiv

Exploring Lorentz Invariance Violation from Ultra-high-energy Gamma Rays Observed by LHAASO

Recently the LHAASO Collaboration published the detection of 12 ultra-high-energy gamma-ray sources above 100 TeV, with the highest energy photon reaching 1.4 PeV. The first detection of PeV gamma rays from astrophysical sources may provide a very sensitive probe of the effect of the Lorentz invariance violation (LIV), which results in decay of high-energy gamma rays in the superluminal scenario and hence a sharp cutoff of the energy spectrum. Two highest energy sources are studied in this work. No signature of the existence of LIV is found in their energy spectra, and the lower limits on the LIV energy scale are derived. Our results show that the first-order LIV energy scale should be higher than about 10^5 times the Planck scale M_{pl} and that the second-order LIV scale is >10^{-3}M_{pl}. Both limits improve by at least one order of magnitude the previous results.

preprint2026arXiv

Extended Very-High-Energy Gamma-Ray Emission Surrounding PSR J0622 + 3749 Observed by LHAASO-KM2A

We report the discovery of an extended very-high-energy (VHE) gamma-ray source around the location of the middle-aged (207.8 kyr) pulsar PSR J0622+3749 with the Large High Altitude Air Shower Observatory (LHAASO). The source is detected with a significance of $8.2σ$ for $E>25$~TeV assuming a Gaussian template. The best-fit location is (R.A., Dec.)$=(95^{\circ}\!.47\pm0^{\circ}\!.11,\,37^{\circ}\!.92 \pm0^{\circ}\!.09)$, and the extension is $0^{\circ}\!.40\pm0^{\circ}\!.07$. The energy spectrum can be described by a power-law spectrum with an index of ${-2.92 \pm 0.17_{\rm stat} \pm 0.02_{\rm sys} }$. No clear extended multi-wavelength counterpart of the LHAASO source has been found from the radio to sub-TeV bands. The LHAASO observations are consistent with the scenario that VHE electrons escaped from the pulsar, diffused in the interstellar medium, and scattered the interstellar radiation field. If interpreted as the pulsar halo scenario, the diffusion coefficient, inferred for electrons with median energies of $\sim160$~TeV, is consistent with those obtained from the extended halos around Geminga and Monogem and much smaller than that derived from cosmic ray secondaries. The LHAASO discovery of this source thus likely enriches the class of so-called pulsar halos and confirms that high-energy particles generally diffuse very slowly in the disturbed medium around pulsars.

preprint2026arXiv

LHAASO Detection of Ultra-High-Energy Gamma-Ray Emission toward the Giant Molecular Clouds

The $γ$-ray from Giant molecular clouds (GMCs) is regarded as the most ideal tool to perform in-situ measurement of cosmic ray (CR) density and spectra in our Galaxy. We report the first detection of $γ$-ray emissions in the very-high-energy (VHE) domain from the five nearby GMCs with a stacking analysis based on a 4.5-year $γ$-ray observation with the Large High Altitude Air Shower Observatory (LHAASO) experiment. The spectral energy distributions derived from the GMCs are consistent with the expected $γ$-ray flux produced via CR interacting with the ISM in the energy interval 1 - 100 $~\rm$ TeV. In addition, we investigate the presence of the CR spectral `knee' by introducing a spectral break in the $γ$-ray data. While no significant evidence for the CR knee is found, the current KM2A measurements from GMCs strongly favor a proton CR knee located above 0.9$~\rm$ PeV, which is consistent with the latest measurement of the CR spectrum by ground-based experiments.

preprint2026arXiv

Measurement of attenuation length of the muon content in extensive air showers from 0.3 to 30 PeV with LHAASO

The attenuation length of the muon content in extensive air showers provides important information regarding the generation and development of air showers. This information can be used not only to improve the description of such showers but also to test fundamental models of hadronic interactions. Using data from the LHAASO-KM2A experiment, the development of the muon content in high-energy air showers was studied. The attenuation length of muon content in the air showers was measured from experimental data in the energy range from 0.3 to 30 PeV using the constant intensity cut method. By comparing the attenuation length of the muon content with predictions from high-energy hadronic interaction models (QGSJET-II-04, SIBYLL 2.3d, and EPOS-LHC), it is evident that LHAASO results are significantly shorter than those predicted by the first two models (QGSJET-II-04 and SIBYLL 2.3d) but relatively close to those predicted by the third model (EPOS-LHC). Thus, the LHAASO data favor the EPOS-LHC model over the other two models. The three interaction models confirmed an increasing trend in the attenuation length as the cosmic-ray energy increases.

preprint2026arXiv

Measurement of transverse $Λ$ and $\barΛ$ hyperon polarization in $p$Pb collisions at $\sqrt{s_{NN}} = 5.02$ TeV

The transverse polarization of $Λ$ and $\barΛ$ hyperons is measured in $p$Pb collisions collected by the LHCb experiment at a nucleon-nucleon center-of-mass energy of $5.02$ TeV. The polarization is averaged over hyperon transverse momentum in the ranges $0.15 < p_{T} < 6.00 $ GeV/$c$, and Feynman-$x$ in the ranges $0.005 < x_{F} < 0.040$ (forward region) and $-0.10 < x_{F} < -0.01$ (backward region) defined relative to the proton beam direction. The transverse polarization is found to be compatible with zero for both $Λ$ and $\barΛ$ hyperons. The results are also measured as a function of $p_{T}$ and $x_{F}$ with no significant dependence on these variables observed. The results are compared with previous experimental measurements at different center-of-mass energies and collision environments.

preprint2026arXiv

Measurement of Very-high-energy Diffuse Gamma-ray Emissions from the Galactic Plane with LHAASO-WCDA

The diffuse Galactic gamma-ray emission is a very important tool used to study the propagation and interaction of cosmic rays in the Milky Way. In this work, we report the measurements of the diffuse emission from the Galactic plane, covering Galactic longitudes from $15^{\circ}$ to $235^{\circ}$ and latitudes from $-5^{\circ}$ to $+5^{\circ}$, in an energy range of 1 TeV to 25 TeV, with the Water Cherenkov Detector Array (WCDA) of the Large High Altitude Air Shower Observatory (LHAASO). After masking the sky regions of known sources, the diffuse emission is detected with $24.6σ$ and $9.1σ$ significance in the inner Galactic plane and outer Galactic plane, respectively. The WCDA spectra in both regions can be well described by a power-law function, with spectral indices of $-2.67\pm0.05_{\rm stat}$ in the inner region and $-2.83\pm0.19_{\rm stat}$ in the outer region, respectively. Combined with the Square Kilometer Array (KM2A) measurements at higher energies, a clear softening of the spectrum is found in the inner region, with change of spectral indices by $\sim0.5$ at a break energy around $30$ TeV. The fluxes of the diffuse emission are higher by a factor of $1.5-2.7$ than the model prediction assuming local CR spectra and the gas column density, which are consistent with those measured by the KM2A. Along Galactic longitude, the spatial distribution of the diffuse emission shows deviation from that of the gas column density. The spectral shape of the diffuse emission are possibly variation in different longitude region. The WCDA measurements bridge the gap between the low-energy measurements by space detectors and the ultra-high-energy observations by LHAASO-KM2A and other experiments. These results suggest that improved modeling of the wide-band diffuse emission is required.

preprint2026arXiv

Search for lepton-number-violating $B^-\to D^{(*)+}μ^-μ^-$ decays

A search is performed for lepton-number-violating $B^-\to D^{(*)+}μ^-μ^-$ decays, using data collected by the LHCb experiment in proton-proton collisions at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 5.4 fb$^{-1}$. No significant signal is observed, and upper limits are set on the branching fractions, ${\cal B}(B^-\to D^{+}μ^-μ^-) < 4.6 \times 10^{-8}$ and ${\cal B}(B^-\to D^{*+}μ^-μ^-) < 5.9 \times 10^{-8}$, at the 95% confidence level.

preprint2026arXiv

Study of charm mixing and CP violation with $D^0\to K^\pmπ^\mpπ^\pmπ^\mp$ decays

A study of charm mixing and CP violation in $D^0\to K^\pmπ^\mpπ^\pmπ^\mp$ decays is performed using data collected by the LHCb experiment in proton-proton collisions from 2015 to 2018, corresponding to an integrated luminosity of 6$\text{fb}^{-1}$. The ratio of promptly produced $D^0\to K^+π^- π^+π^-$ to $D^0\to K^-π^+ π^-π^+$ decay rates is measured as a function of $D^0$ decay time, both inclusive over phase space and in bins of phase space. Taking external inputs for the $D^0 -\overline{D}^0$ mixing parameters $x$ and $y$ allows constraints to be obtained on the hadronic parameters of the charm decay. When combined with previous measurements from charm-threshold experiments and at LHCb, improved knowledge is obtained for these parameters, which is valuable for studies of the angle $γ$ of the Unitarity Triangle. An alternative analysis is also performed, in which external inputs are taken for the hadronic parameters, and the mixing parameters are determined, including $Δx$ and $Δy$, which are nonzero in the presence of CP violation. It is found that $x=\left(0.85^{+0.15}_{-0.24}\right)\%$, $y=\left( 0.21^{+0.29}{-0.27} \right)\%$, $Δx=\left( -0.02\pm {0.04} \right)\% $ and $Δy=\left( 0.02^{+0.04}_{-0.03} \right)\%$. These results are consistent with previous measurements and the hypothesis of CP conservation.

preprint2026arXiv

Transient Large-Scale Anisotropy in TeV Cosmic Rays due to an Interplanetary Coronal Mass Ejection

Large- or medium-scale cosmic ray anisotropy at TeV energies has not previously been confirmed to vary with time. Transient anisotropy changes have been observed below 150 GeV, especially near the passage of an interplanetary shock and coronal mass ejection containing a magnetic flux rope ejected by a solar storm, which can trigger a geomagnetic storm with practical consequences. In such events, cosmic rays provide remote sensing of the magnetic field properties. Here we report the observation of transient large-scale anisotropy in TeV cosmic ray ions using data from the Large High Altitude Air Shower Observatory (LHAASO). We analyze hourly skymaps of the transient cosmic ray intensity excess or deficit, the gradient of which indicates the direction and magnitude of transient large-scale anisotropy across the field of view. We observe enhanced anisotropy above typical hourly fluctuations with $>$5$σ$ significance during some hours of November 4, 2021, in separate data sets for four primary cosmic ray energy ranges of median energy from $E$=0.7 to 3.1 TeV. The gradient varies with energy as $E^γ$, where $γ\approx-0.5$. At a median energy $\leq$1.0 TeV, this gradient corresponds to a dipole anisotropy of at least 1\%, or possibly a weaker anisotropy of higher order. This new type of observation opens the opportunity to study interplanetary magnetic structures using air shower arrays around the world, complementing existing in situ and remote measurements of plasma properties.

preprint2025arXiv

Coherent photoproduction of $ρ^0, ω$ and excited vector mesons in ultraperipheral PbPb collisions

The invariant-mass distribution for the coherent photoproduction of dipions in ultraperipheral PbPb collisions is measured using data, corresponding to an integrated luminosity of $224.6 \pm 9.6 μ$b$^{-1}$, collected by the LHCb experiment in 2018 at a nucleon-nucleon centre-of-mass energy $\sqrt{s_{\rm NN}}=5.02$ TeV. In the mass range from 400 to 1200 MeV, the results are consistent with previous experiments, with the spectrum dominated by the $ρ^0$ meson, which interferes with a nonresonant component, together with a smaller $ω$ meson contribution. In an extended mass range up to 2300 MeV, models previously used do not fit the data and a consistent description requires the introduction of two resonances at masses of $1350\pm20$ MeV and $1790\pm20$ MeV with widths of about 300 MeV. The cross-section for each meson is measured differentially in twelve bins of rapidity from 2.05 to 4.90. The $ρ^0$ cross-section increases with rapidity from about 400 to 600 mb and is measured with a typical precision of 8\%, while the cross-section times branching fraction for the $ω,ρ^\prime$ and $ρ^{\prime\prime}$, with the statistical precision of the data, do not have a pronounced rapidity dependence and are between 0.5 and 1.5 mb, with uncertainties up to 30\%. A large nuclear suppression is observed for the $ρ^0$ meson compared to expectations based on photoproduction on the proton that use the impulse approximation. Significant suppression is also observed compared to that predicted by elastic scattering described in the Glauber approach, or with the addition of inelastic scattering in a Gribov--Glauber model.

preprint2025arXiv

Observation of the $γ$-ray Emission from W43 with LHAASO

In this paper, we report the detection of the very-high-energy (VHE, $ 100{\rm\ GeV} < E < 100{\rm\ TeV} $) and ultra-high-energy (UHE, $E > 100\rm\ TeV$) $γ$-ray emissions from the direction of the young star-forming region W43, observed by the Large High Altitude Air Shower Observation (LHAASO). The extended $γ$-ray source was detected with a significance of ${\sim}16\,σ$ by KM2A and ${\sim}17\,σ$ by WCDA, respectively. The angular extension of this $γ$-ray source is about 0.5 degrees, corresponding to a physical size of about 50 pc. We discuss the origin of the $γ$-ray emission and possible cosmic ray acceleration in the W43 region using multi-wavelength data. Our findings suggest that W43 is likely another young star cluster capable of accelerating cosmic rays (CRs) to at least several hundred TeV.

preprint2025arXiv

Study of Ultra-High-Energy Gamma-Ray Source 1LHAASO J0056+6346u and Its Possible Origins

We report a dedicated study of the newly discovered extended UHE $γ$-ray source 1LHAASO J0056+6346u. Analyzing 979 days of LHAASO-WCDA data and 1389 days of LHAASO-KM2A data, we observed a significant excess of $γ$-ray events with both WCDA and KM2A. Assuming a point power-law source with a fixed spectral index, the significance maps reveal excesses of ${\sim}12.65\,σ$, ${\sim}22.18\,σ$, and ${\sim}10.24\,σ$ in the energy ranges of 1--25 TeV, 25--100 TeV, and $> 100$ TeV, respectively. We use a 3D likelihood algorithm to derive the morphological and spectral parameters, and the source is detected with significances of $12.65\,σ$ by WCDA and $25.27\,σ$ by KM2A. The best-fit positions derived from WCDA and KM2A data are (R.A. = $13.96^\circ\pm0.09^\circ$, Decl. = $63.92^\circ\pm0.05^\circ$) and (R.A. = $14.00^\circ\pm0.05^\circ$, Decl. = $63.79^\circ\pm0.02^\circ$), respectively. The angular size ($r_{39}$) of 1LHAASO J0056+6346u is $0.34^\circ\pm0.04^\circ$ at 1--25 TeV and $0.24^\circ\pm0.02^\circ$ at $> 25$ TeV. The differential flux of this UHE $γ$-ray source can be described by an exponential cutoff power-law function: $(2.67\pm0.25) \times 10^{-15} (E/20\,\text{TeV})^{-1.97\pm0.10} e^{-E/(55.1\pm7.2)\,\text{TeV}} \,\text{TeV}^{-1}\,\text{cm}^{-2}\,\text{s}^{-1}$. To explore potential sources of $γ$-ray emission, we investigated the gas distribution around 1LHAASO J0056+6346u. 1LHAASO J0056+6346u is likely to be a TeV PWN powered by an unknown pulsar, which would naturally explain both its spatial and spectral properties. Another explanation is that this UHE $γ$-ray source might be associated with gas content illuminated by a nearby CR accelerator, possibly the SNR candidate G124.0+1.4.

preprint2023arXiv

Longitudinal and transverse spin transfer to $Λ$ and $\overlineΛ$ hyperons in polarized $p$+$p$ collisions at $\sqrt{s} = 200$ GeV

The longitudinal and transverse spin transfers to $Λ$ ($\overlineΛ$) hyperons in polarized proton-proton collisions are expected to be sensitive to the helicity and transversity distributions, respectively, of (anti-)strange quarks in the proton, and to the corresponding polarized fragmentation functions. We report improved measurements of the longitudinal spin transfer coefficient, $D_{LL}$, and the transverse spin transfer coefficient, $D_{TT}$, to $Λ$ and $\overlineΛ$ in polarized proton-proton collisions at $\sqrt{s}$ = 200 GeV by the STAR experiment at RHIC. The data set includes longitudinally polarized proton-proton collisions with an integrated luminosity of 52 pb$^{-1}$, and transversely polarized proton-proton collisions with a similar integrated luminosity. Both data sets have about twice the statistics of previous results and cover a kinematic range of $|η_{Λ(\overlineΛ)}|$ $<$ 1.2 and transverse momentum $p_{T,{Λ(\overlineΛ)}}$ up to 8 GeV/$c$. We also report the first measurements of the hyperon spin transfer coefficients $D_{LL}$ and $D_{TT}$ as a function of the fractional jet momentum $z$ carried by the hyperon, which can provide more direct constraints on the polarized fragmentation functions.

preprint2022arXiv

AGN feedback duty cycle in Planck SZ selected clusters using Chandra observations

We present a systematic study of X-ray cavities using archival Chandra observations of nearby galaxy clusters selected by their Sunyaev-Zel'dovich (SZ) signature in the Planck survey, which provides a nearly unbiased mass-selected sample to explore the entire AGN feedback duty cycle. Based on X-ray image analysis, we report that 30 of the 164 clusters show X-ray cavities, which corresponds to a detection fraction of 18%. After correcting for spatial resolution to match the high-$z$ SPT-SZ sample, the detection fraction decreases to 9%, consistent with the high-z sample, hinting that the AGN feedback has not evolved across almost 8 Gyrs. Our finding agrees with the lack of evolution of cool-core clusters fraction. We calculate the cavity power, P_{\rm cav}, and find that most systems of our sample have enough AGN heating to offset the radiative losses of the intracluster medium.

preprint2022arXiv

Nonlinear co-generation of graphene plasmons for optoelectronic logic operations

Surface plasmons in graphene provide a compelling strategy for advanced photonic technologies thanks to their tight confinement, fast response and tunability. Recent advances in the field of all optical generation of graphene plasmons in planar waveguides offer a promising method for high speed signal processing in nanoscale integrated optoelectronic devices. Here, we use two counter propagating frequency combs with temporally synchronized pulses to demonstrate deterministic all optical generation and electrical control of multiple plasmon polaritons, excited via difference frequency generation (DFG). Electrical tuning of a hybrid graphene fibre device offers a precise control over the DFG phase matching, leading to tunable responses of the graphene plasmons at different frequencies across a broadband (0 - 50 THz) and provides a powerful tool for high speed logic operations. Our results offer insights for plasmonics on hybrid photonic devices based on layered materials and pave the way to high speed integrated optoelectronic computing circuits.

preprint2022arXiv

Pinch points and half-moons in dipolar-octupolar Nd$_2$Hf$_2$O$_7$

While it is established that the pinch point scattering pattern in spin ice arises from an emergent coulomb phase associated with magnetic moment that is divergence-free, more complex Hamiltonians can introduce a divergence-full part. If these two parts remain decoupled, they give rise to the co-existence of distinct features. Here we show that the moment in ${\rm Nd_2Hf_2O_7}$ forms a static long-range ordered ground state, a flat, gapped pinch point excitation and dispersive excitations. These results confirm recent theories which predict that the dispersive modes, which arise from the divergence-full moment, host a pinch point pattern of their own, observed experimentally as `half-moons'.

preprint2022arXiv

The cycle of metals in the infalling elliptical galaxy NGC 1404

Hot atmospheres pervading galaxy clusters, groups, and early-type galaxies are rich in metals, produced during epochs and diffused via processes that are still to be determined. While this enrichment has been routinely investigated in clusters, metals in lower mass systems are more challenging to probe with standard X-ray exposures and spectroscopy. In this paper, we focus on very deep XMM-Newton ($\sim$350 ks) observations of NGC 1404, a massive elliptical galaxy experiencing ram-pressure stripping of its hot atmosphere while infalling toward the centre of the Fornax cluster, with the aim to derive abundances through its hot gas extent. Importantly, we report the existence of a new fitting bias - the "double Fe bias" - leading to an underestimate of the Fe abundance when two thermal components cannot realistically model the complex temperature structure present in the outer atmosphere of the galaxy. Contrasting with the ''metal conundrum'' seen in clusters, the Fe and Mg masses of NGC 1404 are measured 1-2 orders of magnitude below what stars and supernovae could have reasonably produced and released. In addition, we note the remarkable Solar abundance ratios of the galaxy's halo, different from its stellar counterpart but similar to the chemical composition of the ICM of rich clusters. Completing the clusters regime, all these findings provide additional support toward a scenario of early enrichment, at play over two orders of magnitude in mass. A few peculiar and intriguing features, such as a possible double metal peak as well as an apparent ring of enhanced Si near the galaxy core, are also discussed.

preprint2022arXiv

The Infrared Database of Extragalactic Observables from Spitzer. -- II. The Database & The Diagnostic Power of Crystalline Silicate Features in Galaxy Spectra

We present the Infrared Database of Extragalactic Observables from Spitzer (IDEOS), a homogeneous, publicly available, database of 77 fitted mid-infrared observables in the 5.4-36um range, comprising measurements for 3335 galaxies observed in the low-resolution staring mode of the Infrared Spectrometer onboard the Spitzer Space Telescope. Among the included observables are PAH fluxes and their equivalent widths, the strength of the 9.8um silicate feature, emission line fluxes, solid-state features, rest frame continuum fluxes, synthetic photometry, and a mid-infrared spectral classification. The IDEOS spectra were selected from the Cornell Atlas of Spitzer-IRS Sources. To our surprise we have detected at a >95% confidence level crystalline silicates in the spectra of 786 IDEOS galaxies. The detections range from single band detections to detections of all fitted crystalline bands (16, 19, 23, 28 and 33um). We find the strength of the crystalline silicate bands to correlate with the amorphous silicate strength, and the change from an emission to an absorption feature to occur at higher obscuration as the wavelength of the crystalline silicate band is longer. These observed characteristics are consistent with an origin for the amorphous and crystalline silicate features in a centrally heated dust geometry, either an edge-on disk or a cocoon. We find the 23 and 33um crystalline silicate bands to be well-suited to classify the obscuration level of galactic nuclei, even in the presence of strong circumnuclear star formation. Based on our detection statistics, we conclude that crystalline silicates are a common component of the interstellar medium of galactic nuclei.

preprint2021arXiv

Distinct magnetic ground states of $R_2$ZnIrO$_6$ ($R$ = La and Nd) determined by neutron powder diffraction

Double perovskite iridates $A_2$ZnIrO$_6$ ($A$ = alkaline or lanthanide) show complex magnetic behaviors ranging from weak ferromagnetism to successive antiferromagnetic transitions. Here we report the static ($dc$) and dynamic ($ac$) magnetic susceptibility, and neutron powder diffraction measurements for $A$ = La and Nd compounds to elucidate the magnetic ground state. Below 10~K, the $A$ = La compound is best described as canted iridium moments in an antiferromagnet arrangement with a propagation vector \textbf{k} = 0 and a net ferromagnetic component along the $c$-axis. On the other hand, Nd$_2$ZnIrO$_6$ is described well as an antiferromagnet with a propagation vector \textbf{k} = (1/2~1/2~0) below $T_\mathrm{N} \sim$ 17 K. Scattering from both the Nd and Ir magnetic sublattices were required to describe the data and both were found to lie almost completely within the $ab$-plane. $Dc$ susceptibility revealed a bifurcation between the zero-field-cooled and field-cooled curves below $\sim$13 K in Nd$_2$ZnIrO$_6$. A glassy state was ruled out by $ac$ susceptibility but detailed magnetic isotherms revealed the opening of the loop below 13~K. These results suggest a delicate balance exists between the Dzyaloshinskii-Moriya, crystal field schemes, and $d$-$f$ interaction in this series of compounds.

preprint2020arXiv

Cation controlled wetting properties of vermiculite membranes and its potential for fouling resistant oil-water separation

The surface free energy is one of the most fundamental properties of solids, hence, manipulating the surface energy and thereby the wetting properties of solids, has tremendous potential for various physical, chemical, biological as well as industrial processes. Typically, this is achieved by either chemical modification or by controlling the hierarchical structures of surfaces. Here we report a phenomenon whereby the wetting properties of vermiculite laminates are controlled by the hydrated cations on the surface and in the interlamellar space. We find that by exploiting this mechanism, vermiculite laminates can be tuned from superhydrophillic to hydrophobic simply by exchanging the cations; hydrophilicity decreases with increasing cation hydration free energy, except for lithium. Lithium, which has a higher hydration free energy than potassium, is found to provide a superhydrophilic surface due to its anomalous hydrated structure at the vermiculite surface. Building on these findings, we demonstrate the potential application of superhydrophilic lithium exchanged vermiculite as a thin coating layer on microfiltration membranes to resist fouling, and thus, we address a major challenge for oil-water separation technology.

preprint2019arXiv

Neutron scattering study of commensurate magnetic ordering in single crystal CeSb$_2$

Temperature and field-dependent magnetization $M(H,T)$ measurements and neutron scattering study of a single crystal CeSb$_2$ are presented. Several anomalies in the magnetization curves have been confirmed at low magnetic field, i.e., 15.6 K, 12 K, and 9.8 K. These three transitions are all metamagnetic transitions (MMT), which shift to lower temperatures as the magnetic field increases. The anomaly at 15.6 K has been suggested as paramagnetic (PM) to ferromagnetic (FM) phase transition. The anomaly located at around 12 K is antiferromagnetic-like transition, and this turning point will clearly split into two when the magnetic field $H\geq0.2$ T. Neutron scattering study reveals that the low temperature ground state of CeSb$_2$ orders antiferromagnetically with commensurate propagation wave vectors $\textbf{k}=(-1,\pm1/6,0)$ and $\textbf{k}=(\pm1/6,-1,0)$, with Néel temperature $T_N\sim9.8$ K. This transition is of first-order, as shown in the hysteresis loop observed by the field cooled cooling (FCC) and field cooled warming (FCW) processes.

preprint2016arXiv

Chromospheric evaporation flows and density changes deduced from Hinode/EIS during an M1.6 flare

We analyzed high-cadence sit-and-stare observations acquired with the Hinode/EIS spectrometer and HXR measurements acquired with RHESSI during an M-class flare. During the flare impulsive phase, we observe no significant flows in the cooler Fe XIII line but strong upflows, up to 80-150 km/s, in the hotter Fe XVI line. The largest Doppler shifts observed in the Fe XVI line were co-temporal with the sharp intensity peak. The electron density obtained from a Fe XIII line pair ratio exhibited fast increase (within two minutes) from the pre-flare level of 5.01x10^(9) cm^(-3) to 3.16x10^(10) cm^(-3) during the flare peak. The nonthermal energy flux density deposited from the coronal acceleration site to the lower atmospheric layers during the flare peak was found to be 1.34x10^(10) erg/s/cm^(2) for a low-energy cut-off that was estimated to be 16 keV. During the decline flare phase, we found a secondary intensity and density peak of lower amplitude that was preceded by upflows of 15 km/s that were detected in both lines. The flare was also accompanied by a filament eruption that was partly captured by the EIS observations. We derived Doppler velocities of 250-300 km/s for the upflowing filament material.The spectroscopic results for the flare peak are consistent with the scenario of explosive chromospheric evaporation, although a comparatively low value of the nonthermal energy flux density was determined for this phase of the flare. This outcome is discussed in the context of recent hydrodynamic simulations. It provides observational evidence that the response of the atmospheric plasma strongly depends on the properties of the electron beams responsible for the heating, in particular the steepness of the energy distribution.

preprint2016arXiv

Magnetic polarization of Ir in underdoped, non-superconducting Eu(Fe$_{0.94}$Ir$_{0.06}$)$_{2}$As$_{2}$

Using polarized neutron diffraction and x-ray resonant magnetic scattering (XRMS) techniques, multiple phase transitions were revealed in an underdoped, non-superconducting Eu(Fe$_{1-x}$Ir$_{x}$)$_{2}$As$_{2}$ ($\mathit{x}$ = 0.06) single crystal. Compared with the parent compound EuFe$_{2}$As$_{2}$, the tetragonal-to-orthorhombic structural phase transition and the antiferromagnetic order of the Fe$^{2+}$ moments are significantly suppressed to $\mathit{T_{S}}$ = 111 (2) K and $\mathit{T_{N,Fe}}$= 85 (2) K by 6% Ir doping, respectively. In addition, the Eu$^{2+}$ spins order within the $\mathit{ab}$ plane in the A-type antiferromagnetic structure similar to the parent compound. However, the order temperature is evidently suppressed to $\mathit{T_{N,Eu}}$= 16.0 (5) K by Ir doping. Most strikingly, the XRMS measurements at the Ir $\mathit{L_{3}}$ edge demonstrates that the Ir 5$\mathit{d}$ states are also magnetically polarized, with the same propagation vector as the magnetic order of Fe. With $\mathit{T_{N,Ir}}$ = 12.0 (5) K, they feature a much lower onset temperature compared with $\mathit{T_{N,Fe}}$. Our observation suggests that the magnetism of the Eu sublattice has a considerable effect on the magnetic nature of the 5$\mathit{d}$ Ir dopant atoms and there exists a possible interplay between the localized Eu$^{2+}$ moments and the conduction $\mathit{d}$-electrons on the FeAs layers.

preprint2016arXiv

Molecular clouds and star formation toward the Galactic plane within 216.25$^{\rm o}\le$ l $\le$218.75$^{\rm o}$ and $-$0.75$^{\rm o} \le$ b $\le$1.25$^{\rm o}$

As a part of the Milky Way Imaging Scroll Painting (MWISP) survey, we performed a simultaneous 12CO(1-0), 13CO(1-0), C18O(1-0) mapping toward molecular clouds in a region encompassing 3.75 square degrees. We reveal three molecular clouds, the 15 km/s cloud, the 27 km/s cloud, and the 50 km/s cloud, in the surveyed region. The 50 km/s cloud is resolved with an angular resolution of around 1' for the first time. Investigating their morphology and velocity structures, we find that the 27 km/s cloud is likely affected by feedback from the stellar association Mon OB3 and the 50 km/s cloud is characterised by three large expanding molecular shells. The region is mapped in C18O (1-0) for the first time. We discover seven C18O clumps, which are likely to form massive stars, and 15 dust clumps based on the BGPS archival data. Using infrared color-color diagrams, we find 56 Class I and 107 Class II young stellar object (YSO) candidates. According to the distribution of YSO candidates, an overdensity is found around the HII region S287 and the intersection of two shells, indicative of triggering. The star formation efficiency and rate of the 27~km/s cloud are discussed. Comparing the observed values of the filament S287-main with the models of fragmentation, we suggest that turbulence controls the large scale of fragmentation in the filament while gravitational fragmentation plays an important role in the formation of YSOs on the small scale. We find that star-forming gas tends to have higher excitation temperature, higher 13CO opacity, and higher column density than non-star-forming gas, which is consistent with the point that star formation happens in denser gas and star-forming gas is heated by YSOs. Using the 1.1 mm dust emission to trace dense gas, we obtain a dense gas fraction of 2.7%-10.4% for the 27 km/s cloud.

preprint2016arXiv

SANS Study of Vortex Lattice Structural Transition in Optimally Doped (Ba$_{1-x}$K$_{x}$)Fe$_{2}$As$_{2}$

Small-angle neutron scattering on high quality single crystalline Ba$_{1-x}$K$_{x}$Fe$_{2}$As$_{2}$ reveals the transition from a low--field vortex solid phase with orientational order to a vortex polycrystal at high magnetic field. The vortex order-disorder transition is correlated with the second-peak feature in isothermal hysteresis loops, and is interpreted in terms of the generation of supplementary vortex solid dislocations. The sharp drop of the structure factor above the second peak field is explained by the dynamics of freezing of the vortex ensemble in the high field phase.

preprint2016arXiv

Spin-Wave and Electromagnon Dispersions in Multiferroic MnWO4 as Observed by Neutron Spectroscopy: Isotropic Heisenberg Exchange versus Anisotropic Dzyaloshinskii-Moriya Interaction

High resolution inelastic neutron scattering reveals that the elementary magnetic excitations in multiferroic MnWO4 consist of low energy dispersive electromagnons in addition to the well-known spin-wave excitations. The latter can well be modeled by a Heisenberg Hamiltonian with magnetic exchange coupling extending to the 12th nearest neighbor. They exhibit a spin-wave gap of 0.61(1) meV. Two electromagnon branches appear at lower energies of 0.07(1) meV and 0.45(1) meV at the zone center. They reflect the dynamic magnetoelectric coupling and persist in both, the collinear magnetic and paraelectric AF1 phase, and the spin spiral ferroelectric AF2 phase. These excitations are associated with the Dzyaloshinskii-Moriya exchange interaction, which is significant due to the rather large spin-orbit coupling.

preprint2015arXiv

Superconductivity in Ca-doped graphene

Graphene, a zero-gap semimetal, can be transformed into a metallic, semiconducting or insulating state by either physical or chemical modification. Superconductivity is conspicuously missing among these states despite considerable experimental efforts as well as many theoretical proposals. Here, we report superconductivity in calcium-decorated graphene achieved by intercalation of graphene laminates that consist of well separated and electronically decoupled graphene crystals. In contrast to intercalated graphite, we find that Ca is the only dopant that induces superconductivity in graphene laminates above 1.8 K among intercalants used in our experiments such as potassium, caesium and lithium. Ca-decorated graphene becomes superconducting at ~ 6 K and the transition temperature is found to be strongly dependent on the confinement of the Ca layer and the induced charge carrier concentration. In addition to the first evidence for superconducting graphene, our work shows a possibility of inducing and studying superconductivity in other 2D materials using their laminates.

preprint2015arXiv

The Confined X-class Flares of Solar Active Region 2192

The unusually large NOAA active region 2192, observed in October 2014, was outstanding in its productivity of major two-ribbon flares without coronal mass ejections. On a large scale, a predominantly north-south oriented magnetic system of arcade fields served as a strong, also lateral, confinement for a series of large two-ribbon flares originating from the core of the active region. The large initial separation of the flare ribbons, together with an almost absent growth in ribbon separation, suggests a confined reconnection site high up in the corona. Based on a detailed analysis of the confined X1.6 flare on October 22, we show how exceptional the flaring of this active region was. We provide evidence for repeated energy release, indicating that the same magnetic field structures were repeatedly involved in magnetic reconnection. We find that a large number of electrons was accelerated to non-thermal energies, revealing a steep power law spectrum, but that only a small fraction was accelerated to high energies. The total non-thermal energy in electrons derived (on the order of 10^25 J) is considerably higher than that in eruptive flares of class X1, and corresponds to about 10\% of the excess magnetic energy present in the active-region corona.

preprint2015arXiv

The Relation between Solar Eruption Topologies and Observed Flare Features I: Flare Ribbons

In this paper we present a topological magnetic field investigation of seven two-ribbon flares in sigmoidal active regions observed with Hinode, STEREO, and SDO. We first derive the 3D coronal magnetic field structure of all regions using marginally unstable 3D coronal magnetic field models created with the flux rope insertion method. The unstable models have been shown to be a good model of the flaring magnetic field configurations. Regions are selected based on their pre-flare configurations along with the appearance and observational coverage of flare ribbons, and the model is constrained using pre-flare features observed in extreme ultraviolet and X-ray passbands. We perform a topology analysis of the models by computing the squashing factor, Q, in order to determine the locations of prominent quasi-separatrix layers (QSLs). QSLs from these maps are compared to flare ribbons at their full extents. We show that in all cases the straight segments of the two J-shaped ribbons are matched very well by the flux-rope-related QSLs, and the matches to the hooked segments are less consistent but still good for most cases. In addition, we show that these QSLs overlay ridges in the electric current density maps. This study is the largest sample of regions with QSLs derived from 3D coronal magnetic field models, and it shows that the magnetofrictional modeling technique that we employ gives a very good representation of flaring regions, with the power to predict flare ribbon locations in the event of a flare following the time of the model.

preprint2014arXiv

Coexistence of Superconductivity and Ferromagnetism in P-doped EuFe2As2

The magnetic structure of the Eu2+ moments in the superconducting EuFe2(As1-xPx)2 sample with x = 0.15 has been determined using element specific x-ray resonant magnetic scattering. Combining magnetic, thermodynamic and scattering measurements, we conclude that the long range ferromagnetic order of the Eu2+ moments aligned primarily along the c axis coexists with the bulk superconductivity at zero field. At an applied magnetic field >= 0.6 T, superconductivity still coexists with the ferromagnetic Eu2+ moments which are polarized along the field direction. We propose a spontaneous vortex state for the coexistence of superconductivity and ferromagnetism in EuFe2(As0.85P0.15)2.

preprint2014arXiv

Impermeable Barrier Films and Protective Coatings Based on Reduced Graphene Oxide

Barrier films preventing permeation of gases and moistures are important for many industries ranging from food to medical and from chemical to electronic. From this perspective, graphene has recently attracted particular interest because its defect free monolayers are impermeable to all gases and liquids. However, it has proved challenging to develop large-area defectless graphene films suitable for industrial use. Here we report barrier properties of multilayer graphitic films made by chemical reduction of easily and cheaply produced graphene oxide laminates. They are found to provide a practically perfect barrier that blocks all gases, liquids and aggressive chemicals including, for example, hydrofluoric acid. In particular, if graphene oxide laminates are reduced in hydroiodic acid, no permeation of hydrogen and water could be detected for films as thin as 30 nm, which remain optically transparent. The films thicker than 100 nm become completely impermeable. The exceptional barrier properties are attributed to a high degree of graphitization of the laminates and little structural damage during reduction. This work indicates a close prospect of thin protective coatings with stability and inertness similar to that of graphene and bulk graphite, which can be interesting for numerous applications.

preprint2014arXiv

Incipient Ferromagnetism in Tb2Ge2O7: Application of Chemical Pressure to the Enigmatic Spin Liquid, Tb2Ti2O7

The origin of the spin liquid state in Tb$_2$Ti$_2$O$_7$ has challenged experimentalists and theorists alike for nearly 20 years. To improve our understanding of the exotic magnetism in Tb$_2$Ti$_2$O$_7$, we have synthesized a chemical pressure analog, Tb$_2$Ge$_2$O$_7$. Germanium substitution results in a lattice contraction and enhanced exchange interactions. We have characterized the magnetic ground state of Tb$_2$Ge$_2$O$_7$ with specific heat, ac and dc magnetic susceptibility, and polarized neutron scattering measurements. Akin to Tb$_2$Ti$_2$O$_7$, there is no long-range order in Tb$_2$Ge$_2$O$_7$ down to 20 mK. The Weiss temperature of $-$19.2(1) K, which is more negative than that of Tb$_2$Ti$_2$O$_7$, supports the picture of stronger antiferromagnetic exchange. Polarized neutron scattering of Tb$_2$Ge$_2$O$_7$ reveals that at 3.5 K liquid-like correlations dominate in this system. However, below 1 K, the liquid-like correlations give way to intense short-range ferromagnetic correlations with a length scale related to the Tb-Tb nearest neighbor distance. Despite stronger antiferromagnetic exchange, the ground state of Tb$_2$Ge$_2$O$_7$ has ferromagnetic character, in stark contrast to the pressure-induced antiferromagnetic order observed in Tb$_2$Ti$_2$O$_7$.

preprint2014arXiv

Magnetic and structural transitions in La$_{0.4}$Na$_{0.6}$Fe$_2$As$_2$ single crystals

La$_{0.4}$Na$_{0.6}$Fe$_2$As$_2$ single crystals have been grown out of an NaAs flux in an alumina crucible and characterized by measuring magnetic susceptibility, electrical resistivity, specific heat, as well as single crystal x-ray and neutron diffraction. La$_{0.4}$Na$_{0.6}$Fe$_2$As$_2$ single crystals show a structural phase transition from a high temperature tetragonal phase to a low-temperature orthorhombic phase at T$_s$\,=\,125\,K. This structural transition is accompanied by an anomaly in the temperature dependence of electrical resistivity, anisotropic magnetic susceptibility, and specific heat. Concomitant with the structural phase transition, the Fe moments order along the \emph{a} direction with an ordered moment of 0.7(1)\,$μ_{\textup{B}}$ at \emph{T}\,=\,5 K. The low temperature stripe antiferromagnetic structure is the same as that in other \emph{A}Fe$_{2}$As$_{2}$ (\emph{A}\,=\,Ca, Sr, Ba) compounds. La$_{0.5-x}$Na$_{0.5+x}$Fe$_2$As$_2$ provides a new material platform for the study of iron-based superconductors where the electron-hole asymmetry could be studied by simply varying La/Na ratio.

preprint2014arXiv

Magnetic structure of the Eu2+ moments in superconducting EuFe2(As1-xPx)2 with x = 0.19

The magnetic structure of the Eu2+ moments in the superconducting EuFe2(As1-xPx)2 sample with x = 0.19 has been determined using neutron scattering. We conclude that the Eu2+ moments are aligned along the c direction below T_C = 19.0(1) K with an ordered moment of 6.6(2) mu_B in the superconducting state. An impurity phase similar to the underdoped phase exists within the bulk sample which orders antiferromagnetically below T_N = 17.0(2) K. We found no indication of iron magnetic order, nor any incommensurate magnetic order of the Eu2+ moments in the sample.

preprint2014arXiv

Precise and ultrafast molecular sieving through graphene oxide membranes

There has been intense interest in filtration and separation properties of graphene-based materials that can have well-defined nanometer pores and exhibit low frictional water flow inside them. Here we investigate molecular permeation through graphene oxide laminates. They are vacuum-tight in the dry state but, if immersed in water, act as molecular sieves blocking all solutes with hydrated radii larger than 4.5A. Smaller ions permeate through the membranes with little impedance, many orders of magnitude faster than the diffusion mechanism can account for. We explain this behavior by a network of nanocapillaries that open up in the hydrated state and accept only species that fit in. The ultrafast separation of small salts is attributed to an 'ion sponge' effect that results in highly concentrated salt solutions inside graphene capillaries.

preprint2013arXiv

Magnetic anisotropic energy gap and low energy spin wave excitation in antiferromagnetic block phase of K2Fe4Se5

Neutron scattering experiments were performed to investigate magnetic order and magnetic excitations in ternary iron chalcogenide K2Fe4Se5. The formation of a superlattice structure below 580 K together with the decoupling between the Fe-vacancy order-disorder transition and the antiferromagnetic order transition appears to be a common feature in the A2Fe4Se5 family. The study of spin dynamics of K2Fe4Se5 reveals two distinct energy gaps at the magnetic Brillouin zone center, which indicates the presence of magnetic anisotropy and the decrease of local symmetry due to electronic and orbital anisotropy. The low-energy spin wave excitations of K2Fe4Se5 can be properly described by linear spin wave theory within a Heisenberg model. Compared to iron pnictides, K2Fe4Se5 exhibits a more two-dimensional magnetism as characterized by large differences not only between out-of-plane and in-plane spin wave velocities, but also between out-of-plane and in-plane exchange interactions.

preprint2013arXiv

Magnetic structure of superconducting Eu(Fe0.82Co0.18)2As2 as revealed by single-crystal neutron diffraction

The magnetic structure of superconducting Eu(Fe0.82Co0.18)2As2 is unambiguously determined by single-crystal neutron diffraction. A long-range ferromagnetic order of the Eu2+ moments along the c-direction is revealed below the magnetic phase transition temperature Tc = 17 K. In addition, the antiferromagnetism of the Fe2+ moments still survives and the tetragonal-to-orthorhombic structural phase transition is also observed, although the transition temperatures of the Fe-spin density wave (SDW) order and the structural phase transition are significantly suppressed to Tn = 70 K and Ts = 90 K, respectively, compared to the parent compound EuFe2As2.We present the microscopic evidences for the coexistence of the Eu-ferromagnetism (FM) and the Fe-SDW in the superconducting crystal. The superconductivity (SC) competes with the Fe-SDW in Eu(Fe0.82Co0.18)2As2.Moreover, the comparison between Eu(Fe1-xCox)2As2 and Ba(Fe1-xCox)2As2 indicates a considerable influence of the rare-earth element Eu on the magnetism of the Fe sublattice.

preprint2013arXiv

Magnetization distribution and orbital moment in the non-Superconducting Chalcogenide Compound K0.8Fe1.6Se2

We have used polarized and unpolarized neutron diffraction to determine the spatial distribution of the magnetization density induced by a magnetic field of 9 T in the tetragonal phase of K0.8Fe1.6Se2. The maximum entropy reconstruction shows clearly that most of the magnetization is confined to the region around the iron atoms whereas there is no significant magnetization associated with either Se or K atoms. The distribution of magnetization around the Fe atom is slightly nonspherical with a shape which is extended along the [0 0 1] direction in the projection. Multipolar refinement results show that the electrons which give rise to the paramagnetic susceptibility are confined to the Fe atoms and their distribution suggests that they occupy 3d t2g-type orbitals with around 66% in those of xz/yz symmetry. Detail modeling of the magnetic form factor indicates the presence of an orbital moment to the total paramagnetic moment of Fe2+

preprint2013arXiv

Spin-Lattice Coupling in K0.8Fe1.6Se2 and KFe2Se2: Inelastic Neutron Scattering and ab-initio Phonon Calculations

We report measurements of the temperature dependence of phonon densities of states in K0.8Fe1.6Se2 using inelastic neutron scattering technique. While cooling down to 150 K, a phonon peak splitting around 25 meV is observed and a new peak appears at 31 meV. The measurements support the recent Raman and infra-red measurements indicating a lowering of symmetry of K0.8Fe1.6Se2 upon cooling below 250 K. Ab-initio phonon calculations have been carried out for K0.8Fe1.6Se2 and KFe2Se2. The comparison of the phonon spectra as obtained from the magnetic as well as non magnetic calculations show pronounced differences. We show that in the two calculations the energy range of the vibrational contribution from both Fe and Se are quite different. We conclude that Fe magnetism is correlated to the phonon dynamics and it plays an important role in stabilizing the structure of K0.8Fe1.6Se2 as well as that of KFe2Se2. The calculations highlight the presence of low energy librational modes in K0.8Fe1.6Se2 as compared to KFe2Se2.

preprint2012arXiv

Dynamics aspect of subbarrier fusion reaction in light heavy ion systems

Subbarrier fusion of the 7Li + 12C reaction is studied using an antisymmetrized molecular dynamics model (AMD) with an after burner, GEMINI. In AMD, 7Li shows an α+ t structure at its ground state and it is significantly deformed. Simulations are made near the Coulomb barrier energies, i.e., E_{cm} = 3 - 8 MeV. The total fusion cross section of the AMD + GEMINI calculations as a function of incident energy is compared to the experimental results and both are in good agreement at E_{cm} > 3 MeV. The cross section for the different residue channels of the AMD + GEMINI at E_{cm} = 5 MeV is also compared to the experimental results.

preprint2012arXiv

Low-Altitude Reconnection Inflow-Outflow Observations during a 2010 November 3 Solar Eruption

For a solar flare occurring on 2010 November 3, we present observations using several SDO/AIA extreme-ultraviolet (EUV) passbands of an erupting flux rope followed by inflows sweeping into a current sheet region. The inflows are soon followed by outflows appearing to originate from near the termination point of the inflowing motion - an observation in line with standard magnetic reconnection models. We measure average inflow plane-of-sky speeds to range from ~150-690 km/s with the initial, high-temperature inflows being the fastest. Using the inflow speeds and a range of Alfven speeds, we estimate the Alfvenic Mach number which appears to decrease with time. We also provide inflow and outflow times with respect to RHESSI count rates and find that the fast, high-temperature inflows occur simultaneously with a peak in the RHESSI thermal lightcurve. Five candidate inflow-outflow pairs are identified with no more than a minute delay between detections. The inflow speeds of these pairs are measured to be 10^2 km/s with outflow speeds ranging from 10^2-10^3 km/s - indicating acceleration during the reconnection process. The fastest of these outflows are in the form of apparently traveling density enhancements along the legs of the loops rather than the loop apexes themselves. These flows could either be accelerated plasma, shocks, or waves prompted by reconnection. The measurements presented here show an order of magnitude difference between the retraction speeds of the loops and the speed of the density enhancements within the loops - presumably exiting the reconnection site.

preprint2011arXiv

Antiferromagnetic Order in MnO Spherical Nanoparticles

We have performed unpolarized and polarized neutron diffraction experiments on monodisperse 8 nm and 13 nm antiferromagnetic MnO nanoparticles. For the 8 nm sample, the antiferromagnetic transition temperature $T_N$ (114 K) is suppressed compared to the bulk material (119 K) while for the 13 nm sample $T_N$ (120 K) is comparable to the bulk. The neutron diffraction data of the nanoparticles is well described using the bulk MnO magnetic structure but with a substantially reduced average magnetic moment of 4.2$\pm$0.3 $μ_B$/Mn for the 8 nm sample and 3.9$\pm$0.2 $μ_B$/Mn for the 13 nm sample. An analysis of the polarized neutron data on both samples shows that in an individual MnO nanoparticle about 80$%$ of Mn ions order. These results can be explained by a structure in which the monodisperse nanoparticles studied here have a core that behaves similar to the bulk with a surface layer which does not contribute significantly to the magnetic order.

preprint2011arXiv

Competing Ferri- and Antiferromagnetic Phases in Geometrically Frustrated LuFe2O4

We present a detailed study of magnetism in LuFe2O4, combining magnetization measurements with neutron and soft x-ray diffraction. The magnetic phase diagram in the vicinity of T_N involves a metamagnetic transition separating an antiferro- and a ferrimagnetic phase. For both phases the spin structure is refined by neutron diffraction. Observed diffuse magnetic scattering far above T_N is explained in terms of near degeneracy of the magnetic phases.

preprint2011arXiv

Pressure dependence of the low- temperature crystal structure and phase transition behaviour of CaFeAsF and SrFeAsF: A synchrotron x-ray diffraction study

We report systematic investigation of high pressure crystal structures and structural phase transition upto 46 GPa in CaFeAsF and 40 GPa in SrFeAsF at 40 K using powder synchrotron x-ray diffraction experiments and Rietveld analysis of the diffraction data. We find that CaFeAsF undergoes orthorhombic to monoclinic phase transition at Pc = 13.7 GPa while increasing pressure. SrFeAsF exhibits coexistence of orthorhombic and monoclinic phases over a large pressure range from 9 to 39 GPa. The coexistence of the two phases indicates that the transition is of first order in nature. Unlike in the 122 compounds (BaFe2As2 & CaFe2As2) we do not find any collapse tetragonal transition. The transition to a lower symmetry phase (orthorhombic to monoclinic) in 1111 compounds under pressure is in contrast with the transition to a high symmetry phase (orthorhombic to tetragonal) in 122 type compounds. On heating from 40 K at high pressure, CaFeAsF undergoes monoclinic to tetragonal phase transition around 25 GPa and 200 K. Further, it does not show any post-tetragonal phase transition and remains in the tetragonal phase upto 25 GPa at 300 K. The dPc/dT is found to be positive for the CaFeAsF & CaFe2As2, however the same was not found in case of BaFe2As2. We discuss observations of structural evolution in the context of superconductivity in these and other Fe-based compounds. It appears that the closeness of the Fe-As-Fe bond angle to its ideal tetrahedral value of 109.470 might be associated with occurrence of superconductivity at low temperature.

preprint2011arXiv

Spatial Inhomogeneity in RFeAsO1-xFx (R =Pr, Nd) Determined from Rare Earth Crystal Field Excitations

We report inelastic neutron scattering measurements of crystal field transitions in PrFeAsO, PrFeAsO0.87F0.13, and NdFeAsO0.85F0.15. Doping with fluorine produces additional crystal field excitations, providing evidence that there are two distinct charge environments around the rare earth ions, with probabilities that are consistent with a random distribution of dopants on the oxygen sites. The 4f electrons of the Pr3+ and Nd3+ ions have non-magnetic and magnetic ground states, respectively, indicating that the enhancement of Tc compared to LaFeAsO1-xFx is not due to rare earth magnetism.

preprint2011arXiv

Strong coupling of Sm and Fe magnetism in SmFeAsO as revealed by magnetic x-ray scattering

The magnetic structures adopted by the Fe and Sm sublattices in SmFeAsO have been investigated using element specific x-ray resonant and non-resonant magnetic scattering techniques. Between 110 and 5 K, the Sm and Fe moments are aligned along the c and a directions, respectively according to the same magnetic representation $Γ_{5}$ and the same propagation vector (1, 0, 0.5). Below 5 K, magnetic order of both sublattices change to a different magnetic structure and the Sm moments reorder in a magnetic unit cell equal to the chemical unit cell. Modeling of the temperature dependence for the Sm sublattice as well as a change in the magnetic structure below 5 K provide a clear evidence of a surprisingly strong coupling between the two sublattices, and indicate the need to include anisotropic exchange interactions in models of SmFeAsO and related compounds.

preprint2010arXiv

Diffuse Neutron Scattering Study of Magnetic Correlations in half-doped La0.5Ca0.5-xSrxMnO3 (x = 0.1, 0.3 and 0.4) Manganites

The short range ordered magnetic correlations have been studied in half doped La0.5Ca0.5-xSrxMnO3 (x = 0.1, 0.3 and 0.4) compounds by polarized neutron scattering technique. On doping Sr2+ for Ca2+ ion, these compounds with x = 0.1, 0.3, and 0.4 exhibit CE-type, mixture of CE-type and A-type, and A-type antiferromagnetic ordering, respectively. Magnetic diffuse scattering is observed in all the compounds above and below their respective magnetic ordering temperatures and is attributed to magnetic polarons. The correlations are primarily ferromagnetic in nature above T\_N, although a small antiferromagnetic contribution is also evident. Additionally, in samples x = 0.1 and 0.3 with CE-type antiferromagnetic ordering, superlattice diffuse reflections are observed indicating correlations between magnetic polarons. On lowering temperature below T\_N the diffuse scattering corresponding to ferromagnetic correlations is suppressed and the long range ordered antiferromagnetic state is established. However, the short range ordered correlations indicated by enhanced spin flip scattering at low Q coexist with long range ordered state down to 3K. In x = 0.4 sample with A-type antiferromagnetic ordering, superlattice diffuse reflections are absent. Additionally, in comparison to x = 0.1 and 0.3 sample, the enhanced spin flip scattering at low Q is reduced at 310K, and as temperature is reduced below 200K, it becomes negligibly low. The variation of radial correlation function, g(r) with temperature indicates rapid suppression of ferromagnetic correlations at the first nearest neighbor on approaching TN. Sample x = 0.4 exhibits growth of ferromagnetic phase at intermediate temperatures (~ 200K). This has been further explored using SANS and neutron depolarization techniques.

preprint2010arXiv

Field induced spin reorientation and giant spin-lattice coupling in EuFe2As2

We have studied a EuFe2As2 single crystal by neutron diffraction under magnetic fields up to 3.5 T and temperatures down to 2 K. A field induced spin reorientation is observed in the presence of a magnetic field along both the a and c axes, respectively. Above critical field, the ground state antiferromagnetic configuration of Eu$^{2+}$ moments transforms into a ferromagnetic structure with moments along the applied field direction. The magnetic phase diagram for Eu magnetic sublattice in EuFe2As2 is presented. A considerable strain ($\sim$0.9%) is induced by the magnetic field, caused by the realignment of the twinning structure. Furthermore, the realignment of the twinning structure is found to be reversible with the rebound of magnetic field, which suggested the existence of magnetic shape-memory effect. The Eu moment ordering exhibits close relationship with the twinning structure. We argue that the Zeeman energy in combined with magnetic anisotropy energy is responsible for the observed spin-lattice coupling.

preprint2010arXiv

Helical spin-waves, magnetic order, and fluctuations in the langasite compound Ba3NbFe3Si2O14

We have investigated the spin fluctuations in the langasite compound Ba3NbFe3Si2O14 in both the ordered state and as a function of temperature. The low temperature magnetic structure is defined by a spiral phase characterized by magnetic Bragg peaks at q=(0,0,tau ~ 1/7) onset at TN=27 K as previously reported by Marty et al. The nature of the fluctuations and temperature dependence of the order parameter is consistent with a classical second order phase transition for a two dimensional triangular antiferromagnet. We will show that the physical properties and energy scales including the ordering wavevector, Curie-Weiss temperature, and the spin-waves can be explained through the use of only symmetric exchange constants without the need for the Dzyaloshinskii-Moriya interaction. This is accomplished through a set of ``helical" exchange pathways along the c direction imposed by the chiral crystal structure and naturally explains the magnetic diffuse scattering which displays a strong vector chirality up to high temperatures well above the ordering temperature. This illustrates a strong coupling between magnetic and crystalline chirality in this compound.

preprint2010arXiv

Magnetic correlations in spin ice Ho2-xYxTi2O7 as revealed by neutron polarization analysis

We present single-crystal neutron-diffuse-scattering measurements with polarization analysis on the spin ice Ho2-xYxTi2O7 (x = 0, 0.3, and 1). At 2 K, the spin flip scattering patterns are typical of a nearest-neighbor spin ice and the effects of Y doping are seen in the width of the pinch points. At low-T, the spin flip patterns are characteristic of a dipolar spin ice and are unaffected by Y doping. In the non-spin flip channel, we observe the signature of strong antiferromagnetic correlations at the same T as the dipolar spin ice correlations appear in the spin flip channel. These experimental observations have been well reproduced by our Monte Carlo simulations.

preprint2010arXiv

Magnetic structure and orbital ordering in KCrF3 perovskite

KCrF3 represents another prototypical orbital-ordered perovskite, where Cr2+ possesses the same electronic configuration of 3d4 as that of strongly Jahn-Teller distorted Mn3+ in many CMR manganites. The crystal and magnetic structures of KCrF3 compound are investigated by using polarized and unpolarized neutron powder diffraction methods. The results show that the KCrF3 compound crystallizes in tetragonal structure at room temperature and undergoes a monoclinic distortion with the decrease in temperature. The distortion of the crystal structure indicates the presence of cooperative Jahn-Teller distortion which is driven by orbital ordering. With decreasing temperature, four magnetic phase transitions are observed at 79.5, 45.8, 9.5, and 3.2 K, which suggests a rich magnetic phase diagram. Below T_N = 79.5 K, the Cr2+ moment orders in an incommensurate antiferromagnetic arrangement, which can be defined by the magnetic propagation vector (1/2$\pm\,$$δ\,$, 1/2$\pm\,$$δ\,$, 0). The incommensurate-commensurate magnetic transition occurs at 45.8 K and the magnetic propagation vector locks into (1/2, 1/2, 0) with the Cr moment of 3.34(5) $μ_B \,$ aligned ferromagnetically in (220) plane, but antiferromagnetically along [110] direction. Below 9.5 K, the canted antiferromagnetic ordering and weak ferromagnetism arise from the collinear antiferromagnetic structure, while the Dzyaloshinskii-Moriya interaction and tilted character of the single-ion anisotropy might give rise to the complex magnetic behaviors below 9.5 K.

preprint2010arXiv

Magnetisation distribution in the tetragonal phase of BaFe2As2

We have determined the spatial distribution of the magnetisation induced by a field of 9 T in the tetragonal phase of BaFe2As2 using polarised neutron diffraction. Magnetic structure factors derived from the polarisation dependence of the intensities of Bragg reflections were used to make a maximum entropy reconstruction of the distribution projected on the 110 plane. The reconstruction shows clearly that the magnetisation is confined to the region around the iron atoms and that there is no significant magnetisation associated with either the As or Ba atoms. The distribution of magnetisation around the Fe atom is significantly non-spherical with a shape which is extended in the <111> directions in the projection. These results show that the electrons which give rise to the paramagnetic susceptibility are confined to the Fe atoms their distribution suggests that they occupy 3d t_2g type orbitals with about 60% in those of xy symmetry.

preprint2010arXiv

Phonon dynamics in SrFe2As2 and SrFeAsF parent pnictide compounds: How different/similar when compared to Ca- and Ba-compounds?

We report detailed measurements of temperature dependence and ab initio lattice dynamical calculations of the phonon density-of-states of two Sr based iron pnictide parent compounds; SrFe2As2 and SrFeAsF. Although both compounds undergo a magnetic as well as a tetragonal-to-orthorhombic structural phase transitions, no drastic change of their phonon spectra is found in the temperature range from 2 to 300 K.. In both the systems the Fe-As stretching modes at 32 meV react to the temperature change and shift towards higher energies with decreasing temperature. Furthermore, we found that in the antiferromagnetic phase of SrFe2As2 at 180 K a peak at 25 meV, due to Fe vibrations, moves slightly to higher energies and become narrower on cooling from 300 K. However we did not find any evidence of softening of low energy phonon modes in both SrFe2As2 and SrFeAsF on decreasing the temperature from 300 K in their paramagnetic phase. This is contrary to our earlier observation of phonon softening in CaFe2As2, Ca0.6Na0.4Fe2As2 and Sr0.6K0.4Fe2As2. For the sake of comparison, the previously studied Ca- and Ba-cases are also presented and the role of the magnetic ordering on the calculated phonon spectra in these two families of pnictides is reported.

preprint2010arXiv

Soft X-ray resonant scattering study of single-crystal LaSr$_2$Mn$_2$O$_7$

Soft X-ray resonant scattering studies at the Mn $L_{\texttt{II, III}}$- and the La $M_{\texttt{IV, V}}$- edges of single-crystal LaSr$_2$Mn$_2$O$_7$ are reported. At low temperatures, below $T_\texttt{N} \approx 160$ K, energy scans with a fixed momentum transfer at the \emph{A}-type antiferromagnetic (0 0 1) reflection around the Mn $L_{\texttt{II, III}}$-edges with incident linear $σ$ and $π$ polarizations show strong resonant enhancements. The splitting of the energy spectra around the Mn $L_{\texttt{II, III}}$-edges may indicate the presence of a mixed valence state, e.g., Mn$^{3+}$/Mn$^{4+}$. The relative intensities of the resonance and the clear shoulder-feature as well as the strong incident $σ$ and $π$ polarization dependences strongly indicate its complex electronic origin. Unexpected enhancement of the charge Bragg (0 0 2) reflection at the La $M_{\texttt{IV, V}}$-edges with $σ$ polarization has been observed up to 300 K, with an anomaly appearing around the orbital-ordering transition temperature, $T_{\texttt{OO}} \approx 220$ K, suggesting a strong coupling (competition) between them.

preprint2010arXiv

Synchrotron X-ray Diffraction Study of BaFe2As2 and CaFe2As2 at High Pressures up to 56 GPa: Ambient and Low-Temperatures Down to 33 K

We report high pressure powder synchrotron x-ray diffraction studies on MFe2As2 (M=Ba, Ca) over a range of temperatures and pressures up to about 56 GPa using a membrane diamond anvil cell. A phase transition to a collapsed tetragonal phase is observed in both compounds upon compression. However, at 300 (33) K in the Ba-compound the transition occurs at 26 (29) GPa, which is a much higher pressure than 1.7 (0.3) GPa at 300 (40) K in the Ca-compound, due to its larger volume. It is important to note that the transition in both compounds occurs when they are compressed to almost the same value of the unit cell volume and attain similar ct/at ratios. We also show that the FeAs4 tetrahedra are much less compressible and more distorted in the collapsed tetragonal phase than their nearly regular shape in the ambient pressure phase. We present a detailed analysis of the pressure dependence of the structures as well as equation of states in these important BaFe2As2 and CaFe2As2 compounds.

preprint2009arXiv

Effects of magnetic doping and temperature dependence on phonon dynamics in CaFe\_{1-x}Co\_{x}AsF compounds (x = 0, 0.06, 0.12)

We report detailed measurements of composition as well as temperature dependence of the phonon density-of-states in a new series of FeAs compounds with composition CaFe1\_{1-x}Co\_{x}AsF (x = 0, 0.06, 0.12). The composition as well as temperature dependence of phonon spectra for CaFe\_{1-x}Co\_{x}AsF (x = 0, 0.06, 0.12) compounds have been measured using time of flight IN4C and IN6 spectrometers at ILL, France. The comparison of phonon spectra at 300 K in these compounds shows that acoustic phonon modes up to 12 meV harden in the doped compounds in comparison to the parent CaFeAsF. While intermediate energy phonon modes from 15 meV to 25 meV are also found to shift towards high energies only in the 12 % Co doped CaFeAsF compound. The experimental results for CaFe\_{1-x}Co\_{x}AsF (x = 0, 0.06, 0.12) are quite different from our previous phonon studies on parent and superconducting MFe2As2 (M=Ba, Ca, Sr) where low-energy acoustic phonon modes do not react with doping, while the phonon spectra in the intermediate range from 15 to 25 K are found to soften in these compounds. We argue that stronger spin phonon interaction play an important role for the emergence of superconductivity in these compounds. The lattice dynamics of CaFe\_{1-x}Co\_{x}AsF (x = 0, 0.06, 0.12) compounds is also investigated using the ab-initio as well as shell model phonon calculations. We show that the nature of the interaction between the Ca and the Fe-As layers in CaFeAsF compounds is quite different compared with our previous studies on CaFe2As2.

preprint2009arXiv

Neutron diffraction study on phase transition and thermal expansion of SrFeAsF

The magnetic ordering and crystal structure of iron pnictide SrFeAsF was investigated by using neutron powder diffraction method. With decreasing temperature, the tetragonal to orthorhombic phase transition is found at 180 K, while the paramagnetic to antiferromagnetic phase transition set in at 133 K. Similar to the parent compound of other iron pnictide system, the striped Fe magnetism is confirmed in antiferromagnetic phase and the Fe moment of 0.58(6) uB aligned along long a axis. The thermal expansion of orthorhombic phase of SrFeAsF is also investigated. Based on the Grueneisen approximation and Debye approximation for internal energy, the volume of SrFeAsF can be well fitted with Debye temperature of 347(5) K. The experimental atomic displacement parameters for different crystallographic sites in SrFeAsF are analyzed with Debye model. The results suggested that the expansion of FeAs layers plays an important role in determining the thermal expansion coefficient.

preprint2009arXiv

Phonon spectra in CaFe2As2 and Ca0.6Na0.4Fe2As2: Measurement of the pressure and temperature dependence and comparison with ab-initio and shell model calculations

We report the pressure and temperature dependence of the phonon density-of-states in superconducting Ca0.6Na0.4Fe2As2 (Tc=21 K) and the parent compound CaFe2As2, using inelastic neutron scattering. We observe no significant change in the phonon spectrum for Ca0.6Na0.4Fe2As2 at 295 K up to pressures of 5 kbar. The phonon spectrum for CaFe2As2 shows softening of the low-energy modes by about 1 meV when decreasing the temperature from 300 K to 180 K. There is no appreciable change in the phonon density of states across the structural and anti-ferromagnetic phase transition at 172 K. These results, combined with our earlier temperature dependent phonon density of states measurements for Ca0.6Na0.4Fe2As2, indicate that the softening of low-energy phonon modes in these compounds may be due to the interaction of phonons with electron or short-range spin fluctuations in the normal state of the superconducting compound as well as in the parent compound. The phonon spectra are analyzed with ab-initio and empirical model calculations giving partial densities of states and dispersion relations.

preprint2009arXiv

Pressure dependence of phonon modes across the tetragonal to collapsed tetragonal phase transition in CaFe2As2

The pressure dependence of a large number of phonon modes in CaFe2As2 with energies covering the full range of the phonon spectrum has been studied using inelastic x-ray and neutron scattering. The observed phonon frequency changes are in general rather small despite the sizable changes of the lattice parameters at the phase transition. This indicates that the bonding properties are not profoundly altered by the phase transition. The transverse acoustic phonons propagating along the c-direction are an exception because they stiffen very significantly in response to the large contraction of the c-axis. The lattice parameters are found to change significantly as a function of pressure before, during and after the first-order phase transition. However, the frequencies change nearly uniformly with the change in the lattice parameters due to pressure, with no regard specifically to the first-order phase transition. Density functional theory describes the frequencies in both the zero pressure and in the collapsed phase in a satisfactory way if based on the respective crystal structures.