Search NASA⌕ Search

SEARCH · Search NASA

Results for “first principles simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Double-Atom Catalysts Featuring Inverse Sandwich Structure for CO 2 Reduction Reaction: A Synergetic First-Principles and Machine Learning Investigation

Electrocatalytic CO 2 reduction reactions (CO 2 RR) based on scalable and highly efficient catalysis provide an attractive strategy for reducing CO 2 emissions. Here in this work, we combined first-principles density functional theory (DFT) and machine learning (ML) to comprehensively explore the potential of double-atom catalysts (DACs) featuring an inverse sandwich structure anchored on defective graphene (gra) to catalyze CO 2 RR to generate C 1 products. We started with five homonuclear M 2 ⊥gra (M = Co, Ni, Rh, Ir, and Pt), followed by 127 heteronuclear MM'⊥gra (M = Co, Ni, Rh, Ir, and Pt, M' = Sc–Au). Stable DACs were screened by evaluating their binding energy, formation energy, and dissolution potential of metal atoms, as well as conducting first-principles molecular dynamics simulations with and without solvent water molecules. Based on DFT calculations, Rh 2 ⊥gra DAC was found to outperform the other four homonuclear DACs and the Rh-based single- and double-atom catalysts of noninverse sandwich structures. Out of the 127 heteronuclear DACs, 14 were found to be stable and have good catalytic performance. An ML approach was adopted to correlate key factors with the activity and stability of the DACs, including the sum of radii of metal and ligand atoms (d M–M' , d M–C , and d M'–C ), the sum and difference of electronegativity of two metal atoms (P M + P M' , P M – P M '), the sum and difference of first ionization energy of two metal atoms (I M + I M' , I M – I M '), the sum and difference of electron affinity of two metal atoms (A M + A M' , A M – A M '), and the number of d-electrons of the two metal atoms (Nd). The obtained ML models were further used to predict 154 potential electrocatalysts out of 784 possible DACs featuring the same inverse sandwich configuration. Overall, this work not only identified promising CO 2 RR DACs featuring the reported inverse sandwich structure but also provided insights into key atomic characteristics associated with high CO 2 RR activity.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Density of Fe-Ni-C Liquids at High Pressures and Implications for Liquid Cores of Earth and the Moon

The presence of light elements in the metallic cores of the Earth, the Moon, and other rocky planetary bodies has been widely proposed. Carbon is among the top candidates in light of its high cosmic abundance, siderophile nature, and ubiquity in iron meteorites. It is, however, still controversial whether carbon-rich core compositional models can account for the seismic velocity observations within the Earth and lunar cores. In this paper we report the density and elasticity of Fe 90 Ni 10 -3 wt.% C and Fe 90 Ni 10 -5 wt.% C liquid alloys using synchrotron-based X-ray absorption experiments and first-principles molecular dynamics simulations. Our results show that alloying of 3 wt.% and 5 wt.% C lowers the density of Fe 90 Ni 10 liquid by ~2.9-3.1% at 2 GPa, and ~3.4-3.6% at 9 GPa. More intriguingly, our experiments and simulations both demonstrate that the bulk moduli of the Fe-Ni-C liquids are similar to or slightly higher than those of Fe-Ni liquids. Thus, the calculated compressional velocity ($v_p$) of Fe-Ni-C liquids are higher than that of pure Fe-Ni alloy, promoting carbon as a possible candidate to explain the elevated v p in the Earth's outer core. However, the values and slopes of both density and v p of the studied two Fe-Ni-C liquids do not match the outer core seismic models, suggesting that carbon may not be the sole principal light element in Earth's outer core. The high $v_p$ of Fe-Ni-C liquids does not match the presumptive $v_p$ of the lunar outer core well, indicating that carbon is less likely to be its dominant light element.

58 GEOSCIENCES↗
NASA NTRSDate not supplied

Mixed-Domain Charge Transport in S-Se Alloys as a Li-S Battery Cathode Material

Lithium-sulfur batteries are emerging candidate systems for the next-generation long-range electric aircraft application owing to their potential to deliver high energy density per weight. Their cathodes are to be made primarily with sulfur, but sulfur by itself has much too low electron conductivity to serve as a useful cathode. An idea that has been suggested to overcome this bottleneck is to alloy sulfur with selenium in the hope that electron mobility, which is thought to improve charge transport at some expense of increased weight and reduced energy density. However, understanding of these alloy structures and their transport mechanism is insufficient, and electron mobility at varying degrees of selenium content is not well charted, which are critical for determining the optimum selenium content and testing whether this strategy is generally feasible. The difficulty of computationally characterizing these alloy systems is rooted in the structure. The structures of both sulfur and selenium exhibit eight-atom rings that pack together in various orientations to form their respective crystals. For one, because these crystals (and presumably their alloys as well) feature multiple polymorphs such that their structural coherence is rather unclear. Secondly, these are also relatively large systems (32 atoms per cell) with low-symmetry, which complicate computation both in terms of accuracy and efficiency. Thirdly, such semi-molecular, semi-crystalline structures lead to a combination of band-transport characteristics and hopping-transport characteristics, each of which is a domain with its own physics and set of computational, theoretical challenges. In addition, there is a general dearth of experimental transport data for sulfur-selenium alloys. In this study, we make a comprehensive attempt to tackle this problem using a wide array of first-principles methods. We generate special quasirandom structures to simulate alloy structures, and compute their first-principles Raman spectra for comparison with experimental data to ensure their structural soundness. We then proceed to use recently developed, state-of-the-art tool (AMSET) in order to efficiently compute electronic mobilities under band transport of pure sulfur, pure selenium, as well as their alloys of various compositions at a reasonable accuracy. We then sample numerous dimer configurations of nearest-neighbor eight-atom-ring-pairs and compute electronic hopping rates between them, which yields hopping mobility. A combination of these efforts lead to a general mapping of electron transport behaviors throughout the alloy range. Preliminary results show that introduction of selenium into sulfur initially damages electron mobility due to disorder but eventually improves it beyond that of pure sulfur with additional selenium content, ultimately peaking at the pure-selenium limit. Band transport dominates for holes in sulfur and electrons in selenium, but in all other cases, hopping transport is dominant. Ongoing efforts include determination of charge concentration and conductivity, as well as performing multiphysics modeling to determine the optimum selenium content for best tradeoff between conductivity and energy density for aircraft range.

Junsoo Park↗
DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accurate determination of solvation free energies of neutral organic compounds from first principles

The main goal of molecular simulation is to accurately predict experimental observables of molecular systems. Another long-standing goal is to devise models for arbitrary neutral organic molecules with little or no reliance on experimental data. While separately these goals have been met to various degrees, for an arbitrary system of molecules they have not been achieved simultaneously. For biophysical ensembles that exist at room temperature and pressure, and where the entropic contributions are on par with interaction strengths, it is the free energies that are both most important and most difficult to predict. We compute the free energies of solvation for a diverse set of neutral organic compounds using a polarizable force field fitted entirely to ab initio calculations. The mean absolute errors (MAE) of hydration, cyclohexane solvation, and corresponding partition coefficients are 0.2 kcal/mol, 0.3 kcal/mol and 0.22 log units, i.e. within chemical accuracy. The model (ARROW FF) is multipolar, polarizable, and its accompanying simulation stack includes nuclear quantum effects (NQE). The simulation tools’ computational efficiency is on a par with current state-of-the-art packages. The construction of a wide-coverage molecular modelling toolset from first principles, together with its excellent predictive ability in the liquid phase is a major advance in biomolecular simulation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Sequential Infiltration Synthesis of Cadmium Sulfide Discrete Atom Clusters

Abstract Exposure of soft material templates to alternating volatile chemical precursors can produce inorganic deposition within the permeable template (e.g. a polymer thin film) in a process akin to atomic layer deposition (ALD). While such sequential infiltration synthesis (SIS) processes have now been demonstrated for many metal oxides, we report an SIS process for a transition metal sulfide – CdS. Gas phase dimethyl cadmium and hydrogen sulfide precursors infiltrated into poly(4‐vinylpyridine) thin films result in the 3D‐nucleation of clusters consistent with a cubane‐type Cd 4 S 4 core that are variably terminated with methyl, thiol and hydroxy capping ligands. First principles models and simulation of few‐atom Cd‐based clusters are consistent with electronic and vibrational spectroscopy and grazing‐incidence total X‐ray scattering measurements of 3D‐cluster‐arrays synthesized at 80 °C. The direct synthesis of few‐atom transition metal sulfide clusters within polymer thin films will provide a versatile new route to precision architectures for light‐absorbing materials including solar energy harvesting and conversion applications.

Jayaweera, Nuwanthaka P. [Material Science Divisio↗
DOE OSTI2025

Sequential Infiltration Synthesis of Cadmium Sulfide Discrete Atom Clusters

Exposure of soft material templates to alternating volatile chemical precursors can produce inorganic deposition within the permeable template (e.g. a polymer thin film) in a process akin to atomic layer deposition (ALD). While such sequential infiltration synthesis (SIS) processes have now been demonstrated for many metal oxides, we report an SIS process for a transition metal sulfide – CdS. Gas phase dimethyl cadmium and hydrogen sulfide precursors infiltrated into poly(4-vinylpyridine) thin films result in the 3D-nucleation of clusters consistent with a cubane-type Cd 4 S 4 core that are variably terminated with methyl, thiol and hydroxy capping ligands. First principles models and simulation of few-atom Cd-based clusters are consistent with electronic and vibrational spectroscopy and grazing-incidence total X-ray scattering measurements of 3D-cluster-arrays synthesized at 80 °C. The direct synthesis of few-atom transition metal sulfide clusters within polymer thin films will provide a versatile new route to precision architectures for light-absorbing materials including solar energy harvesting and conversion applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Tailoring the Angular Mismatch in MoS 2 Homobilayers through Deformation Fields

Ultrathin MoS 2 has shown remarkable characteristics at the atomic scale with an immutable disorder to weak external stimuli. Ion beam modification unlocks the potential to selectively tune the size, concentration, and morphology of defects produced at the site of impact in 2D materials. Combining experiments, first-principles calculations, atomistic simulations, and transfer learning, it is shown that irradiation-induced defects can induce a rotation-dependent moiré pattern in vertically stacked homobilayers of MoS 2 by deforming the atomically thin material and exciting surface acoustic waves (SAWs). Furthermore, the direct correlation between stress and lattice disorder by probing the intrinsic defects and atomic environments are demonstrated. The method introduced in this paper sheds light on how engineering defects in the lattice can be used to tailor the angular mismatch in van der Waals (vdW) solids.

2D materials↗
DOE OSTI2025

Speciation and diffusive dynamics in hydrated grain boundaries of complex oxide Gd2Ti2O7

Abstract Grain boundaries in polycrystalline materials significantly affect their properties, such as ionic transport, corrosion, and chemical durability. The pyrochlore compound (Gd 2 Ti 2 O 7 ) is employed as a model for complex oxides and is known for its diverse applications, including nuclear waste immobilization. Density functional theory-based first-principles molecular dynamics simulations were performed at different temperatures on the hydrated grain boundary system. The results show extensive transformations within the grain boundaries among hydrous water species (OH − , H 2 O, and H 3 O + ). The temperature dependence of self-diffusion coefficients follows Arrhenius behavior, with an activation energy of 35.9 kJ/mol for hydrogen and 46.3 kJ/mol for oxygen. The lifetime of OH − is about three to four times longer than that of H 2 O at temperatures from 800 to 2100 K, suggesting the greater stability of OH − over H 2 O, a unique characteristic of the grain boundaries. The estimated lifetime of the hydrous species decreases as the temperature increases, with an activation energy of 9.9 kJ/mol for OH − and 13.4 kJ/mol for H 2 O. While Gd 3 + is more mobile than Ti 4+ , both the Gd 3 + and Ti 4+ cations are orders of magnitude less mobile than the water species. The results suggest that water species are much more mobile within grain boundaries than in the bulk crystal and have the potential to penetrate deep into polycrystalline materials through grain boundaries, leading to grain boundary degradation and dissolution. The different mobilities of cations in complex oxides can lead to leaching of certain cations and incongruent dissolution during the chemical weathering of Earth and industrial materials.

B. Ghosh, Dipta↗
DOE OSTI2023

Dynamics of collision-induced energy transfer

Collisions are of fundamental importance in many relevant physico-chemical phenomena in the gas phase and at the gas-surface interface. Collision-induced energy transfer not only offers a sensitive probe of interatomic interactions, but also impacts an array of processes ranging from transport to reactivity. This perspective highlights some recent advances on collision-induced energy transfer dynamics. In the gas phase, cold collisions revealed remarkable quantum effects and stereodynamics of bimolecular inelastic scattering, which motivated development of new state-to-state quantum scattering methods. For surface processes, experiments of hydrogen atom scattering from various surfaces combined with first-principles molecular dynamics simulations shed valuable light on the relative importance of adiabatic energy transfer to surface phonons and nonadiabatic energy transfer to surface electrons. We outline further challenges in these fields and expect more fruitful interplay between experiment and theory.

Cold collision↗
DOE OSTI2023

Ab initio study of the structure and properties of amorphous silicon hydride from accelerated molecular dynamics simulations

This paper presents a large-scale ab initio simulation study of amorphous silicon hydride (a-Si 1-x H x ) with an emphasis on the structure and properties of the material across a range of hydrogen concentration by combining accelerated molecular dynamics (MD) simulations with first-principles density-functional calculations. The accelerated MD scheme relied on classical metadynamics, which enabled the development of 2500+ high-quality structural models of a-Si 1-x H x , with system sizes ranging from 150 to 6000 atoms and hydrogen concentrations vary from 6 to 20 at. %. The resulting amorphous networks were found to be completely free from any coordination defects and that they all exhibited a pristine band-gap in their electronic spectrum. The microstructural properties of hydrogen distributions were examined with an emphasis on the presence of isolated and clustered environments of hydrogen atoms. The results were compared with experimental data obtained from X-ray diffraction, infrared spectroscopy and nuclear magnetic resonance studies.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Analytic model for neutral penetration and plasma fueling

Neutral atoms recycled from wall interaction interact with confined plasma, thereby refuelling it, most strongly in the region closest to the wall. This occurs near the X-point in diverted configurations, or else near the wall itself in limited configurations. A progression of analytic models is developed for neutral density in the vicinity of a planar or linear source in an ionising domain. First-principles neutral transport simulations with DEGAS2 are used throughout to test the validity and limits of the model when using equivalent sources. The model is further generalised for strong plasma gradients or the inclusion of charge exchange. An important part of the problem of neutral fuelling from recycling is thereby isolated and solved with a closed-form analytic model. A key finding is that charge exchange with the confined plasma can be significantly simplified with a reasonable sacrifice of accuracy by treating it as a loss. The several assumptions inherent to the model (and the simulations with which it is compared) can be adapted according to the particular behaviour of neutrals in the divertor and the manner in which they cross the separatrix.

fusion plasma↗
DOE OSTI2022

How Water Attacks MXene

Two-dimensional (2D) transition metal carbides and nitrides (MXenes) have shown outstanding performances in electrochemical energy storage and many other applications. However, the stability of MXene remains a concern, especially its quick degradation in aqueous solutions under ambient conditions. Here, we report on the water/Ti 3 C 2 O 2 -MXene interfacial chemistry from first-principles molecular dynamics simulations at room temperature. Surprisingly, we find that the water molecules can attack the basal plane of Ti 3 C 2 O 2 and pull the surface Ti atoms out, thereby reconstructing the surface. Further, by tracking close encounters of water molecules and surface Ti atoms on the basal plane of Ti 3 C 2 O 2 , we show that the attack is initiated by the chemisorption of a water molecule on a surface Ti atom, followed by the breaking of Ti–C bonds and deprotonation of the water molecule, leading to the formation of Ti–OH on the Ti 3 C 2 O 2 surface and a hydronium ion in the aqueous phase. Our finding highlights the susceptibility of Ti 3 C 2 O 2 MXene to water attack, supporting recent experimental observations. Furthermore, we demonstrate that preventing close encounters of water molecules and the surface Ti atoms is key to the stability of the basal plane and can be realized by negatively charging the surface (thereby reorienting the O atoms of water away from the surface) or converting the surface O to -OH groups (thereby shifting the water layer further away from the surface). Our insights and approach highlight the importance of the reactivity of water when interfacing with 2D materials such as MXenes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Charge Transport in Solvated Donor–Acceptor Functionalized Peptoids: Molecular Dynamics and Rate Theory

Scalable solar-energy conversion requires photoactive materials that combine the efficiency of natural photosynthetic systems with the stability and processability needed for practical applications. Achieving reliable charge transport in soft, self-assembled organic materials remains challenging, as structural fluctuations and environmental effects strongly influence charge-transfer (CT) rates. Here, we present a broadly applicable computational framework for evaluating CT rates in the condensed phase, combining Fermi’s golden rule rate theory with inputs from all-atom molecular dynamics (MD) simulations and first-principles electronic-structure calculations. The approach does not rely on system-specific parametrization and is applicable to a wide range of soft and disordered materials. We demonstrate the applicability and usefulness of the framework on redox-active peptoids functionalized with iron–porphyrin (Fe–P) complexes, a bioinspired platform with programmable donor–acceptor units and tunable three-dimensional organization. The calculated CT rates exhibit strong sensitivity to molecular conformation, with variations spanning several orders of magnitude. This dependence is shown to arise from the pronounced variation in diabatic electronic coupling with the relative orientations and separations of the Fe–P complexes across the conformational ensemble. The framework provides a consistent route for connecting atomistic structure to CT kinetics in the condensed phase and enables analysis of structure–rate relations in organic semiconducting systems.

Charge transfer↗
DOE OSTI2025

Computational Study of Indium Oxide Photoelectrodes

Using a combination of first-principles molecular dynamics simulations and electronic structure calculations, we characterize the atomistic structure and vibrational properties of a photocatalytic surface of In 2 O 3 , a promising photoelectrode for the production of hydrogen peroxide. We then investigate the surface in contact with water and show that the electronic states of In 2 O 3 are appropriately positioned in energy to facilitate the two-electron water oxidation reaction (WOR) over the competing four-electron oxygen evolution reaction. In conclusion, we further propose that the use of strained thin films interfaced with water is beneficial in decreasing the optical gap of In 2 O 3 and thus in utilizing a wider portion of the solar spectrum for the WOR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

The Local Electronic Structure of Supercritical CO 2 from X-ray Raman Spectroscopy and Atomistic-Scale Modeling

Supercritical CO 2 is encountered in several technical and natural systems related to biology, geophysics, and engineering. While the structure of gaseous CO 2 has been studied extensively, the properties of supercritical CO 2 , particularly close to the critical point, are not well-known. Here, in this work, we combine X-ray Raman spectroscopy, molecular dynamics simulations, and first-principles density functional theory (DFT) calculations to characterize the local electronic structure of supercritical CO 2 at conditions around the critical point. The X-ray Raman oxygen K-edge spectra manifest systematic trends associated with the phase change of CO 2 and the intermolecular distance. Extensive first-principles DFT calculations rationalize these observations on the basis of the 4sσ Rydberg state hybridization. X-ray Raman spectroscopy is found to be a sensitive tool for characterizing electronic properties of CO 2 under challenging experimental conditions and is demonstrated to be a unique probe for studying the electronic structure of supercritical fluids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

ZIF-8 Vibrational Spectra: Peak Assignments and Defect Signals

In this study, Zeolitic Imidazolate Framework-8 (ZIF-8) is a promising material for gas separation applications. It also serves as a prototype for numerous ZIFs, including amorphous ones, with a broader range of possible applications, including sensors, catalysis, and lithography. It consists of zinc coordinated with 2-methylimidazolate (2mIm) and has been synthesized with methods ranging from liquid-phase to solvent-free synthesis, which aim to control its crystal size and shape, film thickness and microstructure, and incorporation into nanocomposites. Depending on the synthesis method and post-synthesis treatments, ZIF-8 materials may deviate from the nominal defect-free ZIF-8 crystal structure due to defects like missing 2mIm, missing zinc, and physically adsorbed 2mIm trapped in the ZIF-8 pores, which may alter its performance and stability. Infra-red (IR) spectroscopy has been used to assess the presence of defects in ZIF-8 and related materials. However, conflicting interpretations by various authors persist in the literature. Here, we systematically investigate ZIF-8 vibrational spectra by combining experimental IR spectroscopy and first-principles molecular dynamics simulations, focusing on assigning peaks and elucidating the spectroscopic signals of putative defects present in ZIF-8 material. We attempt to resolve conflicting assignments from the literature and to provide a comprehensive understanding of the vibrational spectra of ZIF-8 and its defect-induced variations, aiming towards more precise quality control and design of ZIF-8-based materials for emerging applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Synthetic Tuning of Exciton–Phonon Coupling in Janus WS 2(1-x) Se 2x Monolayers Revealed by Resonant Raman Excitation Spectroscopy for Optoelectronic Applications

Janus monolayers, such as WSSe, have broken out-of-plane symmetry and an intrinsic dipole moment, impacting exciton transport, lifetime, and phonon interactions while imbuing piezoelectric, photocatalytic, and Rashba spin-splitting properties to transition metal dichalcogenides (TMDs). The new properties of this atomically thin material can be used for optoelectronic device applications. As TMDs are converted into Janus monolayers, e.g., top selenization of WS2 to WSSe, the bandgap and structure smoothly evolve, impacting not only the formation of excitons but also their complex interactions with different phonon modes. Resonant Raman excitation profiles (REPs) are uniquely well-suited to reveal both excitonic transitions and exciton–phonon coupling. Here, the resonant REPs of $A^{'}_{1}$ WS 2 and A 1 WSSe modes are measured to understand the strength of their coupling with the A, B, and C excitonic bands of a WS2 monolayer throughout its stepwise transformation into Janus WSSe by pulsed laser deposition (PLD) of energetic selenium species. In situ Raman spectroscopy during deposition is used to controllably prepare stable intermediate Janus structures, WS 2(1-x) Se 2x (0 ≤ x ≤ 0.5), for ex situ measurement of their resonant REPs. As x increases, REPs reveal not only pronounced excitonic bands that gradually shift toward lower photon energies but also strong, mode-selective exciton–phonon coupling. First-principles resonant Raman simulations independently predict this spectral behavior and are shown capable of matching the spectrally broadened, experimentally observed REP profiles in this model system, indicating their strong predictive capability for future experiments. The combination of controlled synthesis, REP characterization, and predictive theory employed here demonstrates a powerful pathway to understand and ultimately tune exciton–phonon interactions for future quantum optical devices.

Janus monolayers↗