We perform a three-dimensional numerical simulation based on a one-step chemical reaction model to investigate changes in the mode of H2-Air detonation wave propagation from rotating detonation wave （RDW） mode to st...We perform a three-dimensional numerical simulation based on a one-step chemical reaction model to investigate changes in the mode of H2-Air detonation wave propagation from rotating detonation wave （RDW） mode to standing detonation wave mode. The physical characteristics of an RDW with injection velocity of 500 m/s are analyzed to investigate the physical mechanisms involved. We find that with increasing injection velocity, the detonation wave gradually changes from perpendicular to the head wall to parallel to the head wall. When the injection velocity exceeds the Chapman-Jouguet velocity VCJ （about 1984 m/s）, the detonation wave changes orientation to become perpendicular to the fuel injection direction, and the rotating mode changes accordingly to a standing mode. Finally, the plane detonation characteristic triple-wave structures can be found from the standing mode.展开更多
In traditional finite-temperature Kohn–Sham density functional theory(KSDFT),the partial occupation of a large number of high-energy KS eigenstates restricts the use of first-principles molecular dynamics methods at ...In traditional finite-temperature Kohn–Sham density functional theory(KSDFT),the partial occupation of a large number of high-energy KS eigenstates restricts the use of first-principles molecular dynamics methods at extremely high temperatures.However,stochastic density functional theory(SDFT)can overcome this limitation.Recently,SDFT and the related mixed stochastic–deterministic density functional theory,based on a plane-wave basis set,have been implemented in the first-principles electronic structure software ABACUS[Q.Liu and M.Chen,Phys.Rev.B 106,125132(2022)].In this study,we combine SDFT with the Born–Oppenheimer molecular dynamics method to investigate systems with temperatures ranging from a few tens of eV to 1000 eV.Importantly,we train machine-learning-based interatomic models using the SDFT data and employ these deep potential models to simulate large-scale systems with long trajectories.Subsequently,we compute and analyze the structural properties,dynamic properties,and transport coefficients of warm dense matter.展开更多
To fill the gap between accurate(and expensive)ab initio calculations and efficient atomistic simulations based on empirical interatomic potentials,a new class of descriptions of atomic interactions has emerged and be...To fill the gap between accurate(and expensive)ab initio calculations and efficient atomistic simulations based on empirical interatomic potentials,a new class of descriptions of atomic interactions has emerged and been widely applied;i.e.machine learning potentials(MLPs).One recently developed type of MLP is the deep potential(DP)method.In this review,we provide an introduction to DP methods in computational materials science.The theory underlying the DP method is presented along with a step-by-step introduction to their development and use.We also review materials applications of DPs in a wide range of materials systems.The DP Library provides a platform for the development of DPs and a database of extant DPs.We discuss the accuracy and efficiency of DPs compared with ab initio methods and empirical potentials.展开更多
In this study,we theoretically investigate the feasibility of using laser-wakefield accelerated(LWFA)electrons for the photonuclear measurement of nuclear isomers according to the characteristics of the electrons obta...In this study,we theoretically investigate the feasibility of using laser-wakefield accelerated(LWFA)electrons for the photonuclear measurement of nuclear isomers according to the characteristics of the electrons obtained from LWFA experiments conducted at the Compact Laser-Plasma Accelerator(CLAPA)laboratory.The experiments at the CLAPA show that a stable electron beam with an energy of 78-135 MeV and a charge of 300-600 pC can be obtained.The bremsstrahlung spectra were simulated using Geant4,which suggests that a bremsstrahlung source with a peak intensity of 10^(19) photons/s can be generated.Theoretical calculations of isomer production cross sections from the photonuclear reactions on six target nuclei,^(197)Au,^(180)Hf,^(159)Tb,^(115)In,^(103)Rh,and ^(90)Zr,were performed and compared with the available experimental data in EXFOR,which suggest that further experiments are required for a series of photonuclear reaction channels.Flux-averaged cross sections and isomer ratios(IR)resulting from such bremsstrahlung sources are theoretically deduced.The results suggest that IR measurements can be used to constrain nuclear components,such asγstrength function and optical model potential.In addition,the detection of the decay characteristics was evaluated with Geant4 simulations.The use of the LWFA electron beam and its bremsstrahlung for photonuclear studies involving nuclear isomers is anticipated.展开更多
Ultrashort and powerful laser interactions with a target generate intense wideband electromagnetic pulses(EMPs).In this study,we report EMPs generated by the interactions between petawatt(30 fs,1.4×10^(20) W/cm^(...Ultrashort and powerful laser interactions with a target generate intense wideband electromagnetic pulses(EMPs).In this study,we report EMPs generated by the interactions between petawatt(30 fs,1.4×10^(20) W/cm^(2))femtosecond(fs)lasers with metal flat,plastic flat,and plastic nanowire-array(NWA)targets.Detailed analyses are conducted on the EMPs in terms of their spatial distribution,time and frequency domains,radiation energy,and protection.The results indicate that EMPs from metal targets exhibit larger amplitudes at varying angles than those generated by other types of targets and are enhanced significantly for NWA targets.Using a plastic target holder and increasing the laser focal spot can significantly decrease the radiation energy of the EMPs.Moreover,the composite shielding materials indicate an effective shielding effect against EMPs.The simulation results show that the NWA targets exert a collimating effect on thermal electrons,which directly affects the distribution of EMPs.This study provides guidance for regulating EMPs by controlling the laser focal spot,target parameters,and target rod material and is beneficial for electromagnetic-shielding design.展开更多
Giant electromagnetic pulses(EMPs) induced by high-power laser irradiating solid targets interfere with various experimental diagnoses and even damage equipment,so unveiling the evolution of EMPs inside the laser cham...Giant electromagnetic pulses(EMPs) induced by high-power laser irradiating solid targets interfere with various experimental diagnoses and even damage equipment,so unveiling the evolution of EMPs inside the laser chamber is crucial for designing effective EMP shielding.In this work,the transmission characteristics of EMPs as a function of distances from the target chamber center(TCC) are studied using B-dot probes.The mean EMP amplitude generated by picosecond laser-target interaction reaches 561 kV m^(-1),357 kV m^(-1),395 kV m^(-1),and 341 kV m^(-1)at 0.32 m,0.53 m,0.76 m,and 1 m from TCC,which decreases dramatically from 0.32 m to 0.53 m.However,it shows a fluctuation from 0.53 m to 1 m.The temporal features of EMPs indicate that time-domain EMP signals near the target chamber wall have a wider full width at half maximum compared to that close to TCC,mainly due to the echo oscillation of electromagnetic waves inside the target chamber based on simulation and experimentation.The conclusions of this study will provide a new approach to mitigate strong electromagnetic pulses by decreasing the echo oscillation of electromagnetic waves inside the target chamber during laser coupling with targets.展开更多
The newly built Compact Laser Plasma Accelerator-Therapy facility at Peking University will deliver 60 J/1 Hz laser pulses with 30 fs duration.Driven by this petawatt laser facility,proton beams with energy up to 200 ...The newly built Compact Laser Plasma Accelerator-Therapy facility at Peking University will deliver 60 J/1 Hz laser pulses with 30 fs duration.Driven by this petawatt laser facility,proton beams with energy up to 200 MeV are expected to be generated for tumor therapy.During high-repetition operation,both prompt radiation and residual radiation may cause safety problems.Therefore,human radiological safety assessment before commissioning is essential.In this paper,we simulate both prompt and residual radiation using the Geant4 and FLUKA Monte Carlo codes with reasonable proton and as-produced electron beam parameters.We find that the prompt radiation can be shielded well by the concrete wall of the experimental hall,but the risk from residual radiation is nonnegligible and necessitates adequate radiation cooling.On the basis of the simulation results,we discuss the constraints imposed by radiation safety considerations on the annual working time,and we propose radiation cooling strategies for different shooting modes.展开更多
We further consider the effect of rod strength by employing the compressible penetration model to study the effect of compressibility on hypervelocity penetration.Meanwhile, we define different instances of penetratio...We further consider the effect of rod strength by employing the compressible penetration model to study the effect of compressibility on hypervelocity penetration.Meanwhile, we define different instances of penetration efficiency in various modified models and compare these penetration efficiencies to identify the effects of different factors in the compressible model. To systematically discuss the effect of compressibility in different metallic rod-target combinations, we construct three cases, i.e., the penetrations by the more compressible rod into the less compressible target, rod into the analogously compressible target, and the less compressible rod into the more compressible target. The effects of volumetric strain, internal energy, and strength on the penetration efficiency are analyzed simultaneously. It indicates that the compressibility of the rod and target increases the pressure at the rod/target interface. The more compressible rod/target has larger volumetric strain and higher internal energy. Both the larger volumetric strain and higher strength enhance the penetration or anti-penetration ability. On the other hand, the higher internal energy weakens the penetration or anti-penetration ability. The two trends conflict, but the volumetric strain dominates in the variation of the penetration efficiency, which would not approach the hydrodynamic limit if the rod and target are not analogously compressible. However, if the compressibility of the rod and target is analogous, it has little effect on the penetration efficiency.展开更多
Instability-induced wrinkle patterns of thin sheets are ubiquitous in nature,which often result in origami-like patterns that provide inspiration for the engineering of origami designs.Inspired by instability-induced ...Instability-induced wrinkle patterns of thin sheets are ubiquitous in nature,which often result in origami-like patterns that provide inspiration for the engineering of origami designs.Inspired by instability-induced origami patterns,we propose a computational origami design method based on the nonlinear analysis of loaded thin sheets and topology optimization.The bar-and-hinge model is employed for the nonlinear structural analysis,added with a displacement perturbation strategy to initiate out-of-plane buckling.Borrowing ideas from topology optimization,a continuous crease indicator is introduced as the design variable to indicate the state of a crease,which is penalized by power functions to establish the mapping relationships between the crease indicator and hinge properties.Minimizing the structural strain energy with a crease length constraint,we are able to evolve a thin sheet into an origami structure with an optimized crease pattern.Two examples with different initial setups are illustrated,demonstrating the effectiveness and feasibility of the method.展开更多
Study Objective: The purpose of the study is to present independent laboratory testing for a novel technology in air and on surfaces. Since 2020, public health goals have focused on improving indoor air quality. This ...Study Objective: The purpose of the study is to present independent laboratory testing for a novel technology in air and on surfaces. Since 2020, public health goals have focused on improving indoor air quality. This includes protection from airborne pathogens, such as tuberculosis, RSV, SARS-CoV-2, common cold or influenza viruses, measles, and others. Engineering controls are highly effective at reducing hazardous pathogens found in indoor air and from recontamination of surfaces. This occurs from a continuous cycle of settling of small, sustained airborne pathogens, which may become dehumidified, becoming airborne again, carried by room air currents around indoor spaces, then repeating the cycle. Methods: The novel technology utilizes a catalytic process to produce safe levels of hydrogen peroxide gas that are effective in reducing pathogens in the air and on surfaces. Air testing was performed with the MS2 bacteriophage, the test organism for ASHRAE standard 241, and methicillin-Resistant Staphylococcus aureus (MRSA). Surface testing was performed with SARS-COV-2 (Coronavirus COVID-19) and H1N1 (Influenza). Typical ventilation and filtration does not effectively remove disbursed pathogens from the entire facility, due to inconsistent air circulation and surface deposits of pathogens. Results: MS2 was reduced by 99.9%;MRSA was reduced by 99.9%;SARS-CoV-2 was reduced by 99.9%;H1N1 was reduced by 99.9%. Conclusion: This novel catalytic converter reduces a variety of pathogens in the air (99%) and on surfaces (99%), by actively disinfecting with the introduction of gaseous hydrogen peroxide. This active disinfection provides a strong solution for protecting the entire facility and its occupants.展开更多
X-ray sources with tunable energy spectra have a wide range of applications in different scenarios due to their different penetration depths.However,existing x-ray sources face difficulties in terms of energy regulati...X-ray sources with tunable energy spectra have a wide range of applications in different scenarios due to their different penetration depths.However,existing x-ray sources face difficulties in terms of energy regulation.In this paper,we present a scheme for tuning the energy spectrum of a betatron x-ray generated from a relativistic electron bunch oscillating in a plasma wakefield.The center energy of the x-ray source can be tuned from several keV to several hundred keV by changing the plasma density,thereby extending the control range by an order of magnitude.At different central energies,the brightness of the betatron radiation is in the range of 3.7×10^(22)to 5.5×10^(22)photons/(0.1%BW·s·mm^(2)·mrad^(2))and the photon divergence angle is about 2 mrad.This high-brightness,energy-controlled betatron source could pave the way to a wide range of applications requiring photons of specific energy,such as phase-contrast imaging in medicine,non-destructive testing and material analysis in industry,and imaging in nuclear physics.展开更多
Under the irradiation of ultraintense laser pulses,targets made of gas,solid,or artificial materials can generate high-energy electrons,ions,and X-rays comparable to conventional accelerators or national light source ...Under the irradiation of ultraintense laser pulses,targets made of gas,solid,or artificial materials can generate high-energy electrons,ions,and X-rays comparable to conventional accelerators or national light source facilities.Designing and creating high-performance targets are the core problems for laser acceleration.Nanotechnology and nanomaterials can help to build ideal targets that do not exist in nature.This paper reviews the advances in exploiting carbon nanotubes as outstanding targets for laser-driven particle acceleration in memory of Prof.Sishen Xie,the inventor of the fabrication method.We hope that the successful implementation of such targets in enhanced ion acceleration,high-efficiency electron acceleration,and brilliant X-ray generation could attract more interdiscipline interests and promote the development of this field.展开更多
We propose a new least squares finite element method to solve the Stokes problem with two sequential steps.The approximation spaces are constructed by the patch reconstruction with one unknown per element.For the firs...We propose a new least squares finite element method to solve the Stokes problem with two sequential steps.The approximation spaces are constructed by the patch reconstruction with one unknown per element.For the first step,we reconstruct an approximation space consisting of piecewise curl-free polynomials with zero trace.By this space,we minimize a least squares functional to obtain the numerical approximations to the gradient of the velocity and the pressure.In the second step,we minimize another least squares functional to give the solution to the velocity in the reconstructed piecewise divergence-free space.We derive error estimates for all unknowns under both L 2 norms and energy norms.Numerical results in two dimensions and three dimensions verify the convergence rates and demonstrate the great flexibility of our method.展开更多
The rigorous relations between the propagators of transient Schr¨odinger equations and stationary Green functions are established.Based on the generalized Fourier transform,non-singular transparent boundary condi...The rigorous relations between the propagators of transient Schr¨odinger equations and stationary Green functions are established.Based on the generalized Fourier transform,non-singular transparent boundary condition for transient problem is proposed in a representation of Green functions.The unified framework of Green function method is presented for converting an open boundary problem into a bounded boundary problem.Numerical scheme for time-dependent Schr¨odinger equation with non-singular transparent boundary condition is designed to simulate the propagations of a free Gaussian wave packet and the resonant tunnelling through double barriers.Numerical results validate the effectiveness of non-singular transparent boundary condition.展开更多
The discharged capillary plasma channel has been extensively studied as a high-gradient particle acceleration and transmission medium.A novel measurement method of plasma channel density profiles has been employed,whe...The discharged capillary plasma channel has been extensively studied as a high-gradient particle acceleration and transmission medium.A novel measurement method of plasma channel density profiles has been employed,where the role of plasma channels guiding the advantages of lasers has shown strong appeal.Here,we have studied the high-order transverse plasma density profile distribution using a channel-guided laser,and made detailed measurements of its evolution under various parameters.The paraxial wave equation in a plasma channel with high-order density profile components is analyzed,and the approximate propagation process based on the Gaussian profile laser is obtained on this basis,which agrees well with the simulation under phase conditions.In the experiments,by measuring the integrated transverse laser intensities at the outlet of the channels,the radial quartic density profiles of the plasma channels have been obtained.By precisely synchronizing the detection laser pulses and the plasma channels at various moments,the reconstructed density profile shows an evolution from the radial quartic profile to the quasi-parabolic profile,and the high-order component is indicated as an exponential decline tendency over time.Factors affecting the evolution rate were investigated by varying the incentive source and capillary parameters.It can be found that the discharge voltages and currents are positive factors quickening the evolution,while the electron-ion heating,capillary radii and pressures are negative ones.One plausible explanation is that quartic profile contributions may be linked to plasma heating.This work helps one to understand the mechanisms of the formation,the evolutions of the guiding channel electron-density profiles and their dependences on the external controllable parameters.It provides support and reflection for physical research on discharged capillary plasma and optimizing plasma channels in various applications.展开更多
We propose a numerical method to solve the Monge-Ampère equation which admits a classical convex solution.The Monge-Ampère equation is reformulated into an equivalent first-order system.We adopt a novel reco...We propose a numerical method to solve the Monge-Ampère equation which admits a classical convex solution.The Monge-Ampère equation is reformulated into an equivalent first-order system.We adopt a novel reconstructed discontinuous approximation space which consists of piecewise irrotational polynomials.This space allows us to solve the first-order system in two sequential steps.In the first step,we solve a nonlinear system to obtain the approximation to the gradient.A Newton iteration is adopted to handle the nonlinearity of the system.The approximation to the primitive variable is obtained from the approximate gradient by a trivial least squares finite element method in the second step.Numerical examples in both two and three dimensions are presented to show an optimal convergence rate in accuracy.It is interesting to observe that the approximation solution is piecewise convex.Particularly,with the reconstructed approximation space,the proposed method numerically demonstrates a remarkable robustness.The convergence of the Newton iteration does not rely on the initial values.The dependence of the convergence on the penalty parameter in the discretization is also negligible,in comparison to the classical discontinuous approximation space.展开更多
We propose an efficient scheme that combines density functional theory(DFT)with deep potentials(DPs),to systematically study convergence issues in the computation of the electronic thermal conductivity of warm dense a...We propose an efficient scheme that combines density functional theory(DFT)with deep potentials(DPs),to systematically study convergence issues in the computation of the electronic thermal conductivity of warm dense aluminum(2.7 g/cm^(3)and temperatures ranging from 0.5 eV to 5.0 eV)with respect to the number of k-points,the number of atoms,the broadening parameter,the exchange-correlation functionals,and the pseudopotentials.Furthermore,we obtain the ionic thermal conductivity using the Green–Kubo method in conjunction with DP molecular dynamics simulations,and we study size effects on the ionic thermal conductivity.This work demonstrates that the proposed method is efficient in evaluating both electronic and ionic thermal conductivities of materials.展开更多
The production of broadband,terawatt terahertz(THz) pulses has been demonstrated by irradiating relativistic lasers on solid targets.However,the generation of extremely powerful,narrow-band and frequency-tunable THz p...The production of broadband,terawatt terahertz(THz) pulses has been demonstrated by irradiating relativistic lasers on solid targets.However,the generation of extremely powerful,narrow-band and frequency-tunable THz pulses remains a challenge.Here,we present a novel approach for such THz pulses,in which a plasma wiggler is elaborated by a table-top laser and a near-critical density plasma.In such a wiggler,the laser-accelerated electrons emit THz radiations with a period closely related to the plasma thickness.The theoretical model and numerical simulations predict that a THz pulse with a laser-THz energy conversion of over 2.0%,an ultra-strong field exceeding 80 GV/m,a divergence angle of approximately 200 and a center frequency tunable from 4.4 to 1.5 THz can be generated from a laser of 430 mJ.Furthermore,we demonstrate that this method can work across a wide range of laser and plasma parameters,offering potential for future applications with extremely powerful THz pulses.展开更多
We study the stationary Wigner equation on a bounded, one- dimensional spatial domain with inflow boundary conditions by using the parity decomposition of L. Barletti and P. F. Zweifel [Transport Theory Statist. Phys....We study the stationary Wigner equation on a bounded, one- dimensional spatial domain with inflow boundary conditions by using the parity decomposition of L. Barletti and P. F. Zweifel [Transport Theory Statist. Phys., 2001, 30（4-6）： 507-520]. The decomposition reduces the half-range, two-point boundary value problem into two decoupled initial value problems of the even part and the odd part. Without using a cutoff approximation around zero velocity, we prove that the initial value problem for the even part is well-posed. For the odd part, we prove the uniqueness of the solution in the odd L2-spaee by analyzing the moment system. An example is provided to show that how to use the analysis to obtain the solution of the stationary Wigner equation with inflow boundary conditions.展开更多
文摘We perform a three-dimensional numerical simulation based on a one-step chemical reaction model to investigate changes in the mode of H2-Air detonation wave propagation from rotating detonation wave （RDW） mode to standing detonation wave mode. The physical characteristics of an RDW with injection velocity of 500 m/s are analyzed to investigate the physical mechanisms involved. We find that with increasing injection velocity, the detonation wave gradually changes from perpendicular to the head wall to parallel to the head wall. When the injection velocity exceeds the Chapman-Jouguet velocity VCJ （about 1984 m/s）, the detonation wave changes orientation to become perpendicular to the fuel injection direction, and the rotating mode changes accordingly to a standing mode. Finally, the plane detonation characteristic triple-wave structures can be found from the standing mode.
基金supported by the National Natural Science Foundation of China under Grant Nos.12122401 and 12074007.
文摘In traditional finite-temperature Kohn–Sham density functional theory(KSDFT),the partial occupation of a large number of high-energy KS eigenstates restricts the use of first-principles molecular dynamics methods at extremely high temperatures.However,stochastic density functional theory(SDFT)can overcome this limitation.Recently,SDFT and the related mixed stochastic–deterministic density functional theory,based on a plane-wave basis set,have been implemented in the first-principles electronic structure software ABACUS[Q.Liu and M.Chen,Phys.Rev.B 106,125132(2022)].In this study,we combine SDFT with the Born–Oppenheimer molecular dynamics method to investigate systems with temperatures ranging from a few tens of eV to 1000 eV.Importantly,we train machine-learning-based interatomic models using the SDFT data and employ these deep potential models to simulate large-scale systems with long trajectories.Subsequently,we compute and analyze the structural properties,dynamic properties,and transport coefficients of warm dense matter.
基金T W and D J S gratefully acknowledge the support of the Research Grants Council,Hong Kong SAR,through the Collaborative Research Fund Project No.8730054The work of H W is supported by the National Science Foundation of China under Grant Nos.11871110 and 12122103The work of W E is supported in part by a gift from iFlytek to Princeton University。
文摘To fill the gap between accurate(and expensive)ab initio calculations and efficient atomistic simulations based on empirical interatomic potentials,a new class of descriptions of atomic interactions has emerged and been widely applied;i.e.machine learning potentials(MLPs).One recently developed type of MLP is the deep potential(DP)method.In this review,we provide an introduction to DP methods in computational materials science.The theory underlying the DP method is presented along with a step-by-step introduction to their development and use.We also review materials applications of DPs in a wide range of materials systems.The DP Library provides a platform for the development of DPs and a database of extant DPs.We discuss the accuracy and efficiency of DPs compared with ab initio methods and empirical potentials.
基金supported by the National Natural Science Foundation of China (Nos.11921006,U2230133)Beijing Outstanding Young Scientists Program+2 种基金National Grand Instrument Project (No.2019YFF01014400)National Key R&D Program of China (No.2022YFA1603303)Open Foundation of Key Laboratory of High Power Laser and Physics,Chinese Academy of Sciences (No.SGKF202104)。
文摘In this study,we theoretically investigate the feasibility of using laser-wakefield accelerated(LWFA)electrons for the photonuclear measurement of nuclear isomers according to the characteristics of the electrons obtained from LWFA experiments conducted at the Compact Laser-Plasma Accelerator(CLAPA)laboratory.The experiments at the CLAPA show that a stable electron beam with an energy of 78-135 MeV and a charge of 300-600 pC can be obtained.The bremsstrahlung spectra were simulated using Geant4,which suggests that a bremsstrahlung source with a peak intensity of 10^(19) photons/s can be generated.Theoretical calculations of isomer production cross sections from the photonuclear reactions on six target nuclei,^(197)Au,^(180)Hf,^(159)Tb,^(115)In,^(103)Rh,and ^(90)Zr,were performed and compared with the available experimental data in EXFOR,which suggest that further experiments are required for a series of photonuclear reaction channels.Flux-averaged cross sections and isomer ratios(IR)resulting from such bremsstrahlung sources are theoretically deduced.The results suggest that IR measurements can be used to constrain nuclear components,such asγstrength function and optical model potential.In addition,the detection of the decay characteristics was evaluated with Geant4 simulations.The use of the LWFA electron beam and its bremsstrahlung for photonuclear studies involving nuclear isomers is anticipated.
基金This work was supported by the National Natural Science Foundation of China(Nos.12122501,11975037,61631001,and 11921006)the National Grand Instrument Project(Nos.2019YFF01014400,2019YFF01014404)the Foundation of Science and Technology on Plasma Physics Laboratory(No.6142A04220108).
文摘Ultrashort and powerful laser interactions with a target generate intense wideband electromagnetic pulses(EMPs).In this study,we report EMPs generated by the interactions between petawatt(30 fs,1.4×10^(20) W/cm^(2))femtosecond(fs)lasers with metal flat,plastic flat,and plastic nanowire-array(NWA)targets.Detailed analyses are conducted on the EMPs in terms of their spatial distribution,time and frequency domains,radiation energy,and protection.The results indicate that EMPs from metal targets exhibit larger amplitudes at varying angles than those generated by other types of targets and are enhanced significantly for NWA targets.Using a plastic target holder and increasing the laser focal spot can significantly decrease the radiation energy of the EMPs.Moreover,the composite shielding materials indicate an effective shielding effect against EMPs.The simulation results show that the NWA targets exert a collimating effect on thermal electrons,which directly affects the distribution of EMPs.This study provides guidance for regulating EMPs by controlling the laser focal spot,target parameters,and target rod material and is beneficial for electromagnetic-shielding design.
基金supported by National Grand Instrument Project (No.2019YFF01014404)the National Key Program for S & T Research and Development (No. 2022YFA1603202)+1 种基金National Natural Science Foundation of China (Nos.U2241281 and 11975037)the Foundation of Science and Technology on Plasma Physics Laboratory (No.6142A04220108)。
文摘Giant electromagnetic pulses(EMPs) induced by high-power laser irradiating solid targets interfere with various experimental diagnoses and even damage equipment,so unveiling the evolution of EMPs inside the laser chamber is crucial for designing effective EMP shielding.In this work,the transmission characteristics of EMPs as a function of distances from the target chamber center(TCC) are studied using B-dot probes.The mean EMP amplitude generated by picosecond laser-target interaction reaches 561 kV m^(-1),357 kV m^(-1),395 kV m^(-1),and 341 kV m^(-1)at 0.32 m,0.53 m,0.76 m,and 1 m from TCC,which decreases dramatically from 0.32 m to 0.53 m.However,it shows a fluctuation from 0.53 m to 1 m.The temporal features of EMPs indicate that time-domain EMP signals near the target chamber wall have a wider full width at half maximum compared to that close to TCC,mainly due to the echo oscillation of electromagnetic waves inside the target chamber based on simulation and experimentation.The conclusions of this study will provide a new approach to mitigate strong electromagnetic pulses by decreasing the echo oscillation of electromagnetic waves inside the target chamber during laser coupling with targets.
基金supported by the National Natural Science Foundation of China(Grant No.12205008)the NSFC Innovation Group Project(Grant No.11921006)+1 种基金the National Grand Instrument Project(Grant Nos.2019YFF01014402 and 2019YFF01014403)the National Science Fund for Distinguished Young Scholars(Grant No.12225501).
文摘The newly built Compact Laser Plasma Accelerator-Therapy facility at Peking University will deliver 60 J/1 Hz laser pulses with 30 fs duration.Driven by this petawatt laser facility,proton beams with energy up to 200 MeV are expected to be generated for tumor therapy.During high-repetition operation,both prompt radiation and residual radiation may cause safety problems.Therefore,human radiological safety assessment before commissioning is essential.In this paper,we simulate both prompt and residual radiation using the Geant4 and FLUKA Monte Carlo codes with reasonable proton and as-produced electron beam parameters.We find that the prompt radiation can be shielded well by the concrete wall of the experimental hall,but the risk from residual radiation is nonnegligible and necessitates adequate radiation cooling.On the basis of the simulation results,we discuss the constraints imposed by radiation safety considerations on the annual working time,and we propose radiation cooling strategies for different shooting modes.
文摘We further consider the effect of rod strength by employing the compressible penetration model to study the effect of compressibility on hypervelocity penetration.Meanwhile, we define different instances of penetration efficiency in various modified models and compare these penetration efficiencies to identify the effects of different factors in the compressible model. To systematically discuss the effect of compressibility in different metallic rod-target combinations, we construct three cases, i.e., the penetrations by the more compressible rod into the less compressible target, rod into the analogously compressible target, and the less compressible rod into the more compressible target. The effects of volumetric strain, internal energy, and strength on the penetration efficiency are analyzed simultaneously. It indicates that the compressibility of the rod and target increases the pressure at the rod/target interface. The more compressible rod/target has larger volumetric strain and higher internal energy. Both the larger volumetric strain and higher strength enhance the penetration or anti-penetration ability. On the other hand, the higher internal energy weakens the penetration or anti-penetration ability. The two trends conflict, but the volumetric strain dominates in the variation of the penetration efficiency, which would not approach the hydrodynamic limit if the rod and target are not analogously compressible. However, if the compressibility of the rod and target is analogous, it has little effect on the penetration efficiency.
基金National Key Research and Development Program of China(2020YFE0204200,2022YFB4701900)National Natural Science Foundation of China(11988102,12202008)Experiments for Space Exploration Program and the Qian Xuesen Laboratory,China Academy of Space Technology(TKTSPY-2020-03-05).
文摘Instability-induced wrinkle patterns of thin sheets are ubiquitous in nature,which often result in origami-like patterns that provide inspiration for the engineering of origami designs.Inspired by instability-induced origami patterns,we propose a computational origami design method based on the nonlinear analysis of loaded thin sheets and topology optimization.The bar-and-hinge model is employed for the nonlinear structural analysis,added with a displacement perturbation strategy to initiate out-of-plane buckling.Borrowing ideas from topology optimization,a continuous crease indicator is introduced as the design variable to indicate the state of a crease,which is penalized by power functions to establish the mapping relationships between the crease indicator and hinge properties.Minimizing the structural strain energy with a crease length constraint,we are able to evolve a thin sheet into an origami structure with an optimized crease pattern.Two examples with different initial setups are illustrated,demonstrating the effectiveness and feasibility of the method.
文摘Study Objective: The purpose of the study is to present independent laboratory testing for a novel technology in air and on surfaces. Since 2020, public health goals have focused on improving indoor air quality. This includes protection from airborne pathogens, such as tuberculosis, RSV, SARS-CoV-2, common cold or influenza viruses, measles, and others. Engineering controls are highly effective at reducing hazardous pathogens found in indoor air and from recontamination of surfaces. This occurs from a continuous cycle of settling of small, sustained airborne pathogens, which may become dehumidified, becoming airborne again, carried by room air currents around indoor spaces, then repeating the cycle. Methods: The novel technology utilizes a catalytic process to produce safe levels of hydrogen peroxide gas that are effective in reducing pathogens in the air and on surfaces. Air testing was performed with the MS2 bacteriophage, the test organism for ASHRAE standard 241, and methicillin-Resistant Staphylococcus aureus (MRSA). Surface testing was performed with SARS-COV-2 (Coronavirus COVID-19) and H1N1 (Influenza). Typical ventilation and filtration does not effectively remove disbursed pathogens from the entire facility, due to inconsistent air circulation and surface deposits of pathogens. Results: MS2 was reduced by 99.9%;MRSA was reduced by 99.9%;SARS-CoV-2 was reduced by 99.9%;H1N1 was reduced by 99.9%. Conclusion: This novel catalytic converter reduces a variety of pathogens in the air (99%) and on surfaces (99%), by actively disinfecting with the introduction of gaseous hydrogen peroxide. This active disinfection provides a strong solution for protecting the entire facility and its occupants.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.11921006 and 12175058)the Beijing Distinguished Young Scientist Program and National Grand Instrument Project (Grant No.SQ2019YFF01014400)+1 种基金the Beijing Municipal Science&Technology Commission,Administrative Commission of Zhongguancun Science Park (Grant No.Z231100006023003)in part funded by United Kingdom EPSRC (Grant Nos.EP/G054950/1,EP/G056803/1,EP/G055165/1,and EP/M022463/1)。
文摘X-ray sources with tunable energy spectra have a wide range of applications in different scenarios due to their different penetration depths.However,existing x-ray sources face difficulties in terms of energy regulation.In this paper,we present a scheme for tuning the energy spectrum of a betatron x-ray generated from a relativistic electron bunch oscillating in a plasma wakefield.The center energy of the x-ray source can be tuned from several keV to several hundred keV by changing the plasma density,thereby extending the control range by an order of magnitude.At different central energies,the brightness of the betatron radiation is in the range of 3.7×10^(22)to 5.5×10^(22)photons/(0.1%BW·s·mm^(2)·mrad^(2))and the photon divergence angle is about 2 mrad.This high-brightness,energy-controlled betatron source could pave the way to a wide range of applications requiring photons of specific energy,such as phase-contrast imaging in medicine,non-destructive testing and material analysis in industry,and imaging in nuclear physics.
基金supported by the following projects:the National Natural Science Foundation of China Innovation Group Project(No.11921006)National Grand Instrument Project(No.2019YFF01014402)National Science Fund for Distinguished Young Scholars(No.12225501).
文摘Under the irradiation of ultraintense laser pulses,targets made of gas,solid,or artificial materials can generate high-energy electrons,ions,and X-rays comparable to conventional accelerators or national light source facilities.Designing and creating high-performance targets are the core problems for laser acceleration.Nanotechnology and nanomaterials can help to build ideal targets that do not exist in nature.This paper reviews the advances in exploiting carbon nanotubes as outstanding targets for laser-driven particle acceleration in memory of Prof.Sishen Xie,the inventor of the fabrication method.We hope that the successful implementation of such targets in enhanced ion acceleration,high-efficiency electron acceleration,and brilliant X-ray generation could attract more interdiscipline interests and promote the development of this field.
基金supported by the Science Challenge Project(No.TZ2016002)the National Natural Science Foundation in China(No.11971041 and 11421101).
文摘We propose a new least squares finite element method to solve the Stokes problem with two sequential steps.The approximation spaces are constructed by the patch reconstruction with one unknown per element.For the first step,we reconstruct an approximation space consisting of piecewise curl-free polynomials with zero trace.By this space,we minimize a least squares functional to obtain the numerical approximations to the gradient of the velocity and the pressure.In the second step,we minimize another least squares functional to give the solution to the velocity in the reconstructed piecewise divergence-free space.We derive error estimates for all unknowns under both L 2 norms and energy norms.Numerical results in two dimensions and three dimensions verify the convergence rates and demonstrate the great flexibility of our method.
基金supported in part by the National Natural Science Foundation of China(Grant Nos.12171035,11671038).
文摘The rigorous relations between the propagators of transient Schr¨odinger equations and stationary Green functions are established.Based on the generalized Fourier transform,non-singular transparent boundary condition for transient problem is proposed in a representation of Green functions.The unified framework of Green function method is presented for converting an open boundary problem into a bounded boundary problem.Numerical scheme for time-dependent Schr¨odinger equation with non-singular transparent boundary condition is designed to simulate the propagations of a free Gaussian wave packet and the resonant tunnelling through double barriers.Numerical results validate the effectiveness of non-singular transparent boundary condition.
基金supported by the National Natural Science Foundation of China(Nos.11975037,61631001 and 11921006)the National Grand Instrument Project(Nos.2019YFF01014400 and 2019YFF01014404).
文摘The discharged capillary plasma channel has been extensively studied as a high-gradient particle acceleration and transmission medium.A novel measurement method of plasma channel density profiles has been employed,where the role of plasma channels guiding the advantages of lasers has shown strong appeal.Here,we have studied the high-order transverse plasma density profile distribution using a channel-guided laser,and made detailed measurements of its evolution under various parameters.The paraxial wave equation in a plasma channel with high-order density profile components is analyzed,and the approximate propagation process based on the Gaussian profile laser is obtained on this basis,which agrees well with the simulation under phase conditions.In the experiments,by measuring the integrated transverse laser intensities at the outlet of the channels,the radial quartic density profiles of the plasma channels have been obtained.By precisely synchronizing the detection laser pulses and the plasma channels at various moments,the reconstructed density profile shows an evolution from the radial quartic profile to the quasi-parabolic profile,and the high-order component is indicated as an exponential decline tendency over time.Factors affecting the evolution rate were investigated by varying the incentive source and capillary parameters.It can be found that the discharge voltages and currents are positive factors quickening the evolution,while the electron-ion heating,capillary radii and pressures are negative ones.One plausible explanation is that quartic profile contributions may be linked to plasma heating.This work helps one to understand the mechanisms of the formation,the evolutions of the guiding channel electron-density profiles and their dependences on the external controllable parameters.It provides support and reflection for physical research on discharged capillary plasma and optimizing plasma channels in various applications.
基金This research was supported by the National Natural Science Foundation in China(Nos.12201442,and 11971041).
文摘We propose a numerical method to solve the Monge-Ampère equation which admits a classical convex solution.The Monge-Ampère equation is reformulated into an equivalent first-order system.We adopt a novel reconstructed discontinuous approximation space which consists of piecewise irrotational polynomials.This space allows us to solve the first-order system in two sequential steps.In the first step,we solve a nonlinear system to obtain the approximation to the gradient.A Newton iteration is adopted to handle the nonlinearity of the system.The approximation to the primitive variable is obtained from the approximate gradient by a trivial least squares finite element method in the second step.Numerical examples in both two and three dimensions are presented to show an optimal convergence rate in accuracy.It is interesting to observe that the approximation solution is piecewise convex.Particularly,with the reconstructed approximation space,the proposed method numerically demonstrates a remarkable robustness.The convergence of the Newton iteration does not rely on the initial values.The dependence of the convergence on the penalty parameter in the discretization is also negligible,in comparison to the classical discontinuous approximation space.
基金supported by the Strategic Priority Research Program of the Chinese Academy of Sciences under Grant No.XDC01040100M.C.is supported by the National Science Foundation of China under Grant No.12074007.The numerical simulations were performed on the High Performance Computing Platform of CAPT.
文摘We propose an efficient scheme that combines density functional theory(DFT)with deep potentials(DPs),to systematically study convergence issues in the computation of the electronic thermal conductivity of warm dense aluminum(2.7 g/cm^(3)and temperatures ranging from 0.5 eV to 5.0 eV)with respect to the number of k-points,the number of atoms,the broadening parameter,the exchange-correlation functionals,and the pseudopotentials.Furthermore,we obtain the ionic thermal conductivity using the Green–Kubo method in conjunction with DP molecular dynamics simulations,and we study size effects on the ionic thermal conductivity.This work demonstrates that the proposed method is efficient in evaluating both electronic and ionic thermal conductivities of materials.
基金supported by the National Natural Science Foundation of China(Grant Nos.11921006 and 12175058)Beijing Distinguished Young Scientist Program and National Grand Instrument Project(Grant No.SQ2019YFF01014400)in part funded by UK EPSRC(Grant Nos.EP/G054950/1,EP/G056803/1,EP/G055165/1 and EP/M022463/1).
文摘The production of broadband,terawatt terahertz(THz) pulses has been demonstrated by irradiating relativistic lasers on solid targets.However,the generation of extremely powerful,narrow-band and frequency-tunable THz pulses remains a challenge.Here,we present a novel approach for such THz pulses,in which a plasma wiggler is elaborated by a table-top laser and a near-critical density plasma.In such a wiggler,the laser-accelerated electrons emit THz radiations with a period closely related to the plasma thickness.The theoretical model and numerical simulations predict that a THz pulse with a laser-THz energy conversion of over 2.0%,an ultra-strong field exceeding 80 GV/m,a divergence angle of approximately 200 and a center frequency tunable from 4.4 to 1.5 THz can be generated from a laser of 430 mJ.Furthermore,we demonstrate that this method can work across a wide range of laser and plasma parameters,offering potential for future applications with extremely powerful THz pulses.
文摘We study the stationary Wigner equation on a bounded, one- dimensional spatial domain with inflow boundary conditions by using the parity decomposition of L. Barletti and P. F. Zweifel [Transport Theory Statist. Phys., 2001, 30（4-6）： 507-520]. The decomposition reduces the half-range, two-point boundary value problem into two decoupled initial value problems of the even part and the odd part. Without using a cutoff approximation around zero velocity, we prove that the initial value problem for the even part is well-posed. For the odd part, we prove the uniqueness of the solution in the odd L2-spaee by analyzing the moment system. An example is provided to show that how to use the analysis to obtain the solution of the stationary Wigner equation with inflow boundary conditions.
基金financially supported by the National Natural Science Foundation of China(11921006)the National Key Research and Development Program of China(2020YFC1909604)+1 种基金Shenzhen Key Projects of Technological Research(JSGG20200925145800001)Shenzhen Science and Technology Program(CJGJZD20210408092801005)。
