2025-05-15 | | Total: 78
Collisions between ultracold calcium monohydroxide (CaOH) molecules are realized and studied. Inelastic collision rate constants are measured for CaOH prepared in ground and excited vibrational states, and the electric field dependence of these rates is measured for molecules in single quantum states of the parity-doubled bending mode. Theoretical calculations of collision rate coefficients are performed and found to agree with measured values. The lowest collisional loss rates are for states with repulsive long-range potentials that shield ultracold molecules from loss channels at short distance. These results unveil the collisional behavior of parity doublet molecules in the ultracold regime, and lay the foundation for future experiments to evaporatively cool polyatomic molecules to quantum degeneracy.
An all-silica-fiber thulium-doped fiber laser emitting at 0.82 um on the transition from 3H4 to the ground state 3H6 outputs 105 W continuous-wave (CW) power and 555 W quasi-continuous-wave (QCW) instantaneous power with 0.96% duty cycle in 240-us rectangular pulses. The TDFL comprises a double-clad thulium-doped fiber (TDF) which is designed and fabricated in-house and is incorporated into an all-fiber cavity and cladding-pumped by five pigtailed diode lasers at 0.79 um. Co-lasing at 1.9 um counteracts population trapping in 3F4. The slope efficiency relative to absorbed pump power reaches 64% QCW and 77.5% CW. QCW, the beam quality M2 becomes 2.2 (beam parameter product BPP 0.57 mm mrad) and 2.45 (BPP 0.64) in orthogonal directions at ~250 W of instantaneous output power. Additionally, a modified QCW setup is continuously wavelength-tunable from 812 nm to 835 nm. We believe this is the first reported demonstration of high-power laser operation of the 3H4 to 3H6 transition in a TDF. Given also the simplicity and other attractions of an all-silica-fiber laser with direct-diode cladding-pumping, we believe our demonstration is valuable for applications ranging from laser machining of aluminum (benefitting from an absorption peak at 0.83 um) to scientific applications including strontium-based atomic clocks and cesium-based quantum metrology.
Relativistic calculations for the g factor of the lowest excited states 2p1/2 and 2p3/2 of lithium-like ions over a wide range of the nuclear charge numbers Z=10−92 are presented. Interelectronic interaction is considered within the perturbation theory up to the second order. One-loop QED contributions are calculated to all orders in αZ. Leading contributions of the nuclear recoil effects are taken into account. The quadratic and cubic terms in magnetic field are considered as well. A set of screening potentials is used in all calculations to estimate the unknown higher-order correlation effects.
Nonlinear pulse propagation in multimode fibers (MMFs) has attracted significant attention recently due to the rich spatiotemporal nonlinearities and promising applications. In practical scenarios, random birefringence in MMFs cannot be neglected, affecting the polarization-dependent nonlinear pulse propagation. This paper investigates the influence of random birefringence in MMFs on nonlinear ultrashort pulse propagation using a modified generalized multimode nonlinear Schrödinger equation. Two scenarios, spatial beam self-cleaning and multimode soliton propagation, are specifically examined. It is found that while random birefringence typically weakens nonlinearity in MMFs, certain nonlinear processes such as soliton self-frequency shift caused by intra-pulse Raman effect exhibit a complex relationship with random birefringence. Moreover, the study reveals that beam self-cleaning can endure random birefringence at high input peak powers. This research provides guidance for practical applications that utilize the nonlinear transmission of ultrashort pulses in MMFs.
Hexagonal boron nitride (hBN) has emerged as an excellent host material for generating room temperature single photons exhibiting high brightness and spin-photon entanglement. However, challenges in improving purity, stability, and scalability limit its use in quantum technologies. Here, we demonstrate highly pure and stable single photon emitters (SPEs) in hBN by directly growing carbon-doped, centimeter-scale hBN thin films using the pulsed laser deposition (PLD) method. These SPEs exhibit room temperature operation with polarized emission, achieving a g(2)(0) value of 0.015, which is among the lowest reported for room temperature SPEs and the lowest achieved for hBN SPEs. It also exhibits high brightness ( around 0.5 million counts per second), remarkable stability during continuous operation (> 15 minutes), and a Debye-Waller factor of 45%. First-principles calculations reveal unique carbon defects responsible for these properties, enabled by PLD's low-temperature synthesis and in-situ doping. Our results demonstrate an effective method for large-scale production of high-purity, stable SPEs in hBN, enabling robust quantum optical sources for various applications.
The single sideband (SSB) framework of analytical electron ptychography can account for the presence of residual geometrical aberrations induced by the probe-forming lens. However, the accuracy of this aberration correction method is highly sensitive to noise, in part due to the necessity of phase unwrapping. In this work, we thus propose two strategies to improve aberration correction performance in low-dose conditions: confining phase unwrapping within the sidebands and selecting only well-unwrapped sidebands for calculating aberration coefficients. These strategies are validated through SSB reconstructions of both simulated and experimental 4D-STEM datasets of monolayer tungsten diselenide (WSe2). A comparison of results demonstrates significant improvements in Poisson noise tolerance, making aberration correction more robust and reliable for low-dose imaging.
We present "torchGDM", a numerical framework for nano-optical simulations based on the Green's Dyadic Method (GDM). This toolkit combines a hybrid approach, allowing for both fully discretized nano-structures and structures approximated by sets of effective electric and magnetic dipoles. It supports simulations in three dimensions and for infinitely long, two-dimensional structures. This capability is particularly suited for multi-scale modeling, enabling accurate near-field calculations within or around a discretized structure embedded in a complex environment of scatterers represented by effective models. Importantly, torchGDM is entirely implemented in PyTorch, a well-optimized and GPU-enabled automatic differentiation framework. This allows for the efficient calculation of exact derivatives of any simulated observable with respect to various inputs, including positions, wavelengths or permittivity, but also intermediate parameters like Green's tensors. We anticipate that this toolkit will be valuable for applications merging nano-photonics and machine learning, as well as for solving nano-photonic optimization and inverse problems, such as the global design and characterization of metasurfaces, where optical interactions between structures are critical.
Motivated by the development of reusable software for computational chemistry [Lehtola, J. Chem. Phys. 159, 180901 (2023)], we aim to establish novel open source infrastructure for calculations with numerical atomic orbital (NAO) basis sets. Soft confinement potentials are typically used in this context to force the NAO radial basis functions {\psi}nl (r) to vanish smoothly in increasing r and to generate localized unoccupied states; we review such potentials in this work as a follow-up to our recent study on hard-wall confinement [Åström and Lehtola, J. Phys. Chem. A 129, 2791 (2025)]. We also note that in addition to their use in NAO generation, confinement potentials are also employed to simulate environmental effects in other research areas, such as studies of (i) atoms in solids, (ii) quantum dots, and (iii) high-pressure chemistry. As in our earlier work, we perform fully numerical density functional theory (DFT) calculations with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional, and employ spherically averaged, spin-restricted densities, as is usual in NAO studies. Our calculations employ the HelFEM program, which implements the finite element method (FEM), yielding variational energies and enabling the use of various boundary conditions. We consider four families of potentials to study the Mg and Ca atoms, which are textbook examples of extended electronic structures. We show that the resulting ground-state orbitals are surprisingly insensitive to the employed form of the confinement potential, and that the orbitals decay quickly under confinement. We study increasingly steep potentials and point out how they approach the hard-wall limit in a systematic and smooth manner. Finally, we assess NAO basis set truncation errors arising from various parameter choices for the singular potentials that are now broadly used in the NAO literature.
We adapt recent ideas for many-body chaos in nonlinear, Hamiltonian fluids [Murugan \textit{et al.}, Phys. Rev. Lett. 127, 124501 (2021)] to revisit the question of the Reynolds number Re dependence of the Lyapunov exponent λ∝Reα in fully developed turbulence. The use of such decorrelators allow us to investigate the interplay of the competing effects of viscous dissipation and nonlinearity. We obtain a precise value of α=0.59±0.04 and show that departure from the Kolmogorov mean field result λ∝√Re is a consequence of the intermittent fluctuations in the velocity-gradient tensor. The robustness of our results are further confirmed in a local, dynamical systems model for turbulence.
We carry out extensive pseudospectral direct numerical simulations (DNSs) of decaying three-dimensional (3D) Hall magnetohydrodynamics (3D HMHD) plasma turbulence at three magnetic Prandtl numbers Prm=0.1, 1.0 and 10.0. Our DNSs have been designed to uncover the dependence of the statistical properties of 3D HMHD turbulence on Prm and to bring out the subtle interplay between three lengths, the kinetic and magnetic dissipation length scales ηu, and ηb and the ion-inertial scale di, below which we see the manifestations of the Hall term. This interplay, qualitatively apparent from isosurface plots of the moduli of the vorticity and the current density, is exposed clearly by the kinetic-energy and magnetic-energy spectra, Eu(k) and Eb(k), respectively. We find two different inertial regions, In the first inertial region k<ki∼1/di, both the kinetic-energy and magnetic-energy spectra, Eu(k) and Eb(k), respectively, display power-law regions with an exponent that is consistent with Kolmogorov-type −5/3 scaling, at all values of Prm. In the second inertial region k>ki, the scaling of Eb(k) depends upon PrM: At Prm=0.1, the spectral-scaling exponent is −17/3, but for Prm=1 and 10 this exponent is −11/3. We then show theoretically that Eu(k)∼k2Eb(k) for Prm≪1 and Eb(k)∼k2Eu(k) for Prm≫1; our DNS results are consistent with our theoretical predictions. We examine, furthermore, left- and right-polarised fluctuations of the fields that lead, respectively, to the dominance of ion-cyclotron or whistler waves.
Classifying a crystalline solid's phase using X-ray diffraction (XRD) is a challenging endeavor, first because this is a poorly constrained problem as there are nearly limitless candidate phases to compare against a given experimental spectrum, and second because experimental signals are confounded by overlapping peaks, preferred orientations, and phase mixtures. To address this challenge, we develop Chem-XRD, a multimodal framework based on vision transformer (ViT) architecture with two defining features: (i) crystallographic data is constrained by elemental priors and (ii) phase is classified according to probabilities (not absolutes). Elemental information can be extracted from pre-synthesis precursors or post-synthesis elemental analysis. By combining structural and elemental modalities, Chem-XRD simultaneously predicts both the number and identity of phases of lead-halide perovskite materials and their decomposition products. Through integrated gradient calculations, we show that the model can adeptly adjust the contributions of structural and elemental modalities toward the final prediction of phase identification, achieving a level of interpretability beyond what self-attention mechanisms can achieve.
This paper examines "Free Creations of the Human Mind: The Worlds of Albert Einstein" by Diana Kormos Buchwald and Michael D. Gordin. The authors seek to dispel the long-standing myths of Einstein as the "lone genius" of Bern and the "stubborn sage" of Princeton, drawing on newly uncovered archival materials to illuminate his intellectual networks and collaborative engagements. By exploring the authors' reasoning, this paper engages with their interpretations, highlighting the strengths of their archival revelations and areas where alternative perspectives may enrich the understanding of Einstein's intellectual development.
Subskin modes are distinct from conventional skin modes as they localize not at the system's edge but rather below the edge. Unlike skin modes, where a substantial number of them can accumulate at the boundaries of a system due to the non-Hermitian skin effect, subskin modes are typically limited to one or a few and emerge only when the size of the system and the couplings are related in a very specific way. The nonlinear interaction can lift the restrictions on the couplings for the formation of these modes. The findings reveal the potential of subskin modes for practical applications and new research in non-Hermitian systems.
Silicon dominates electronics, optoelectronics, photovoltaics and photonics thanks to its suitable properties, abundance, and well-developed cost-effective manufacturing processes. Recently, crystalline silicon has been demonstrated to be an appealing alternative plasmonic material, both for the infrared where free-carrier plasmons are enabled by heavy doping, and for the ultraviolet where plasmonic effects are induced by interband transitions. Herein, we demonstrate that nanostructured amorphous silicon exhibits such so-called interband plasmonic properties in the ultraviolet, as opposed to the expectation that they would only arise in crystalline materials. We report optical plasmon resonances in the 100-to-300 nm wavelength range in ultrathin nanostructures. These resonances shift spectrally with the nanostructure shape and the nature of the surrounding matrix, while their field enhancement properties turn from surface plasmonic to epsilon-near-zero plasmonic. We present a vacuum ultraviolet wavelength- and polarization-selective ultrathin film absorber design based on deeply-subwavelength anisotropically-shaped nanostructures. These findings reveal amorphous silicon as a promising material platform for ultracompact and room-temperature-processed ultraviolet plasmonic devices operating down to vacuum ultraviolet wavelengths, for applications including anticounterfeiting, data encryption and storage, sensing and detection. Furthermore, these findings raise a fundamental question on how plasmonics can be based on amorphous nanostructures.
Meta-optics, the application of metasurfaces into optical systems, is seeing an accelerating development owing to advantages in size, weight and cost and the ability to program optical functions beyond traditional refractive optics. The transition of meta-optics from the laboratory into applications is enabled by scalable production methods based on highly reproducible semiconductor process technology. Here, we introduce a novel method for fabrication of double-sided metasurfaces through deep-UV lithography as a production-ready method for achieving high-quality meta-optics. We achieve patterning of a silicon wafer on both sides with mutual alignment of around 25 μm based on tool accuracy, without requiring through-wafer alignment markers other than the wafer notch. A first novel application highlighting the benefits of double-sided design is demonstrated in the form of a dual-band metalens with independent control over focal lengths in mid- and long-wave infrared bands. Using multi-reticle stitching we demonstrate a 40 mm diameter, large-area metalens with excellent broadband imaging performance, showing partial cancelling of chromatic dispersion when used in a hybrid configuration with a BaF2 refractive lens. Our work opens new avenues for infrared meta-optics designs and double-side meta-optics fabrication through a production-ready technique which can be directly translated into scalable technology for real-world applications.
We have considered a quantum droplet of two components of Bose-Einstein condensate (BEC) inside the electron of a Rydberg atom to study the surface mode of collective excitation using the Bogoliubov theory of excitation. We have calculated the surface excitation spectrum for various Rydberg electron-atom interaction strengths. From the energy spectrum, we calculated the surface tension of the droplet as a function of Rydberg electron-atom interaction strength. Our study shows that the electron-atom interaction enhances the surface energy; hence, the droplet will be more stable inside the electron of a Rydberg atom.
Metasurfaces (MSs) have recently become an efficient solution for controlling a spatial field distribution according to the demand of an electromagnetic device. In particular, anomalous reflectors based on impenetrable MSs can be realized by introducing a periodic spatial modulation of a local surface impedance defined on a flat or curvilinear surface. Two main methods to derive the profile of a purely reactive impedance that enable perfect anomalous reflection at an arbitrary angle are the power-flow conformal metasurface (PFCM) and auxiliary-field synthesis (AFS) methods. Both methods have been recently implemented with practical MS structures and have been successfully demonstrated in the experiments. However, the two methods assume drastically different synthesis algorithms and result in different impedance profiles. While the PFCM method is based on a direct analytical solution and does not require bounded wave excitation, the AFS method involves non-linear optimization of auxiliary bounded fields (surface waves) with rational selection of proper solutions. From a practical standpoint, the question of which method provides the most frequency-stabile solution becomes important to answer. Here we present a comprehensive analytical and numerical comparison of both methods applied to corrugated impedance surfaces in terms of the co-polarized anomalous reflection efficiency in a wide range of frequencies and corresponding scan angles. It is revealed that even though both methods can achieve imperceptible operational frequency bandwidths difference, the angular coverage of the frequency scan at larger reflection angles is up to 23deg wider for the PFCM method.
In this article, we studied quantum droplet with the Pöschl-Teller (PT) interaction potential between the Bose atoms. The Gross-Pitaevskii (GP) equation governs the system. The range and strength of the PT interaction can be adjusted. First, we studied the quantum droplet's density variation for various PT interaction parameters by the imaginary-time split-step Crank-Nicolson (CN) method. We then used the Bogoliubov theory to examine the collective excitation spectra. We observed that sharp roton forms and phonon modes are missing during long-range interactions. There is a gap at the zero momentum zone due to the long-range PT interaction, which increases with the range and strength of the interaction.
Ultrafast, time-resolved spectroscopies enable the direct observation of non-equilibrium processes in condensed-phase systems and have revealed key insights into energy transport, hydrogen-bond dynamics, and vibrational coupling. While ab initio molecular dynamics (AIMD) provides accurate, atomistic resolution of such dynamics, it becomes prohibitively expensive for non-equilibrium processes that require many independent trajectories to capture the stochastic nature of excitation and relaxation. To address this, we implemented a machine learning potential that incorporates time-dependent electric fields in a perturbative fashion, retaining AIMD-level accuracy. Using this approach, we simulate the time-dependent response of liquid water to a 12.3 THz Gaussian pump pulse (1.3 ps width), generating 32 ns of total trajectory data. With access to ab initio-quality electronic structure, we compute absorption coefficients and frequency-dependent refractive indices before, during, and after the pulse. The simulations reproduce key experimental observables, including transient birefringence and relaxation times. We observe energy transfer from the excited librational modes into translational and intramolecular vibrations, accompanied by a transient, nonlinear response of the hydrogen-bond network and a characteristic 0.7 ps timescale associated with librational energy dissipation. These findings demonstrate the method's ability to capture essential non-equilibrium dynamics with theoretical time-dependent IR spectroscopy and establish a broadly applicable framework for studying field-driven processes in complex molecular systems.
Large language models (LLMs) are now widely accessible, reaching learners at all educational levels. This development has raised concerns that their use may circumvent essential learning processes and compromise the integrity of established assessment formats. In physics education, where problem solving plays a central role in instruction and assessment, it is therefore essential to understand the physics-specific problem-solving capabilities of LLMs. Such understanding is key to informing responsible and pedagogically sound approaches to integrating LLMs into instruction and assessment. This study therefore compares the problem-solving performance of a general-purpose LLM (GPT-4o, using varying prompting techniques) and a reasoning-optimized model (o1-preview) with that of participants of the German Physics Olympiad, based on a set of well-defined Olympiad problems. In addition to evaluating the correctness of the generated solutions, the study analyzes characteristic strengths and limitations of LLM-generated solutions. The findings of this study indicate that both tested LLMs (GPT-4o and o1-preview) demonstrate advanced problem-solving capabilities on Olympiad-type physics problems, on average outperforming the human participants. Prompting techniques had little effect on GPT-4o's performance, while o1-preview almost consistently outperformed both GPT-4o and the human benchmark. Based on these findings, the study discusses implications for the design of summative and formative assessment in physics education, including how to uphold assessment integrity and support students in critically engaging with LLMs.
The software package DIALECT is introduced, which provides the capability of calculating excited-state properties and nonadiabatic dynamics of large molecular systems and can be applied to simulate energy and charge-transfer processes in molecular materials. To this end, we employ the FMO-LC-TDDFTB methodology, which combines the use of the fragment molecular orbital approach with the density-functional tight-binding method and an excitonic Hamiltonian including local and charge-transfer excitations. In this work, we present the features and capabilities of the DIALECT software package in simulating the excited state dynamics of molecules and molecular aggregates using exemplary trajectory surface hopping as well as decoherence corrected Ehrenfest dynamics calculations in the framework of LC-TDDFTB and FMO-LC-TDDFTB. In addition, the capability of simulating the polaritonic excited state properties is highlighted by the calculation of the polariton dispersion of an aggregate of naphthalene molecules. The development of the DIALECT program will facilitate the investigation of exciton and charge transport in large and complex molecular systems, such as biological aggregates, nanomaterials and other complex organic molecular systems.
To address the challenges of bidirectional tidal energy utilization efficiency and operational reliability of tidal turbines under low-flow conditions, this paper presents a novel hollow adaptive variable-pitch tidal energy generator based on symmetric airfoil design.In this document, the hydrodynamic performance of the device is analyzed by CFD method, and the energy capture and thrust load characteristics of the turbine are analyzed and discussed. The simulation results show that the power coefficient of the three-bladed turbine is better than that of the four-bladed and five-bladed turbine in the range of 2-4 tip speed ratios at a pitch angle of 10{\deg}, and the optimum power coefficient is 36.8%, which is higher than that of the ordinary axisymmetric wing turbine in terms of its energy acquisition efficiency. Its optimal power coefficient is 36.8%, and its energy efficiency is higher than that of ordinary axial symmetrical wing turbines.
Interplanetary coronal mass ejections (ICMEs) are major drivers of space weather disturbances, posing risks to both technological infrastructure and human activities. Automatic detection of ICMEs in solar wind in situ data is essential for early warning systems. While several methods have been proposed to identify these structures in time series data, robust real-time detection remains a significant challenge. In this work, we present ARCANE - the first framework explicitly designed for early ICME detection in streaming solar wind data under realistic operational constraints, enabling event identification without requiring observation of the full structure. Our approach evaluates the strengths and limitations of detection models by comparing a machine learning-based method to a threshold-based baseline. The ResUNet++ model, previously validated on science data, significantly outperforms the baseline, particularly in detecting high-impact events, while retaining solid performance on lower-impact cases. Notably, we find that using real-time solar wind (RTSW) data instead of high-resolution science data leads to only minimal performance degradation. Despite the challenges of operational settings, our detection pipeline achieves an F1 score of 0.53, with an average detection delay of 21.5% of the event's duration while only seeing a minimal amount of data. As more data becomes available, the performance increases significantly. These results mark a substantial step forward in automated space weather monitoring and lay the groundwork for enhanced real-time forecasting capabilities.
In the framework of the LHC Injectors Upgrade project, the CERN Proton Synchrotron underwent hardware changes during the Long Shutdown 2. These changes required an update of the CERN Proton Synchrotron impedance model and an assessment of their impact on beam quality. In the following, the updated CERN Proton Synchrotron impedance model is presented and then validated with beam-based measurements, comparing impedance-induced tune shifts and instability growth rates with simulations performed with a macroparticle tracking code, PyHEADTAIL, which includes the effects of beam-coupling impedance, betatron motion, non-linear synchrotron motion, and space charge. The interplay between first and second-order chromaticity, space-charge forces, and their effect on beam observables is thus studied in detail. Their respective impacts have been assessed and benchmarked with measurements to validate the PS transverse impedance model.
Over the past 7 decades, the classical Monte Carlo method has played a huge role in the fields of rarefied gas flow and micro/nano scale heat transfer, but it also has shortcomings: the time step and cell size are limited by the relaxation time and mean free path, making it difficult to efficiently simulate multi-scale heat and mass transfer problems from the ballistic to diffusion limit. To overcome this drawback, a unified gas-kinetic wave-particle (UGKWP) method is developed for solving the phonon Boltzmann transport equation (BTE) in all regimes covering both ballistic and diffusive limits. This method is built upon the space-time coupled evolution model of the phonon BTE, which provides the framework for constructing a multi-scale flux at the cell interfaces. At the same time, in order to capture non-equilibrium transport efficiently, the multi-scale flux comprises two distinct components: a deterministic part for capturing the near-equilibrium or diffusive transport and a statistical particle part for recovering non-equilibrium or ballistic transport phenomena. The UGKWP method exhibits remarkable multi-scale adaptability and versatility, seamlessly bridging the gap between the diffusive and ballistic transport phenomena. In the diffusive limit, the present method naturally converges to the Fourier's law, with the diminishing particle contribution, whereas in the ballistic limit, the non-equilibrium flux is fully described by the free-streaming particles. This inherent adaptability not only allows for precise capturing of both equilibrium and non-equilibrium heat transfer processes but also guarantees that the model adheres strictly to the underlying physical laws in each phonon transport regime.