2025-05-15 | | Total: 10
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.
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.
Physics-based models often involve large systems of parametrized partial differential equations, where design parameters control various properties. However, high-fidelity simulations of such systems on large domains or with high grid resolution can be computationally expensive for the accurate evaluation of a large number of parameters. Reduced-order modeling has emerged as a solution to reduce the dimensionality of such problems. This work focuses on a nonlinear compression technique using a convolutional autoencoder for accelerating the solution of transport in porous media problems. The model demonstrates successful training, achieving a mean square error (MSE) on the order of \num{1e-3} for the validation data. For an unseen parameter set, the model exhibits mixed performance; it achieves acceptable accuracy for larger time steps but shows lower performance for earlier times. This issue could potentially be resolved by fine-tuning the network architecture.
The group SU(3) has applications in several branches of physics. Many of these applications depend on availability of SU(3) coupling and recoupling coefficients. We have developed a modern Fortran library for calculation of the coupling coefficients, for both the SU(3)⊃U(1)×SU(2) and SU(3)⊃SO(3) group chains, and the recoupling coefficients. The library implements the algorithms of Draayer, Akiyama, and Millener, which are laid out in the paper. Performance of the library has been tested and compared to the Akiyama-Draayer (AD) library implementing the same algorithms as well as to a more recent implementation. Our library works for a larger range of SU(3) quantum numbers and provides more accurate coupling coefficients with large quantum numbers than the AD library.
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.
Surfactants have important effects on the dynamics of droplets on solid surfaces, which has inspired many industrial applications. Phase-field surfactant model with moving contact lines (PFS-MCL) has been employed to investigate the complex droplet dynamics with surfactants, while its numerical simulation remains challenging due to the coupling of gradient flows with respect to transport distances involving nonlinear and degenerate mobilities. We propose a novel structure-preserving variational scheme for PFS-MCL model with the dynamic boundary condition based on the minimizing movement scheme and optimal transport theory for Wasserstein gradient flows. The proposed scheme consists of a series of convex minimization problems and can be efficiently solved by our proposed primal-dual splitting method and its accelerated versions. By respecting the underlying PDE's variational structure with respect to the transport distance, the proposed scheme is proved to inherits the desirable properties including original energy dissipation, bound-preserving, and mass conservation. Through a suite of numerical simulations, we validate the performance of the proposed scheme and investigate the effects of surfactants on the droplet dynamics.
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.
SpectraMatcher is a cross-platform desktop application for interactive comparison of experimental and computed vibronic spectra, designed to assist in the recognition and assignment of spectral patterns. It provides an intuitive graphical interface -- with no coding or scripting required -- for importing experimental spectra, visualizing them alongside the corresponding theoretical spectra constructed from Gaussian frequency calculations, and adjusting key parameters such as peak width, intensity scaling factors, and vibration-type-specific anharmonic corrections. SpectraMatcher features an automated peak-matching algorithm that assigns experimental and computed peaks based on their intensity ratio and proximity. Assignments and spectra can be exported in multiple formats for publication or for further analysis. The software remains responsive even for large datasets, and supports efficient and reproducible interpretation of vibronic spectra.
Quantifying the uncertainty in predictive models is critical for establishing trust and enabling risk-informed decision making for personalized medicine. In contrast to one-size-fits-all approaches that seek to mitigate risk at the population level, digital twins enable personalized modeling thereby potentially improving individual patient outcomes. Realizing digital twins in biomedicine requires scalable and efficient methods to integrate patient data with mechanistic models of disease progression. This study develops an end-to-end data-to-decisions methodology that combines longitudinal non-invasive imaging data with mechanistic models to estimate and predict spatiotemporal tumor progression accounting for patient-specific anatomy. Through the solution of a statistical inverse problem, imaging data inform the spatially varying parameters of a reaction-diffusion model of tumor progression. An efficient parallel implementation of the forward model coupled with a scalable approximation of the Bayesian posterior distribution enables rigorous, but tractable, quantification of uncertainty due to the sparse, noisy measurements. The methodology is verified on a virtual patient with synthetic data to control for model inadequacy, noise level, and the frequency of data collection. The application to decision-making is illustrated by evaluating the importance of imaging frequency and formulating an optimal experimental design question. The clinical relevance is demonstrated through a model validation study on a cohort of patients with publicly available longitudinal imaging data.
We present a novel artificial intelligence approach to explore beyond Standard Model parameter spaces by leveraging a multi-objective optimisation algorithm. We apply this methodology to a non-minimal scotogenic model which is constrained by Higgs mass, anomalous magnetic moment of the muon, dark matter relic density, dark matter direct detection, neutrino masses and mixing, and lepton flavour violating processes. Our results successfully expand on the phenomenological realisations presented in previous work. We compare between multi- and single-objective algorithms and we observe more phenomenologically diverse solutions and an improved search capacity coming from the former. We use novelty detection to further explore sparsely populated regions of phenomenological interest. These results suggest a powerful search strategy that combines the global exploration of multi-objective optimisation with the exploitation of single-objective optimisation.