Search

With sizable volatile envelopes but smaller radii than the solar system ice giants, sub-Neptunes have been revealed as one of the most common types of planet in the galaxy. While the spectroscopic characterization of larger sub-Neptunes (2.5-4R$_\oplus$) has revealed hydrogen-dominated atmospheres, smaller sub-Neptunes (1.6--2.5R$_\oplus$) could either host thin, rapidly evaporating hydrogen-rich atmospheres or be stable metal-rich "water worlds" with high mean molecular weight atmospheres and a fundamentally different formation and evolutionary history. Here, we present the 0.6--2.8$\mu$m JWST NIRISS/SOSS transmission spectrum of GJ 9827 d, the smallest (1.98 R$_\oplus$) warm (T$_\mathrm{eq, A_B=0.3} \sim 620$K) sub-Neptune where atmospheric absorbers have been detected to date. Our two transit observations with NIRISS/SOSS, combined with the existing HST/WFC3 spectrum, enable us to break the clouds-metallicity degeneracy. We detect water in a highly metal-enriched "steam world" atmosphere (O/H of $\sim 4$ by mass and H$_2$O found to be the background gas with a volume mixing ratio of >31%). We further show that these results are robust to stellar contamination through the transit light source effect. We do not detect escaping metastable He, which, combined with previous nondetections of escaping He and H, supports the steam atmosphere scenario. In water-rich atmospheres, hydrogen loss driven by water photolysis happens predominantly in the ionized form which eludes observational constraints. We also detect several flares in the NIRISS/SOSS light-curves with far-UV energies of the order of 10$^{30}$ erg, highlighting the active nature of the star. Further atmospheric characterization of GJ 9827 d probing carbon or sulfur species could reveal the origin of its high metal enrichment., Comment: 37 pages, 18 figures, Accepted for publication in ApJL

Ultrafast all-optical logic devices based on nonlinear light-matter interactions hold the promise to overcome the speed limitations of conventional electronic devices. Strong coupling of excitons and photons inside an optical resonator enhances such interactions and generates new polariton states which give access to unique nonlinear phenomena, such as Bose-Einstein condensation, used for all-optical ultrafast polariton transistors. However, the pulse energies required to pump such devices range from tens to hundreds of pJ, making them not competitive with electronic transistors. Here we introduce a new paradigm for all-optical switching based on the ultrafast transition from the strong to the weak coupling regime in microcavities embedding atomically thin transition metal dichalcogenides. Employing single and double stacks of hBN-encapsulated MoS$_2$ homobilayers with high optical nonlinearities and fast exciton relaxation times, we observe a collapse of the 55-meV polariton gap and its revival in less than one picosecond, lowering the threshold for optical switching below 4 pJ per pulse, while retaining ultrahigh switching frequencies. As an additional degree of freedom, the switching can be triggered pumping either the intra- or the interlayer excitons of the bilayers at different wavelengths, speeding up the polariton dynamics, owing to unique interspecies excitonic interactions. Our approach will enable the development of compact ultrafast all-optical logical circuits and neural networks, showcasing a new platform for polaritonic information processing based on manipulating the light-matter coupling.

In CrX$_2$ (X=Br, I), a Jahn-Teller effect distorts the octahedral configuration of the anions about Cr, resulting in a "ribbon chain" structure. We previously observed helimagnetism in CrI$_2$ propagating along the ribbon chains, with an 80-90$^{\circ}$ rotation per Cr ion. Via neutron scattering, we report that CrBr$_2$ and the solid solution CrIBr are also helimagnetic, with N\'{e}el temperatures of 17 and 12 K, respectively, and signs of intermediate magnetic transitions in both compounds. The helical angle between spins on consecutive Cr ions along the ribbon chains increases with increasing Br substitution, from 89.7$^{\circ}$ (CrI$_2$) to 116$^{\circ}$ (CrIBr) to 147$^{\circ}$ (CrBr$_2$), possibly as a result of decreasing longer-range intrachain spin interactions with the increased electron localization of the lighter anions.

In the study of van der Waals-layered magnetic materials, the properties of CrCl$_3$ continue to attract attention. This compound is reported to undergo antiferromagnetic (AFM) ordering below $\sim$14 K, with a ferromagneticlike region proposed to exist between 14 and 17 K. Ideally, the crystal structure is rhombohedral (R) below $\sim$235 K, separated from a higher-temperature monoclinic (M) phase by a layer-sliding structural phase transition. However, the structural transition is often inhibited even in bulk single crystals, allowing M-type layer stacking, reported to have a tenfold greater interlayer magnetic coupling than R-type stacking, to be present at low temperature. To clarify the effect of stacking defects on CrCl$_3$, we report magnetization measurements on samples of varying crystalline quality. At low applied magnetic field, some crystals predominantly show the $T_N=14$ K peak, but other crystals show hysteretic behavior and a magnetization enhancement at a slightly higher temperature ($14 < T \lesssim 17$ K.) Samples with anomalous behavior exhibit a transition around $\sim$2 T in isothermal magnetization-field data, providing evidence that M-type stacking defects are the source of these anomalies. Ground powder samples are especially likely to show strongly anomalous behavior. We suggest that the anomalous behavior arises from few-layer magnetic domains that form just above $T_N$ in an environment of mixed interlayer magnetic coupling strength. We argue that the influence of M-type stacking boundaries on sublattice magnetization is already observable in reported neutron scattering data, and may be responsible for a certain feature in reported specific heat data., Comment: Supplement is part of pdf file

The recent observation of unconventional superconductivity in thin films of LaNiO2 (critical temperature, Tc~10 K) and in bulk single crystals of La3Ni2O7 (327) under pressure (Tc~80 K), has brought to light a long sought-after class of superconducting nickelates. Through structural measurements in the 327-system, it was shown that the absence of superconductivity is related to bending of the O-Ni-O bonds. Similarly, the bond bending may be linked to the absence of superconductivity in bulk LaNiO2. Neutron diffraction was used on bulk non-superconducting La(1-x)Sr(x)NiO(2) samples to show that the layers are naturally buckled, creating a Ni-O-Ni bond angle of 177 at 2 K and ambient pressure. The buckling angle increases to 170 on warming to room temperature. Furthermore, a broad paramagnetic continuum is observed that decreases in intensity on cooling from room temperature signaling a possible transition to a coherent state. However, no antiferromagnetic (AFM) peaks are detected, although enhancement of ferromagnetic (FM) correlations cannot be excluded., Comment: 5 pages, 4 figures

CrI$_{2}$ is a quasi-one dimensional (1D) van der Waals (vdW) system that exhibits helimagnetism that propagates along the ribbons. This was determined from neutron time-of-flight diffraction measurements. Below $T_N=17$ K, a screw-like helimagnetic order develops with an incommensurate wavevector of $\mathbf{q} \approx (0.2492,0,0)$ at 8 K. Using density functional theory (DFT)$+U$ calculations, the $J_{1}$-$J_{2}$ model was leveraged to describe the helimagnetism, where $J_{1} (> 0)$ and $J_2 (< 0)$ correspond, respectively, to a ferromagnetic nearest neighbor (NN) and antiferromagnetic next-nearest neighbor (NNN) intrachain interaction. The DFT$+U$ calculations predict that bulk CrI$_2$ in the orthorhombic $Cmc2_1$ crystal structure satisfies the $|J_2| > |J_1|/4$ condition, which favors formation of helimagnetic order., Comment: main pdf file includes supplement

In 1T-TaS$_{2-x}$Se$_x$, the charge density wave (CDW) state features a star of David lattice that expands across layers as the system becomes commensurate on cooling. The layers can also order along the c-axis and different stacking orders have been proposed. Using neutron scattering on powder samples, we compared the stacking order previously observed in 1T-TaS$_2$ as the system is doped with Se. While at low temperature, a 13c layer sequence stacking was observed in TaS$_2$, this type of ordering was not evident with doping. Doping with Se results in a nearly commensurate state with the Mott state suppressed which may be linked to the absence of the layer stacking., Comment: 5 pages, 3 figures

Bandstructure engineering using alloying is widely utilised for achieving optimised performance in modern semiconductor devices. While alloying has been studied in monolayer transition metal dichalcogenides, its application in van der Waals heterostructures built from atomically thin layers is largely unexplored. Here, we fabricate heterobilayers made from monolayers of WSe$_2$ (or MoSe$_2$) and Mo$_x$W$_{1-x}$Se$_2$ alloy and observe nontrivial tuning of the resultant bandstructure as a function of concentration $x$. We monitor this evolution by measuring the energy of photoluminescence (PL) of the interlayer exciton (IX) composed of an electron and hole residing in different monolayers. In Mo$_x$W$_{1-x}$Se$_2$/WSe$_2$, we observe a strong IX energy shift of $\approx$100 meV for $x$ varied from 1 to 0.6. However, for $x<0.6$ this shift saturates and the IX PL energy asymptotically approaches that of the indirect bandgap in bilayer WSe$_2$. We theoretically interpret this observation as the strong variation of the conduction band K valley for $x>0.6$, with IX PL arising from the K-K transition, while for $x<0.6$, the bandstructure hybridization becomes prevalent leading to the dominating momentum-indirect K-Q transition. This bandstructure hybridization is accompanied with strong modification of IX PL dynamics and nonlinear exciton properties. Our work provides foundation for bandstructure engineering in van der Waals heterostructures highlighting the importance of hybridization effects and opening a way to devices with accurately tailored electronic properties., Comment: Supporting Information can be found downloading and extracting the gzipped tar source file listed under "Other formats"

The recent observations on WASP-39 b by JWST have revealed hints of high metallicity within the atmosphere compared to its host star (Feinstein et al. 2022; Ahrer et al. 2023; Alderson et al. 2023; Rustamkulov et al. 2023; Tsai et al. 2023). There are various theories on how these high metallic atmospheres emerge. In this study, we closely investigate the impact of extreme escape in the form of hydrodynamic escape to see its impact on atmospheric metallicity and spectral features such as CH$_4$, CO$_2$, and SO$_2$. We perform a grid simulation, with an adapted version of MESA that includes hydrodynamic escape (Kubyshkina et al. 2018; 2020), to fully evolve planets with similar masses and radii to the currently observed WASP-39 b estimates. By making use of (photo-)chemical kinetics and radiative transfer codes, we evaluate the transmission spectra at various time intervals throughout the simulation. Our results indicate that the massive size of WASP-39 b limits the metal enhancement to a maximum of ~1.23 the initial metallicity. When incorporating metal drag, this enhancement factor is repressed to an even greater degree, resulting in an enrichment of at most ~0.4%. As a consequence, when assuming an initial solar metallicity, metal-enriched spectral features like SO$_2$ are still missing after ~9 Gyr into the simulation. This paper, thus, demonstrates that hydrodynamic escape cannot be the primary process behind the high metallicity observed in the atmosphere of WASP-39 b, suggesting instead that a metal-enhanced atmosphere was established during its formation., Comment: 11 pages, 3 figures, accepted for publication in ApJL

Bibliometric studies based on the Web of Science (WOS) database have become an increasingly popular method for analysing the structure of scientific research. So do network approaches, which, based on empirical data, make it possible to characterize the emergence of topological structures over time and across multiple research areas. Our paper is a contribution to interweaving these two lines of research that have progressed in separate ways but whose common applications have been increasingly more frequent. Among other attributes, Author Keywords and Keywords Plus are used as units of analysis that enable us to identify changes in the topics of interest and related bibliography. By considering the co-occurrence of those keywords with the Author Keyword \texttt{Complexity}, we provide an overview of the evolution of studies on Complexity Sciences, and compare this evolution in seven scientific fields. The results show a considerable increase in the number of papers dealing with complexity, as well as a general tendency across different disciplines for this literature to move from a more foundational, general and conceptual to a more applied and specific set of co-occurring keywords. Moreover, we provide evidence of changing topologies of networks of co-occurring keywords, which are described through the computation of some topological coefficients. In so doing, we emphasize the distinguishing structures that characterize the networks of the seven research areas., Comment: 25 pages, 7 figures

In this era of exoplanet characterisation with JWST, the need for a fast implementation of classical forward models to understand the chemical and physical processes in exoplanet atmospheres is more important than ever. Notably, the time-dependent ordinary differential equations to be solved by chemical kinetics codes are very time-consuming to compute. In this study, we focus on the implementation of neural networks to replace mathematical frameworks in one-dimensional chemical kinetics codes. Using the gravity profile, temperature-pressure profiles, initial mixing ratios, and stellar flux of a sample of hot-Jupiters atmospheres as free parameters, the neural network is built to predict the mixing ratio outputs in steady state. The architecture of the network is composed of individual autoencoders for each input variable to reduce the input dimensionality, which is then used as the input training data for an LSTM-like neural network. Results show that the autoencoders for the mixing ratios, stellar spectra, and pressure profiles are exceedingly successful in encoding and decoding the data. Our results show that in 90% of the cases, the fully trained model is able to predict the evolved mixing ratios of the species in the hot-Jupiter atmosphere simulations. The fully trained model is ~1000 times faster than the simulations done with the forward, chemical kinetics model while making accurate predictions., Comment: 13 pages, 9 figures, accepted for publication at MNRAS

Ground-based high-resolution spectroscopy (HRS) has detected numerous chemical species and atmospheric dynamics in exoplanets, most notably ultra-hot Jupiters (UHJs). However, quantitative estimates on abundances have been challenging but are essential for accurate comparative characterisation and to determine formation scenarios. In this work we retrieve the atmospheres of six UHJs (WASP-76~b, MASCARA-4~b, MASCARA-2~b, WASP-121~b, HAT-P-70~b and WASP-189~b) with ESPRESSO and HARPS-N/HARPS observations, exploring trends in eleven neutral species and dynamics. While Fe abundances agree well with stellar values, Mg, Ni, Cr, Mn and V show more variation, highlighting the difficulty in using a single species as a proxy for metallicity. We find that Ca, Na, Ti and TiO are under-abundant, potentially due to ionisation and/or night-side rain-out. Our retrievals also show that relative abundances between species are more robust, consistent with previous works. We perform spatially- and phase-resolved retrievals for WASP-76~b and WASP-121~b given their high signal-to-noise observations, and find the chemical abundances in each of the terminator regions are broadly consistent. We additionally constrain dynamics for our sample through Doppler shifts and broadening of the planetary signals during the primary eclipse, with median blue shifts between $\sim$0.9-9.0~km/s due to day-night winds. Furthermore, we constrain spectroscopic masses for MASCARA-2~b and HAT-P-70~b consistent with their known upper limits, but we note that these may be biased due to degeneracies. This work highlights the importance of future HRS studies to further probe differences and trends between exoplanets., Comment: 26 pages, 11 figures, 5 tables, published in AJ

Van der Waals (vdW) Dirac magnon system CrI$_3$, a potential host of topological edge magnons, orders ferromagnetically (FM) (T$_C=61$ K) in the bulk, but antiferromagnetic (AFM) order has been observed in nanometer thick flakes, attributed to monoclinic (M) type stacking. We report neutron scattering measurements on a powder sample where the usual transition to the rhombohedral (R) phase was inhibited for a majority of the structure. Elastic measurements (and the opening of a hysteresis in magnetization data on a pressed pellet) showed that an AFM transition is clearly present below $\sim$50 K, coexisting with the R-type FM order. Inelastic measurements showed a decrease in magnon energy compared to the R phase, consistent with a smaller interlayer magnetic coupling in M-type stacking. A gap remains at the Dirac point, suggesting that the same nontrivial magnon topology reported for the R phase may be present in the M phase as well., Comment: includes supplement

A key first step to constrain the impact of energetic particles in exoplanet atmospheres is to detect the chemical signature of ionisation due to stellar energetic particles and Galactic cosmic rays. We focus on GJ$\,$436, a well-studied M dwarf with a warm Neptune-like exoplanet. We demonstrate how the maximum stellar energetic particle momentum can be estimated from the stellar X-ray luminosity. We model energetic particle transport through the atmosphere of a hypothetical exoplanet at orbital distances between $a=0.01-0.2\,$au from GJ$\,$436, including GJ$\,$436$\,$b's orbital distance (0.028$\,$au). For these distances we find that, at top-of-atmosphere, stellar energetic particles ionise molecular hydrogen at a rate of $\zeta_{\rm StEP,H_2} \sim 4\times10^{-10}-2\times10^{-13}\,\mathrm{s^{-1}}$. In comparison, Galactic cosmic rays alone lead to $\zeta_{\rm GCR, H_2}\sim2\times 10^{-20}-10^{-18} \,\mathrm{s^{-1}}$. At 10au we find that ionisation due to Galactic cosmic rays equals that of stellar energetic particles: $\zeta_{\rm GCR,H_2} = \zeta_{\rm StEP,H_2} \sim 7\times10^{-18}\,\rm{s^{-1}}$ for the top-of-atmosphere ionisation rate. At GJ$\,$436$\,$b's orbital distance, the maximum ion-pair production rate due to stellar energetic particles occurs at pressure $P\sim 10^{-3}\,$bar while Galactic cosmic rays dominate for $P>10^2\,$bar. These high pressures are similar to what is expected for a post-impact early Earth atmosphere. The results presented here will be used to quantify the chemical signatures of energetic particles in warm Neptune-like atmospheres., Comment: 14 pages, 9 figures, accepted for publication in MNRAS

Topological quantum magnets FeSn and Fe$_{3}$Sn$_{2}$ were studied using neutron scattering and first-principles calculations. Both materials are metallic but host dispersionless flat bands with Dirac nodes at the $K$ point in reciprocal space. The local structure determined from the pair density function analysis of the neutron diffraction data provided no evidence for electron localization in both compounds, consistent with their metallic nature. At the same time, in FeSn, an anomalous suppression in the $c$-axis lattice constant coupled with changes in the phonon spectra were observed across T$_{N}$ indicating the presence of magneto-elastic coupling and spin-phonon interactions. In addition, it was observed that spin waves persisted well above T$_{N}$, suggesting that the in-plane ferromagnetic spin correlations survive at high temperatures. In contrast, no lattice anomaly was observed in Fe$_{3}$Sn$_{2}$. The inelastic signal could be mostly accounted for by phonons, determined from density functional theory, showing typical softening on warming., Comment: 10 pages, 7 figures

This manuscript outlines a concept that would leverage the existing proton accelerator at the Spallation Neutron Source of Oak Ridge National Laboratory to enable transformative science via one world-class facility serving two missions: Single Event Effects (SEE) and Muon Spectroscopy ($\mu$SR). The $\mu$SR portion would deliver the world's highest flux and highest resolution pulsed muon beams for materials characterization purposes, with precision and capabilities well beyond comparable facilities. The SEE capabilities deliver neutron, proton and muon beams for aerospace industries that are facing an impending challenge to certify equipment for safe and reliable behavior under bombardment from atmospheric radiation originating from cosmic and solar rays. With negligible impact on the primary neutron scattering mission of the SNS, the proposed facility will have enormous benefit for science and industry alike. We have designated this facility 'SEEMS'., Comment: 11 pages, 7 figures, to be published in Review of Scientific Instruments

Photochemistry is a fundamental process of planetary atmospheres that regulates the atmospheric composition and stability. However, no unambiguous photochemical products have been detected in exoplanet atmospheres to date. Recent observations from the JWST Transiting Exoplanet Early Release Science Program found a spectral absorption feature at 4.05 $\mu$m arising from SO$_2$ in the atmosphere of WASP-39b. WASP-39b is a 1.27-Jupiter-radii, Saturn-mass (0.28 M$_J$) gas giant exoplanet orbiting a Sun-like star with an equilibrium temperature of $\sim$1100 K. The most plausible way of generating SO$_2$ in such an atmosphere is through photochemical processes. Here we show that the SO$_2$ distribution computed by a suite of photochemical models robustly explains the 4.05 $\mu$m spectral feature identified by JWST transmission observations with NIRSpec PRISM (2.7$\sigma$) and G395H (4.5$\sigma$). SO$_2$ is produced by successive oxidation of sulphur radicals freed when hydrogen sulphide (H$_2$S) is destroyed. The sensitivity of the SO$_2$ feature to the enrichment of the atmosphere by heavy elements (metallicity) suggests that it can be used as a tracer of atmospheric properties, with WASP-39b exhibiting an inferred metallicity of $\sim$10$\times$ solar. We further point out that SO$_2$ also shows observable features at ultraviolet and thermal infrared wavelengths not available from the existing observations., Comment: 39 pages, 14 figures, accepted to be published in Nature

Transmission spectroscopy of exoplanets has revealed signatures of water vapor, aerosols, and alkali metals in a few dozen exoplanet atmospheres. However, these previous inferences with the Hubble and Spitzer Space Telescopes were hindered by the observations' relatively narrow wavelength range and spectral resolving power, which precluded the unambiguous identification of other chemical species$-$in particular the primary carbon-bearing molecules. Here we report a broad-wavelength 0.5-5.5 $\mu$m atmospheric transmission spectrum of WASP-39 b, a 1200 K, roughly Saturn-mass, Jupiter-radius exoplanet, measured with JWST NIRSpec's PRISM mode as part of the JWST Transiting Exoplanet Community Early Release Science Team program. We robustly detect multiple chemical species at high significance, including Na (19$\sigma$), H$_2$O (33$\sigma$), CO$_2$ (28$\sigma$), and CO (7$\sigma$). The non-detection of CH$_4$, combined with a strong CO$_2$ feature, favours atmospheric models with a super-solar atmospheric metallicity. An unanticipated absorption feature at 4$\mu$m is best explained by SO$_2$ (2.7$\sigma$), which could be a tracer of atmospheric photochemistry. These observations demonstrate JWST's sensitivity to a rich diversity of exoplanet compositions and chemical processes., Comment: 41 pages, 4 main figures, 10 extended data figures, 4 tables. Under review in Nature

In quasi-two-dimensional 1T-TaS2, a charge density wave (CDW) prototype, the transition occurs in two steps, from incommensurate (ICDW) and nearly commensurate (NCCDW) and from NCDW to commensurate (CCDW), locked in the resistivity step-wise behavior. The hysteresis observed in the resistivity across the NCDW-to-CCDW transition has roots to a structural hysteresis, where local distortions of the sqrt(13)a x sqrt(13)a superstructure are revealed by neutron and X-ray diffraction. The structural hysteresis is due to faulty stars of David (SODs) because of Ta displacements away from the perfect trigonal geometry as well as out of plane S distortions. Furthermore, the superstructure exhibits a 3co layer stacking order that weakens on warming and fully disappears in the ICDW state., Comment: 5 pages including 1 page of References, 3 figures

The layer stacking order in 2D materials strongly affects functional properties and holds promise for next generation electronic devices. In bulk, octahedral MoTe$_2$ possesses two stacking arrangements, the Weyl semimetal T$_d$ phase, and the higher-order topological insulator 1T' phase; however, it remains unclear if thin exfoliated flakes of MoTe$_2$ follow the T$_d$, 1T', or an alternative stacking sequence. Here, we resolve this debate using atomic-resolution imaging within the transmission electron microscope. We find that the layer stacking in thin flakes of MoTe$_2$ is highly disordered and pseudo-random, which we attribute to intrinsic confinement effects. Conversely, WTe$_2$, which is isostructural and isoelectronic to MoTe$_2$, displays ordered stacking even for thin exfoliated flakes. Our results are important for understanding the quantum properties of MoTe$_2$ devices, and suggest that thickness may be used to alter the layer stacking in other 2D materials., Comment: 6 figures, 2 tables

Numerous optical phenomena and applications have been enabled by nanophotonic structures. Their current fabrication from high refractive index dielectrics, such as silicon or gallium phosphide, pose restricting fabrication challenges, while metals, relying on plasmons and thus exhibiting high ohmic losses, limit the achievable applications. Here, we present an emerging class of layered so-called van der Waals (vdW) crystals as a viable nanophotonics platform. We extract the dielectric response of 11 mechanically exfoliated thin-film (20-200 nm) van der Waals crystals, revealing high refractive indices up to n = 5, pronounced birefringence up to $\Delta$n = 3, sharp absorption resonances, and a range of transparency windows from ultraviolet to near-infrared. We then fabricate nanoantennas on SiO$_2$ and gold utilizing the compatibility of vdW thin films with a variety of substrates. We observe pronounced Mie resonances due to the high refractive index contrast on SiO$_2$ leading to a strong exciton-photon coupling regime as well as largely unexplored high-quality-factor, hybrid Mie-plasmon modes on gold. We demonstrate further vdW-material-specific degrees of freedom in fabrication by realizing nanoantennas from stacked twisted crystalline thin-films, enabling control of nonlinear optical properties, and post-fabrication nanostructure transfer, important for nano-optics with sensitive materials.

M-dwarfs are thought to be hostile environments for exoplanets. Stellar events are very common on such stars. These events might cause the atmospheres of exoplanets to change significantly over time. It is not only the major stellar flare events that contribute to this disequilibrium, but the smaller flares might also affect the atmospheres in an accumulating manner. In this study, we aim to investigate the effects of time-dependent stellar activity on the atmospheres of known exoplanets. We simulate the chemistry of GJ876c, GJ581c, and GJ832c that go from H$_2$-dominated to N$_2$-dominated atmospheres using observed stellar spectra from the MUSCLES-collaboration. We make use of the chemical kinetics code VULCAN and implement a flaring routine that stochastically generates synthetic flares based on observed flare statistics. Using the radiative transfer code petitRADTrans we also simulate the evolution of emission and transmission spectra. We investigate the effect of recurring flares for a total of 11 days covering 515 flares. Results show a significant change in abundance for some relevant species such as H, OH, and CH$_4$, with factors going up to 3 orders of magnitude difference with respect to the preflare abundances. We find a maximum change of $\sim$12 ppm for CH$_4$ in transmission spectra on GJ876c. These changes in the spectra remain too small to observe. We also find that the change in abundance and spectra of the planets accumulate throughout time, causing permanent changes in the chemistry. We conclude this small but gradual change in chemistry arises due to the recurring flares., Comment: Accepted for publication in MNRAS

Nonlinear interactions between excitons strongly coupled to light are key for accessing quantum many-body phenomena in polariton systems. Atomically-thin two-dimensional semiconductors provide an attractive platform for strong light-matter coupling owing to many controllable excitonic degrees of freedom. Among these, the recently emerged exciton hybridization opens access to unexplored excitonic species, with a promise of enhanced interactions. Here, we employ hybridized interlayer excitons (hIX) in bilayer MoS2 to achieve highly nonlinear excitonic and polaritonic effects. Such interlayer excitons possess an out-of-plane electric dipole as well as an unusually large oscillator strength allowing observation of dipolar polaritons(dipolaritons) in bilayers in optical microcavities. Compared to excitons and polaritons in MoS2 monolayers, both hIX and dipolaritons exhibit about 8 times higher nonlinearity, which is further strongly enhanced when hIX and intralayer excitons, sharing the same valence band, are excited simultaneously. This gives rise to a highly nonlinear regime which we describe theoretically by introducing a concept of hole crowding. The presented insight into many-body interactions provides new tools for accessing few-polariton quantum correlations.