2026-04-17 | | Total: 31
By applying our methodology, we propose a defect in 4H-SiC which combines a Si vacancy and a C atom substituted with S (VSiSC) to have a spin-triplet ground state with the spin qubit functionality. Our calculations confirm that all configurations of the defect have a dynamically and thermodynamically stable triplet ground state and higher energy singlet states, essential for the spin-qubit polarization cycle. From GW calculations, we found that the electronic states associated with the defect form sharp and isolated peaks within the band gap for both triplet and singlet states. Further Bethe-Salpeter-equation calculations show that all considered configurations have intense optical excitations in the near infrared spectrum range. Analysis of the excitation energies and rates indicate that the VSiSC defect can be an excellent optically controlled spin qubit. Crucially, the host elements and the dopant have high-abundance isotopes with zero nuclear spin ensuring high spin-coherence time of the qubit.
In this study, we present a comprehensive analysis of the Raman active phonon modes in orthorhombic LaInO$_3$ based on a combination of polarization-angle resolved Raman spectroscopy and density functional theory calculations. By using backscattering from multiple crystallographic surface orientations and employing a full symmetry analysis, we identify and assign most of the Raman-active $Γ$-point phonons to their irreducible representations of the D$_{\rm{2h}}$ point group. A multidimensional hyperspectral fitting procedure allows us to extract the relative Raman tensor elements from the angular dependence of the scattering intensities, even for strongly overlapping modes. First-principles calculations yield the phonon dispersion along high-symmetry directions, the phonon densities of states, and atomic displacement patterns, which are found to be in good agreement with the experimental mode frequencies.
A limited remanent polarization (Pr) in HfO2-based FeRAM remains a key obstacle to density scaling and reliability, while material and process optimizations offer only incremental improvements. This limitation fundamentally originates from the thickness-constrained switchable polarization and the intrinsic polarization ceiling of HfO2-based ferroelectrics. Here, we propose an all-ALD-grown vertical complementary FeRAM (VCF) architecture, in which the top and bottom stacked FeRAM cells maintain complementary polarization. This complementary dipole configuration converts the readout from a single-layer polarization response into a differential polarization summation, thereby amplifying the effective charge window without increasing the switching field of each individual layer or incurring area overhead. Viewed from top to bottom, an "up-down" polarization pair stores logic '1', whereas a "down-up" pair stores logic '0'. Using a complementary polarization write-read scheme, the VCF achieves an effective differential polarization above 100 uC/cm^2 and retains above 90 uC/cm^2 after 1e10 switching cycles without electrical breakdown. Robust retention (longer than 1e4 s at 85 degC) and strong disturb immunity are demonstrated, with an effective differential polarization above 80 uC/cm^2 under a V/3 scheme after 1e6 disturb pulses. Array-level operation is validated in a 5 x 5 selector-free crosspoint array. The performance enhancement of the VCF arises from the co-optimization of the all-ALD-grown process, device architecture, and operation scheme, enabling high density, a wide memory window, and strong reliability for scalable FeRAM integration.
Since its discovery, ferroelectric hafnia has been extensively studied due to its CMOS-compatibility and ability to remain polarized at sub-10 nm thicknesses. The ferroelectric behaviour is generally attributed to a polar orthorhombic (OIII) phase. However, a second polar phase with rhombohedral symmetry (R-phase) has also been reported in epitaxial films. The nature of the R-phase remains disputed due to the subtle differences with the OIII-phase when probed by standard thin film characterisation techniques. Given the functional properties of ferroelectrics are crucially determined by the crystal symmetry, resolving this matter is imperative. In this work, we settle the controversy through extensive 3D reciprocal space surveys made possible via synchrotron-based grazing incidence diffraction from epitaxial films of both phases. These experiments, together with direct comparison of their temperature dependence and electrical responses, conclusively establish them as two distinct phases and provide insight into their key characteristics.
Magnetic skyrmions are emerging as promising candidates for next-generation information technologies, while the realization of scalable skyrmion lattices with tailored configurations is essential for advancing fundamental skyrmion physics and developing future applications. Here we achieved the controllable generation and regulation of a large-area, highly oriented skyrmion track array (STA) in ferromagnetic Fe3GaTe2 using a vector magnetic field manipulation technique. The orientation and ordering of STA, along with the types and density of skyrmions, are precisely controlled by modulating parameters during the manipulation. The critical roles of in-plane magnetic fields and Dzyaloshinskii-Moriya interaction in STA generation is further confirmed by micromagnetic simulation. Our findings develop a strategy for engineering large-area and highly-oriented skyrmion configurations, offering a new pathway for the future application of next-generation spintronic and information technologies.
Non-van der Waals two-dimensional (2D) materials derived from strongly bonded non-layered crystals have recently emerged as a novel and rising platform for nanoscale research. While uncovering and tuning their (opto-)electronic, catalytic, and magnetic properties has been the focus of intense research, the impact of spin-orbit coupling (SOC) onto their electronic structure has not yet been explored in detail. Studying these effects is, however, particularly relevant due to their surface cation termination and the presence of heavy elements in several representative compounds. Here, we investigate the effect of SOC onto the electronic structure of 2D AgBiO3, NaBiO3, and SbTlO3. While the first two systems show negligible band renormalization upon inclusion of relativistic effects around the band gap, SbTlO3 showcases a large SOC induced splitting (229meV) for the lowest conduction bands associated with a band inversion. Substitution of Tl with Pb forming SbPbO3 brings the band-inverted feature to the Fermi level. Analysis of topological invariants and investigation of edge states of zig-zag and armchair ribbons within the 200meV gap confirms the topological nature of the band splitting. Our work thus establishes a foundation for the systematic study of robust non-van der Waals 2D topological insulators.
Altermagnets are a new class of magnets accompanying global time-reversal symmetry breaking (TRSB) without net magnetization. The TRSB results in formation of novel altermagnetic domains. Features of altermagnetic domains, in particular their responses to external stimuli, are essentially important but yet unexplored. Here, we report visualization of bulk altermagnetic domains in MnTe based on scanning magneto-optical Kerr-effect microscopy using telecom infrared wavelength. We found two distinct TRSB domains with large Kerr rotations that do not scale with its tiny bulk magnetization. We also revealed controllability and stability of domains against magnetic or thermal perturbations. Our first observation of altermagnetic domains using a laboratory-scale simple optical technique showing their movable nature provide firm bases for future fundamental and application studies of altermagnets.
We investigate optically induced phase transitions in the two-dimensional (2D) ferroelectric (FE) material Nb2O2I4 using real-time time-dependent density functional theory (rt-TDDFT). Our results demonstrate that tailored laser pulses can activate specific coherent phonon modes. Specifically, the anharmonic atomic distortions of the A1-1 and A1-2 modes at the Γ-point facilitate the reversal of in-plane polarization. By fine-tuning laser parameters, additional phonon modes at both the Y and Γ points are excited. The resulting nonequilibrium atomic dynamics enable the formation of previously unreported ferroic phases, including three antiferroelectric (AFE) phases and one ferrielectric (FiE) phase. Notably, these optically induced phases can be reverted to the initial FE state using appropriate techniques. This controllable reversibility among multiple ferroic phases positions 2D Nb2O2I4 as a highly promising candidate for next-generation electronic storage applications.
Magnetoacoustic devices that harness the strong coupling between acoustic waves and magnons have emerged as a promising platform for energy-efficient spintronics. Laser-generated pulsed surface acoustic waves (SAWs) are particularly attractive for such applications, offering broadband frequency content up to the gigahertz (GHz) range, remote excitation without lithographic patterning, and surface localization for efficient on-chip integration. In this work, we present a comprehensive experimental study of laser-generated SAW pulses in the Fe49Rh51/MgO(001) system. A thin film of the near-equiatomic FeRh alloy serves both as an opto-acoustic transducer and as a mechanical load that modulates SAW propagation. The antiferromagnetic to ferromagnetic phase transition in FeRh, occurring slightly above room temperature, is accompanied by abrupt changes in its elastic properties, enabling controlled modification of the SAW excitation efficiency and dispersion characteristics by tuning the sample temperature and laser fluence. Using 160 fs laser pulses for excitation and time-resolved Sagnac interferometry for detection, we evaluated key SAW parameters, including amplitude, spectral content, phase and group velocities, and their in-plane anisotropy. Particular emphasis is placed on the dispersion relation and its anisotropy, which govern the coherent interaction between phonons and magnons and are determined primarily by the FeRh film.
In this systematic density functional theory study, we compare a standard gradient corrected functional (PBE) with a long-range hybrid functional (HSE06), with and without correction for the dispersion forces, relatively to their ability to correctly reproduce structural and electronic properties of different bulk 3D C3N4 phases, encompassing diamond- and graphitic-like models. Corrugation is found to provide further stabilization to the layered structures with all methods. We observe that HSE06-D3 method provides results in good agreement with experimental data and with more sophisticated G0W0 calculations. Based on that, we exploited the method to investigate the nature of the bulk triplet excitons in these C3N4 structures to evaluate the S0-T1 energy difference, the selftrapping triplet exciton energy and the photoluminescence emission energy, since this is a promising vis-light photocatalyst. Nanostructuring (0D and 2D) is another relevant aspect of these materials in practical applications, therefore we have considered the effect of single or multilayer exfoliation or space confinement in nanoparticles. Finally, we also discuss how the introduction of extrinsic dopants (e.g. S atoms) in the nanostructures modifies the atomic and electronic structure.
We present a study of the melting dynamics of a two-phase eutectic solid. In situ, thin-sample experiments using a transparent eutectic alloy and two-dimensional phase field simulations calibrated for the very same alloy are combined to assess pattern formation during directional melting in a temperature gradient. Depending on the melting velocity $V_m$ and the spacing $λ$ of the pre-solidified lamellar microstructure, an unexpectedly rich diversity of melting patterns is observed, with good agreement between experiments and simulations. We unravel the different physical mechanisms leading to this diversity, and establish the scaling behaviors of (i) the penetration of the liquid along the solid-solid interface at large $V_m$, (ii) the thickening of the primary-phase fingers at low $V_m$, and (iii) a period-doubling instability for small $λ$ values. Our study provides a fundamental basis for further investigations of eutectic melting, including additive manufacturing during which melting/solidification cycles take place.
Altermagnet-based heterojunctions have demonstrated magnetoresistive effects in experiments, however, a predictive theoretical model for non-ferromagnetic structures has remained elusive. In this work, we develop a tunneling-based spin-transport theory that explicitly incorporates the transverse-wavevector ($\bf{k}_\|$)-dependent spin polarization of an altermagnet's transport channels, enabling the prediction of giant tunneling magnetoresistance (TMR). Based on the theory, we predict that the altermagnet KV$_2$Se$_2$O can reach the extreme limit of magnetoresistance. By performing first-principles transport calculations, we verify that magnetic tunnel junctions using the metallic KV$_2$Se$_2$O as the electrodes and few-layer MgO as the spacer exhibit zero-bias magnetoresistance larger than $7.57\times10^7$\%, which is robust against the bias and thickness of the spacer. Our research provides a quantitative design principle for next-generation spin-electronic devices and establishes KV$_2$Se$_2$O/MgO/KV$_2$Se$_2$O as a leading candidate material system for room-temperature ultra-high-density non-volatile memory.
Rutile RuO$_2$ has attracted great interest recently, as its magnetic ground state remains controversial. Experimental studies have reported either nonmagnetic or altermagnetic (AM) ground states in different crystalline samples of RuO$_2$, highlighting the need for a reasonable explanation to resolve this contradiction. In this study, density functional theory calculations are performed to reveal the correlation-sensitive and strain-dependent magnetism of bulk RuO$_2$. On one hand, multiple AM phases with different magnitudes of the spin magnetic moment are identified in the Hubbard parameter space for RuO$_2$. On the other hand, when appropriate strains which significantly change the crystal cell volume are applied, the ground state of RuO$_2$ can undergo transitions between the nonmagnetic state with no spin splitting and the magnetic states with spin splitting in the band structure. These findings not only demonstrate intriguing physics in 4d-electron-correlated RuO$_2$, but also retain its potential for spintronic applications.
Mound formation on flat and miscut crystal surfaces exhibits distinct growth behaviors. While mound structures are the predominant feature on flat surfaces, miscut surfaces display a smooth transition from meandered patterns to three-dimensional mounds, depending on both internal and external conditions. We investigate this morphological evolution-from meander-like surface patterns to faceted pyramidal structures-using a vicinal Cellular Automata modeling framework. The transition is shown to be governed by the competition between the Ehrlich-Schwoebel barrier and adatom mobility on terraces. Under moderate barrier strengths and sufficiently high terrace diffusivity, the system demonstrates a reversible transition from mounded configurations to regular step meandered patterns. This reveals a complex interplay between kinetic barriers and mass transport. Our simulations cover a wide range of growth conditions, including variations in deposition flux, surface diffusion rates, temperature, and miscut angle. By applying the height-height correlation function, we calculate the correlation lengths along and across the steps and analyze their scaling behavior. These results offer insight into the continuum pathways that connect distinct classes of surface structures and provide a unified framework for describing pattern evolution across different crystal growth regimes.
Electroforming of metal-oxide-metal memristors is generally attributed to the creation of oxygen-vacancy filaments within the oxide, with noble metal electrodes such as Pt and Au remaining chemically inert. Here, we demonstrate that electroforming and subsequent operation of Pt/NbOx/Nb2O5/Pt devices can induce an unexpected and highly correlated redistribution of both oxygen and platinum. Time-of-flight secondary ion mass spectrometry reveals a filamentary pathway characterized by micrometer-scale oxygen enrichment extending from the Nb2O5 layer through NbOx and deep into the Pt top electrode. Surprisingly, this is accompanied by the formation of a Pt-rich filament penetrating the oxide stack along the same filamentary path. Finite-element and lumped-element modelling show that current-controlled negative-differential-resistance operation produces localized Joule heating and high-frequency thermal cycling, which strongly enhances oxygen migration and enables thermally assisted Pt diffusion along vacancy-rich pathways. These findings reveal a previously unrecognized metal-ion transport mechanism in NbOx memristors and highlight the critical role of post-forming electrical dynamics in determining filament chemistry, stability, and device reliability.
Mg-IV-N\textsubscript{2} compounds with IV=Si, Ge, Sn are ultra-wide band gap semiconductors with various potential electronic and optoelectronic applications. They share the \textit{Pna}2\textsubscript{1} space group crystal structure. Here we present Density Function Perturbation Theory (DFPT) calculations of the vibrational modes of these materials. We focus on the vibrational modes at the zone center to establish the relation between vibrational modes and their corresponding point-group symmetries, which determine the Raman and infrared spectra but also report the full Brillouin zone phonon dispersions and density of states. We also determine the piezoelectric tensor and the elastic compliance tensor.
As the number of theoretically predicted materials continues to grow, it becomes increasingly important to assess not only their thermodynamic stability but also their kinetic viability under realistic synthesis conditions. In this study, we investigate the hypothesis that the kinetic persistence of a metastable polymorph is related to the topography of the potential energy landscape separating it from lower energy phases. To accomplish this, we develop a new method for identifying diffusionless transformation pathways between metastable polymorphs and their ground-state counterparts and discuss the energetics of those pathways with respect to the experimental observation of each phase. This algorithm is underpinned by the recently developed Crystal Normal Form, which provides a graph representation of crystal configuration space and supplies the substrate for our pathfinding algorithm. We apply this method to the titanium dioxide system which contains the well-known anatase, rutile, and brookite phases in addition to a number of hypothetical metastable polymorphs.
\textit{Ab initio} pseudo-atomic orbital (PAO) Hamiltonians express the electronic structure of a solid in a compact, localized basis that spans the same Hilbert space as a conventional Slater--Koster tight-binding model, thereby providing an exact \textit{ab initio} representation without any loss of accuracy. Building on this correspondence, we develop an environment-dependent tight-binding (EDTB) framework in which Slater--Koster hopping integrals are augmented with bond-screening functions that capture the local coordination environment. All parameters are determined by fitting to the PAO eigenvalue spectrum across multiple atomic configurations simultaneously, which breaks the degeneracy between screening and hopping parameters and yields physically meaningful, transferable models capable of generating Hamiltonians for large systems with \textit{ab initio} precision. We demonstrate the efficiency and accuracy of the approach on four prototypical systems: bulk platinum, silicon surfaces, Si/Ge~[001] superlattices, and twisted bilayer graphene with up to $4{,}324$ atoms. The method is implemented in the \paoflow{} code and integrates seamlessly with its full post-processing suite, enabling the evaluation of a broad range of electronic, optical, and transport properties.
The large-area synthesis of high-crystalline-quality two-dimensional (2D) materials is at the core of novel material integration for semiconductor technology. This effort relies on developing fabrication and characterization techniques that can uncover the material's intrinsic properties by preserving its pristine conditions. In this article, we present an all ultra-high-vacuum cluster for the growth using molecular beam epitaxy of 2D semiconductors that are unstable under ambient conditions and optical spectroscopy using low temperature (20 K) photoluminescence and Raman scattering. The optical chamber of the setup provides micrometer scale spatial resolution and the ability to scan the entire wafer. The performance of its setup regarding spatial resolution, temperature control over a temperature range of 20-300 K using a closed-cycle cryostat and long-term preservation are demonstrated using as-grown post-transition metal monochalcogenides. Furthermore, we introduce a deconvolution-based algorithm to recover spatial information under vibration using a system-specific point-spread function. This enables in situ analysis of the structural and optoelectronic properties of as-grown materials in their pristine form, providing rich and reproducible feedback for both fundamental studies and the optimization of scalable 2D material growth toward integration in advanced devices.
Here we analyze multiple symmetry-inequivalent atomic configurations across the entire composition range of the isovalent and isostructural Mo$_x$W$_{1-x}$S$_2$ alloy using density-functional theory and Monte Carlo simulations. Our results show that although structural stability and energetics are largely composition-driven, the electronic and optical properties exhibit configuration dependence, with local atomic arrangements critically shaping band-edge splitting, valley structure, effective-mass anisotropy, and optical selection rules. In contrast to the pristine monolayers, even in the absence of spin-orbit coupling (SOC), splitting of the band edges at the $K$ point is observed across the entire composition range. In particular, while the valence-band maximum (VBM) remains largely robust, the conduction-band minimum (CBM) shows strong configuration-dependent splitting from few meV up to hundredths of meV. This behaviour leads to a non-trivial dependence of the valley energetics. Configurations with well-separated conduction bands support additional optically active transitions beyond the conventional A and B excitons in MoS$_2$ and WS$_2$ monolayers, whereas nearly degenerate cases exhibit a reduced number of allowed transitions, observed for specific configurations at $x = 1/3$ and $x = 2/3$. These results demonstrate that the number and character of optically active transitions are governed not only by composition, but also by the microscopic arrangement of atoms. Moreover, we found the hole effective masses at the VBM show configuration-dependent anisotropy, reflecting sensitivity to local symmetry breaking and implying direction-dependent transport.
Shape adaptive shape memory alloy hybrid composites (SMAHCs) are composites that incorporate shape memory alloys (SMAs) to realize shape transformation. Despite the availability of numerous analytical and finite element models for predicting the transient response of SMAHCs, many approaches exhibit limitations with respect to the thermomechanical coupling and comprehensive experimental validation. Therefore, this paper presents a coupled, multiphysics, 3D finite element approach for the simulation of a SMAHC actuator, integrating mechanical, thermal and electromagnetic solvers in the Finite Element Code ANSYS LS-DYNA. The proposed approach employs a micromechanical constitutive model implemented in ANSYS LS-DYNA, to accurately capture the complex thermomechanical phase transformation of SMAs. A key feature of the model is the ability to prescribe a defined martensitic pre-strain through a preceding simulation step, in which an initially scaled SMA wire is mechanically loaded and stretched to its nominal length. This procedure enables partial detwinning of the martensitic microstructure and provides a physically motivated initialization of the material state. Joule heating of the SMA wires, as well as varying mechanical loads and ambient temperature conditions, are explicitly considered. The simulation results are validated against experimental data and a fully coupled transient staggered scheme model to assess the predictive capability of the 3D approach. The results show good qualitative agreement, reproducing the characteristic hysteresis of actuator deflection as a function of temperature. Quantitatively, the predicted deflections are of the correct order of magnitude, although marginally outside the 95 % experimental confidence interval. Overall, a consistent trend between simulation and experiment is observed, giving rise to possibility of simulating more complex SMAHC systems.
Actively tunable photonic devices are vital for next-generation optoelectronics requiring rapid switching and high bandwidth. Although organic optoelectronic devices have found wide application, their use as optical modulators has been limited by low absorption in the critical near-infrared (NIR) region, slow response time, and weak nonlinearities. To address these limitations, we developed a scheme based on intermediate exciton-photon coupling in a NIR absorbing squaraine-dye based photonic structure. Using energy-momentum resolved pump-probe spectroscopy, we show that the sign and magnitude of the optical response of our system depends strongly on the energy detuning between the excitonic and photonic modes. These data are analyzed using temporal coupled-mode theory to show that near resonance, a distinct energy exchange process emerges in the cross-over regime between strong and weak light-matter coupling. This effect enables dynamical control over the photoinduced response, providing a pathway for broadband optical signal modulation extending into the NIR spectral region.
Near-infrared (NIR) emitters operating in the 650-900 nm range are highly attractive for imaging and sensing in turbid media; however, cadmium-free InP-based quantum dots (QDs) often suffer from limited brightness due to nonradiative pathways and inefficient photon outcoupling. In particular, heterostructured InP QDs can exhibit band alignments that induce partial spatial separation of charge carriers, leading to reduced electron-hole wavefunction overlap. This modifies intrinsic recombination dynamics and enhances the sensitivity of their emission to the surrounding photonic environment. Here, we investigate silver toroidal plasmonic nanoantenna dimers (Ag TPNDs) through finite-difference-time-domain (FDTD) simulations as a geometry-tunable platform for enhancing NIR emission of heterostructured InP-based QDs. The coupled toroidal geometry supports strongly confined bonding modes that generate intense nanogap hotspots, while its resonance can be systematically tuned through the toroid aspect ratio. By spectrally aligning the antenna response with QD emission bands (675-845 nm), we achieve large Purcell enhancements together with high quantum efficiencies, demonstrating efficient conversion of enhanced decay rates into radiative emission. We further show that nanometer-scale variations in emitter-antenna separation strongly modulate the radiative rates and spectral response. These results establish toroidal plasmonic nanodimers as a topology-driven platform for controlling emission in NIR quantum emitters and for advancing NIR nanophotonic applications.
Understanding the thermodynamics of radiation and the quantum-mechanical interactions between light and matter is important both for theoretical purposes and for technological advances, such as determining the limits of key processes like light-to-usable-energy conversion efficiencies. In this report, we discuss the physics of these two aspects, considering spontaneous emission as a pathway, and highlight the limitations of such descriptions in assessing energy-harvesting efficiency. In view of these limitations, we adopt a simplified approach to evaluate the free energy of radiation, providing a framework to assess various aspects of light-to-usable-energy conversion efficiencies. Our approach allows a theoretical estimate of the thermodynamic maximum limit for light-to-usable-energy conversion, which is approximately 74%. We validate this free energy estimate by modeling and accurately reproducing the Shockley-Queisser limit (~ 33%), which imposes a practical constraint on solar-to-usable-energy conversion efficiency. Beyond free-energy considerations, our model incorporates various processes, such as spontaneous emission, nonradiative thermal losses, and photon upconversion, allowing us to evaluate their roles. The model further suggests that, under certain conditions, the maximum conversion efficiency can reach approximately 48%, for example with multijunction solar cells or via photon upconversion. These findings further suggest that the true thermodynamic limit for light-to-usable-energy conversion may be much higher (approximately 74%). However, accurately estimating this limit requires a more complete understanding of the thermodynamics of light, light-matter interactions, and the connection between them.
Active colloidal microparticles require reliable actuation to sustain directed motion. Light-based propulsion is particularly attractive as it provides persistent energy supply and enables direct spatiotemporal control. Here, we introduce 3D-printable particles with symmetry-broken refractive index profiles (SBRIP particles) that achieve propulsion through direct momentum transfer from asymmetric light refraction. Internal refractive-index gradients provide optical symmetry breaking independent of external shape, fundamentally decoupling propulsion from particle geometry. Geometrically symmetry-broken particles with a homogeneous refractive index are another special case, where propulsion originates from refractive contrast at the boundary instead of within the particle. Unlike conventional systems relying on absorption or reflection, this transparency-based mechanism minimizes heating and mitigates shadowing in bulk suspensions. We present a theoretical framework for refractive propulsion as well as numerical simulations of the SBRIP particles using raytracing and the finite volume method. This is complemented by experiments, validating the momentum transfer mechanism using particles with geometric symmetry breaking. The high transparency of our particles ensures deep light penetration, enabling the realization of volumetric active matter. This opens pathways toward adaptive nonlinear optical materials where light-driven particle reorganization modulates the local refractive index, establishing a dynamic feedback loop between the optical field and the material structure.