Entropy production in quasi-isentropic compression (QIC) is critically important for understanding the properties of materials under extremeconditions. However, the origin and accurate quantification of entropy in thi...Entropy production in quasi-isentropic compression (QIC) is critically important for understanding the properties of materials under extremeconditions. However, the origin and accurate quantification of entropy in this situation remain long-standing challenges. In this work, a framework is established for the quantification of entropy production and partition, and their relation to microstructural change in QIC. Cu50Zr50is taken as a model material, and its compression is simulated by molecular dynamics. On the basis of atomistic simulation-informed physicalproperties and free energy, the thermodynamic path is recovered, and the entropy production and its relation to microstructural change aresuccessfully quantified by the proposed framework. Contrary to intuition, entropy production during QIC of metallic glasses is relativelyinsensitive to the strain rate ˙γ when ˙γ ranges from 7.5 × 10^(8) to 2 × 10^(9)/s, which are values reachable in QIC experiments, with a magnitudeof the order of 10^(−2)kB/atom per GPa. However, when ˙γ is extremely high (>2 × 10^(9)/s), a notable increase in entropy production rate with˙γ is observed. The Taylor–Quinney factor is found to vary with strain but not with strain rate in the simulated regime. It is demonstrated thatentropy production is dominated by the configurational part, compared with the vibrational part. In the rate-insensitive regime, the increase inconfigurational entropy exhibits a linear relation to the Shannon-entropic quantification of microstructural change, and a stretched exponential relation to the Taylor–Quinney factor. The quantification of entropy is expected to provide thermodynamic insights into the fundamentalrelation between microstructure evolution and plastic dissipation.展开更多
For a long time,there have been huge discrepancies between different models and experiments concerning the liquid-liquid phase transition(LLPT)in dense hydrogen.We present the results of extensive calculations of the ...For a long time,there have been huge discrepancies between different models and experiments concerning the liquid-liquid phase transition(LLPT)in dense hydrogen.We present the results of extensive calculations of the LLPT in dense hydrogen using the most expensive first-principle path-integral molecular dynamics simulations available.The nonlocal density functional rVV10 and the hybrid functional PBEO are used to improve the description of the electronic structure of hydrogen.Of all the density functional theory calculations available,we report the most consistent results through quantum Monte Carlo simulations and coupled electron-ion Monte Carlo simulations of the LLPT in dense hydrogen.The critical point of the first-order LLPT is estimated to be above 2000 K according to the equation of state.Moreover,the metallization pressure obtained from the jump of dc electrical conductivity almost coincides with the plateau of equation of state.展开更多
Accurate knowledge of the equation of state(EOS)of deuterium–tritium(DT)mixtures is critically important for inertial confinement fusion(ICF).Although the study of EOS is an old topic,there is a longstanding lack of ...Accurate knowledge of the equation of state(EOS)of deuterium–tritium(DT)mixtures is critically important for inertial confinement fusion(ICF).Although the study of EOS is an old topic,there is a longstanding lack of global accurate EOS data for DT within a unified theoretical framework.DT fuel goes through very wide ranges of density and temperature from a cold condensed state to a hot dense plasma where ions are in a moderately or even strongly coupled state and electrons are in a partially or strongly degenerate state.The biggest challenge faced when using first-principles methods for obtaining accurate EOS data for DT fuel is the treatment of electron–ion interactions and the extremely high computational cost at high temperatures.In the present work,we perform extensive state-of-the-art ab initio quantum Langevin molecular dynamics simulations to obtain EOS data for DT mixtures at densities from 0.1 g/cm3 to 2000 g/cm3 and temperatures from 500 K to 2000 eV,which are relevant to ICF processes.Comparisons with average-atom molecular dynamics and orbital-free molecular dynamics simulations show that the ionic strong-coupling effect is important for determining the whole-range EOS.This work can supply accurate EOS data forDTmixtures within a unified ab initio framework,as well as providing a benchmark for various semiclassical methods.展开更多
Lattice thermal conductivity(κlat)of MgSiO_(3) perovskite and post-perovskite is an important parameter for the thermal dynamics in the Earth.Here,we develop a deep potential of density functional theory quality unde...Lattice thermal conductivity(κlat)of MgSiO_(3) perovskite and post-perovskite is an important parameter for the thermal dynamics in the Earth.Here,we develop a deep potential of density functional theory quality under entire thermodynamic conditions in the lower mantle,and calculate theκlatby the Green-Kubo relation.Deep potential molecular dynamics captures full-order anharmonicity and considers ill-defined phonons in low-κlatmaterials ignored in the phonon gas model.Theκlatshows negative temperature dependence and positive linear pressure dependence.Interestingly,theκlatundergos an increase at the phase boundary from perovskite to post-perovskite.We demonstrate that,along the geotherm,theκlatincreases by 18.2% at the phase boundary.Our results would be helpful for evaluating Earth’s thermal dynamics and improving the Earth model.展开更多
Superionic ices with highly mobile protons within stable oxygen sub-lattices occupy an important proportion of the phase diagram of ice and widely exist in the interior of icy giants and throughout the Universe.Unders...Superionic ices with highly mobile protons within stable oxygen sub-lattices occupy an important proportion of the phase diagram of ice and widely exist in the interior of icy giants and throughout the Universe.Understanding the thermal transport in superionic ice is vital for the thermal evolution of icy planets.However,it is highly challenging due to the extreme thermodynamic conditions and dynamical nature of protons,beyond the capability of the traditional lattice dynamics and empirical potential molecular dynamics approaches.By utilizing the deep potential molecular dynamics approach,we investigate the thermal conductivity of ice-Ⅶ and superionic ice-Ⅶ’’ along the isobar of P = 30 GPa.A non-monotonic trend of thermal conductivity with elevated temperature is observed.Through heat flux decomposition and trajectory-based spectra analysis,we show that the thermally activated proton diffusion in ice-Ⅶ and superionic ice-Ⅶ′′contribute significantly to heat convection,while the broadening in vibrational energy peaks and significant softening of transverse acoustic branches lead to a reduction in heat conduction.The competition between proton diffusion and phonon scattering results in anomalous thermal transport across the superionic transition in ice.This work unravels the important role of proton diffusion in the thermal transport of high-pressure ice.Our approach provides new insights into modeling the thermal transport and atomistic dynamics in superionic materials.展开更多
The equation of states,diffusions,and viscosities of strongly coupled Fe at 80 and 240 eV with densities from 1.6 to 40 g/cm^(3) are studied by orbital-free molecular dynamics,classical molecular dynamics with a corre...The equation of states,diffusions,and viscosities of strongly coupled Fe at 80 and 240 eV with densities from 1.6 to 40 g/cm^(3) are studied by orbital-free molecular dynamics,classical molecular dynamics with a corrected Yukawa potential and compared with the results from average atom model.A new local pseudopotential is generated for orbital free calculations.For low densities,the Yukawa model captures the correct ionic interaction behavior around the first peak of the radial distribution function(RDF),thus it gives correct RDFs and transport coefficients.For higher densities,the scaled transformation of the Yukawa potential or adding a short range repulsion part to the Yukawa potential can give correct RDFs and transport coefficients.The corrected potentials are further validated by the force matching method.展开更多
We consider a steady-state(but transient)situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature T_(e),ions at T_(i),and T_(e)≠T_(i).Such states are achievable ...We consider a steady-state(but transient)situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature T_(e),ions at T_(i),and T_(e)≠T_(i).Such states are achievable by pump–probe experiments.For warm dense hydrogen in such a twotemperature situation,we investigate nuclear quantum effects(NQEs)on structure and thermodynamic properties,thereby delineating the limitations of ordinary ab initio molecular dynamics.We use path integral molecular dynamics(PIMD)simulations driven by orbital-free density functional theory(OFDFT)calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence.We calibrate the OFDFT calculations against conventional(explicit orbitals)Kohn–Sham DFT.We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30,the ionic radial distribution function is meaningfully affected by the inclusion of NQEs.Moreover,NQEs induce a substantial increase in both the ionic and electronic pressures.This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen.For Te>20 kK,increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered.In addition,we confirm that compared with thermostatted ring-polymer molecular dynamics,the primitive PIMD algorithm overestimates electronic pressures,a consequence of the overly localized ionic description from the primitive scheme.展开更多
High-Zmaterials exhibit a broad range of variation of the charge state in the hot dense regime,and so ionic structures becomecomplexwith increasing density and temperature owing to ionization.Taking high-Z uranium as ...High-Zmaterials exhibit a broad range of variation of the charge state in the hot dense regime,and so ionic structures becomecomplexwith increasing density and temperature owing to ionization.Taking high-Z uranium as example,we study its electronic and ionic structures in the hot dense regime by combining an average-atommodelwith the hypernetted chain approximation.The electronic structure is described by solving theDirac equation,taking account of relativistic effects,including broadening of the energy levels,and the effect of other ions via correlation functions.On the basis of the electronic distribution around a nucleus,the ion pair potential is constructed using the modified Gordon–Kim model in the frame of temperaturedependent density functional theory.Because of the presence of ion–ion strong coupling,the bridge function is included in the hypernetted chain approximation,whichis usedto calculate the correlation functions.To take account of the influenceon transportpropertiesof the strong correlation of electrons with highly charged ions,we perform both classical and Langevin molecular dynamics simulations to determine ion self-diffusion coefficients and the shear viscosity,using theGreen–Kubo relation and an ion–ion pair potential with good convergence.We show that the influence of electron–ion collisions on transport properties becomes more important as the free electron density increases owing to thermal ionization.展开更多
Two-dimensional(2D)van der Waals materials have attracted tremendous attention due to their versatile physical properties and flexible manipulation approaches.Among the various types of van der Waals materials,α-In_(...Two-dimensional(2D)van der Waals materials have attracted tremendous attention due to their versatile physical properties and flexible manipulation approaches.Among the various types of van der Waals materials,α-In_(2)Se_(3)is remarkable for its intrinsic 2D ferroelectricity and high-performance opto-electronic properties.However,the study of theα-In_(2)Se_(3)system in terahertz(THz)radiation is scarce,although it is promising for electrically controlled THz field manipulation.We investigate theα-In_(2)Se_(3)in different thicknesses and report that the THz generation efficiency induced by femtosecond laser pulses can be largely improved by reducing the thickness from the bulk.Furthermore,we reveal the surge current in thin film coupled with THz emission exhibits a different Auger recombination mode,which is helpful in understanding the mechanism and provides insights into the design of 2D highly efficient THz devices.展开更多
The coupling of excited states and ionic dynamics is the basic and challenging point for the materials response at extreme conditions.In the laboratory,the intense laser produces transient nature and complexity with h...The coupling of excited states and ionic dynamics is the basic and challenging point for the materials response at extreme conditions.In the laboratory,the intense laser produces transient nature and complexity with highly nonequilibrium states,making it extremely difficult and interesting for both experimental measurements and theoretical methods.With the inclusion of laser-excited states,we extend an ab initio method into the direct simulations of whole laser-driven microscopic dynamics from solid to liquid.We construct the framework of combining the electron-temperature-dependent deep neural-network potential energy surface with a hybrid atomistic-continuum approach,controlling non-adiabatic energy exchange and atomistic dynamics,which enables consistent interpretation of experimental data.By large-scale ab initio simulations,we demonstrate that the nonthermal effects introduced by hot electrons play a dominant role in modulating the lattice dynamics,thermodynamic pathway,and structural transformation.We highlight that the present work provides a path to realistic computational studies of laser-driven processes,thus bridging the gap between experiments and simulations.展开更多
The internal solar structure predicted by the standard solar model disagrees with the helioseismic observations even by utilizing the most updated physical inputs, such as the opacity and element abundances. By increa...The internal solar structure predicted by the standard solar model disagrees with the helioseismic observations even by utilizing the most updated physical inputs, such as the opacity and element abundances. By increasing the Rosseland mean, the decade-old open problem of the missing opacity can be resolved. Herein, we propose that the continuum electrons in the radiative processes lose phases and coherence as matter waves, giving rise to a phenomenon of transient spatial localization. It not only enhances the continuum opacity but also increases the line widths of the bound-bound transitions. We demonstrate our theoretical formulation by investigating the opacity of solar mixtures in the interior. The Rosseland mean demonstrates an increase of 10%-26% in the range of 0.3 R⊙-0.75 R⊙. The results are compared with the recent experimental data and the existing theoretical models. Our findings provide novel clues to the open problem of the missing opacity in the solar interior and new insight on the radiative opacity in the hot dense-plasma regime.展开更多
基金supported by the NSAF under Grant No.U1830206,the National Key R&D Program of China under Grant No.2017YFA0403200the National Natural Science Foundation of China under Grant Nos.11874424 and 12104507the Science and Technology Innovation Program of Hunan Province under Grant No.2021RC4026.
文摘Entropy production in quasi-isentropic compression (QIC) is critically important for understanding the properties of materials under extremeconditions. However, the origin and accurate quantification of entropy in this situation remain long-standing challenges. In this work, a framework is established for the quantification of entropy production and partition, and their relation to microstructural change in QIC. Cu50Zr50is taken as a model material, and its compression is simulated by molecular dynamics. On the basis of atomistic simulation-informed physicalproperties and free energy, the thermodynamic path is recovered, and the entropy production and its relation to microstructural change aresuccessfully quantified by the proposed framework. Contrary to intuition, entropy production during QIC of metallic glasses is relativelyinsensitive to the strain rate ˙γ when ˙γ ranges from 7.5 × 10^(8) to 2 × 10^(9)/s, which are values reachable in QIC experiments, with a magnitudeof the order of 10^(−2)kB/atom per GPa. However, when ˙γ is extremely high (>2 × 10^(9)/s), a notable increase in entropy production rate with˙γ is observed. The Taylor–Quinney factor is found to vary with strain but not with strain rate in the simulated regime. It is demonstrated thatentropy production is dominated by the configurational part, compared with the vibrational part. In the rate-insensitive regime, the increase inconfigurational entropy exhibits a linear relation to the Shannon-entropic quantification of microstructural change, and a stretched exponential relation to the Taylor–Quinney factor. The quantification of entropy is expected to provide thermodynamic insights into the fundamentalrelation between microstructure evolution and plastic dissipation.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11774429,11874424 and U1830206the Science Challenge Project under Grant No TZ2016001+2 种基金the National Key R&D Program of China under Grant No 2017YFA0403200the Science and Technology Project of Hunan Province under Grant No 2017RS3038the Advanced Research Foundation of National University of Defense Technology under Grant No JQ14-02-01
文摘For a long time,there have been huge discrepancies between different models and experiments concerning the liquid-liquid phase transition(LLPT)in dense hydrogen.We present the results of extensive calculations of the LLPT in dense hydrogen using the most expensive first-principle path-integral molecular dynamics simulations available.The nonlocal density functional rVV10 and the hybrid functional PBEO are used to improve the description of the electronic structure of hydrogen.Of all the density functional theory calculations available,we report the most consistent results through quantum Monte Carlo simulations and coupled electron-ion Monte Carlo simulations of the LLPT in dense hydrogen.The critical point of the first-order LLPT is estimated to be above 2000 K according to the equation of state.Moreover,the metallization pressure obtained from the jump of dc electrical conductivity almost coincides with the plateau of equation of state.
基金This work was supported by the Science Challenge Project under Grant No.TZ2016001the National Key R&D Program of China under Grant No.2017YFA0403200+1 种基金the National Natural Science Foundation of China under Grant Nos.11874424 and 11774429the NSAF under Grant No.U1830206.All calculations were carried out at the Research Center of Supercomputing Applications at NUDT.
文摘Accurate knowledge of the equation of state(EOS)of deuterium–tritium(DT)mixtures is critically important for inertial confinement fusion(ICF).Although the study of EOS is an old topic,there is a longstanding lack of global accurate EOS data for DT within a unified theoretical framework.DT fuel goes through very wide ranges of density and temperature from a cold condensed state to a hot dense plasma where ions are in a moderately or even strongly coupled state and electrons are in a partially or strongly degenerate state.The biggest challenge faced when using first-principles methods for obtaining accurate EOS data for DT fuel is the treatment of electron–ion interactions and the extremely high computational cost at high temperatures.In the present work,we perform extensive state-of-the-art ab initio quantum Langevin molecular dynamics simulations to obtain EOS data for DT mixtures at densities from 0.1 g/cm3 to 2000 g/cm3 and temperatures from 500 K to 2000 eV,which are relevant to ICF processes.Comparisons with average-atom molecular dynamics and orbital-free molecular dynamics simulations show that the ionic strong-coupling effect is important for determining the whole-range EOS.This work can supply accurate EOS data forDTmixtures within a unified ab initio framework,as well as providing a benchmark for various semiclassical methods.
基金supported by the National Natural Science Foundation of China(Grant No.U1830206)the National Key R&D Program of China(Grant No.2017YFA0403200)+1 种基金the National Natural Science Foundation of China(Grant Nos.11874424,11904401,11974423,and 12104507)the Science and Technology Innovation Program of Hunan Province(Grant No.2021RC4026)。
文摘Lattice thermal conductivity(κlat)of MgSiO_(3) perovskite and post-perovskite is an important parameter for the thermal dynamics in the Earth.Here,we develop a deep potential of density functional theory quality under entire thermodynamic conditions in the lower mantle,and calculate theκlatby the Green-Kubo relation.Deep potential molecular dynamics captures full-order anharmonicity and considers ill-defined phonons in low-κlatmaterials ignored in the phonon gas model.Theκlatshows negative temperature dependence and positive linear pressure dependence.Interestingly,theκlatundergos an increase at the phase boundary from perovskite to post-perovskite.We demonstrate that,along the geotherm,theκlatincreases by 18.2% at the phase boundary.Our results would be helpful for evaluating Earth’s thermal dynamics and improving the Earth model.
基金supported by the National Natural Science Foundation of China(Grant Nos.11874424,12122103,and 12304307)the Science and Technology Innovation Program of Hunan Province(Grant No.2021RC4026)。
文摘Superionic ices with highly mobile protons within stable oxygen sub-lattices occupy an important proportion of the phase diagram of ice and widely exist in the interior of icy giants and throughout the Universe.Understanding the thermal transport in superionic ice is vital for the thermal evolution of icy planets.However,it is highly challenging due to the extreme thermodynamic conditions and dynamical nature of protons,beyond the capability of the traditional lattice dynamics and empirical potential molecular dynamics approaches.By utilizing the deep potential molecular dynamics approach,we investigate the thermal conductivity of ice-Ⅶ and superionic ice-Ⅶ’’ along the isobar of P = 30 GPa.A non-monotonic trend of thermal conductivity with elevated temperature is observed.Through heat flux decomposition and trajectory-based spectra analysis,we show that the thermally activated proton diffusion in ice-Ⅶ and superionic ice-Ⅶ′′contribute significantly to heat convection,while the broadening in vibrational energy peaks and significant softening of transverse acoustic branches lead to a reduction in heat conduction.The competition between proton diffusion and phonon scattering results in anomalous thermal transport across the superionic transition in ice.This work unravels the important role of proton diffusion in the thermal transport of high-pressure ice.Our approach provides new insights into modeling the thermal transport and atomistic dynamics in superionic materials.
基金This work is supported by the National Basic Research Program of China(973 Program)under grant no.2013CB922203the National NSFC under grant Nos.11422432 and 11774429+1 种基金Science Challenge Project under grant no.JCKY2016212A501the Advanced Research Foundation of National University of Defense Technology under grant no.JQ14-02-01.
文摘The equation of states,diffusions,and viscosities of strongly coupled Fe at 80 and 240 eV with densities from 1.6 to 40 g/cm^(3) are studied by orbital-free molecular dynamics,classical molecular dynamics with a corrected Yukawa potential and compared with the results from average atom model.A new local pseudopotential is generated for orbital free calculations.For low densities,the Yukawa model captures the correct ionic interaction behavior around the first peak of the radial distribution function(RDF),thus it gives correct RDFs and transport coefficients.For higher densities,the scaled transformation of the Yukawa potential or adding a short range repulsion part to the Yukawa potential can give correct RDFs and transport coefficients.The corrected potentials are further validated by the force matching method.
基金The majority of this work was done while D.K.was a visitor at the University of Florida.He was supported by the Science Challenge Project of China under Grant No.TZ2016001the NSFC under Grant No.11874424+3 种基金the National Key R&D Program of China under Grant No.2017YFA0403200He also acknowledges support by the China Scholarship Council.K.L.(for the majority of the work done while at the University of Florida)S.B.T.were supported by U.S.Department of Energy Grant No.DE-SC0002139Most of the computations were performed on the HiPerGator-II system at the University of Florida.
文摘We consider a steady-state(but transient)situation in which a warm dense aggregate is a two-temperature system with equilibrium electrons at temperature T_(e),ions at T_(i),and T_(e)≠T_(i).Such states are achievable by pump–probe experiments.For warm dense hydrogen in such a twotemperature situation,we investigate nuclear quantum effects(NQEs)on structure and thermodynamic properties,thereby delineating the limitations of ordinary ab initio molecular dynamics.We use path integral molecular dynamics(PIMD)simulations driven by orbital-free density functional theory(OFDFT)calculations with state-of-the-art noninteracting free-energy and exchange-correlation functionals for the explicit temperature dependence.We calibrate the OFDFT calculations against conventional(explicit orbitals)Kohn–Sham DFT.We find that when the ratio of the ionic thermal de Broglie wavelength to the mean interionic distance is larger than about 0.30,the ionic radial distribution function is meaningfully affected by the inclusion of NQEs.Moreover,NQEs induce a substantial increase in both the ionic and electronic pressures.This confirms the importance of NQEs for highly accurate equation-of-state data on highly driven hydrogen.For Te>20 kK,increasing Te in the warm dense hydrogen has slight effects on the ionic radial distribution function and equation of state in the range of densities considered.In addition,we confirm that compared with thermostatted ring-polymer molecular dynamics,the primitive PIMD algorithm overestimates electronic pressures,a consequence of the overly localized ionic description from the primitive scheme.
文摘High-Zmaterials exhibit a broad range of variation of the charge state in the hot dense regime,and so ionic structures becomecomplexwith increasing density and temperature owing to ionization.Taking high-Z uranium as example,we study its electronic and ionic structures in the hot dense regime by combining an average-atommodelwith the hypernetted chain approximation.The electronic structure is described by solving theDirac equation,taking account of relativistic effects,including broadening of the energy levels,and the effect of other ions via correlation functions.On the basis of the electronic distribution around a nucleus,the ion pair potential is constructed using the modified Gordon–Kim model in the frame of temperaturedependent density functional theory.Because of the presence of ion–ion strong coupling,the bridge function is included in the hypernetted chain approximation,whichis usedto calculate the correlation functions.To take account of the influenceon transportpropertiesof the strong correlation of electrons with highly charged ions,we perform both classical and Langevin molecular dynamics simulations to determine ion self-diffusion coefficients and the shear viscosity,using theGreen–Kubo relation and an ion–ion pair potential with good convergence.We show that the influence of electron–ion collisions on transport properties becomes more important as the free electron density increases owing to thermal ionization.
基金supported by the National Key R&D Program of China(No.2017YFA0403200)the NSAF(No.U1830206)+1 种基金the Science and Technology Innovation Program of Hunan Province(No.2021RC4026)the National Natural Science Foundation of China(NSFC)(No.62005058)。
文摘Two-dimensional(2D)van der Waals materials have attracted tremendous attention due to their versatile physical properties and flexible manipulation approaches.Among the various types of van der Waals materials,α-In_(2)Se_(3)is remarkable for its intrinsic 2D ferroelectricity and high-performance opto-electronic properties.However,the study of theα-In_(2)Se_(3)system in terahertz(THz)radiation is scarce,although it is promising for electrically controlled THz field manipulation.We investigate theα-In_(2)Se_(3)in different thicknesses and report that the THz generation efficiency induced by femtosecond laser pulses can be largely improved by reducing the thickness from the bulk.Furthermore,we reveal the surge current in thin film coupled with THz emission exhibits a different Auger recombination mode,which is helpful in understanding the mechanism and provides insights into the design of 2D highly efficient THz devices.
基金supported by the National Natural Science Foundation of China under Grant Nos.11874424,11904401,12104507,12304307the Science and Technology Innovation Program of Hunan Province under Grant No.2021RC4026.
文摘The coupling of excited states and ionic dynamics is the basic and challenging point for the materials response at extreme conditions.In the laboratory,the intense laser produces transient nature and complexity with highly nonequilibrium states,making it extremely difficult and interesting for both experimental measurements and theoretical methods.With the inclusion of laser-excited states,we extend an ab initio method into the direct simulations of whole laser-driven microscopic dynamics from solid to liquid.We construct the framework of combining the electron-temperature-dependent deep neural-network potential energy surface with a hybrid atomistic-continuum approach,controlling non-adiabatic energy exchange and atomistic dynamics,which enables consistent interpretation of experimental data.By large-scale ab initio simulations,we demonstrate that the nonthermal effects introduced by hot electrons play a dominant role in modulating the lattice dynamics,thermodynamic pathway,and structural transformation.We highlight that the present work provides a path to realistic computational studies of laser-driven processes,thus bridging the gap between experiments and simulations.
基金supported by the Science Challenge Project(Grant No.TZ2018005)the National Key R&D Program of China(Grant Nos.2019YFA0307700,and 2017YFA0403202)the National Natural Science Foundation of China(Grant Nos.12174343,and 11774322)。
文摘The internal solar structure predicted by the standard solar model disagrees with the helioseismic observations even by utilizing the most updated physical inputs, such as the opacity and element abundances. By increasing the Rosseland mean, the decade-old open problem of the missing opacity can be resolved. Herein, we propose that the continuum electrons in the radiative processes lose phases and coherence as matter waves, giving rise to a phenomenon of transient spatial localization. It not only enhances the continuum opacity but also increases the line widths of the bound-bound transitions. We demonstrate our theoretical formulation by investigating the opacity of solar mixtures in the interior. The Rosseland mean demonstrates an increase of 10%-26% in the range of 0.3 R⊙-0.75 R⊙. The results are compared with the recent experimental data and the existing theoretical models. Our findings provide novel clues to the open problem of the missing opacity in the solar interior and new insight on the radiative opacity in the hot dense-plasma regime.