In the framework of the Gibbs statistical theory, the question of the size of the particles forming the statistical system is investigated. This task is relevant for a wide variety of applications. The distribution for particle sizes and the moments of this quantity are determined from probabilistic considerations. The results are compared with different models and approximations. The particle size depends on the interactions in the system, on the compressibility factor, on the number of interacting particles. The cases of the absence and presence of a rigid incompressible particle core are considered.
We introduce a simple ansatz for the wavefunction of a many-body system based on coupled forward and backward-propagating semiclassical trajectories. This method is primarily aimed at, but not limited to, treating nonequilibrium dynamics in electron-phonon systems. The time-evolution of the system is obtained from the Euler-Lagrange variational principle, and we show that this ansatz yields Ehrenfest mean field theory in the limit that the forward and backward trajectories are orthogonal, and in the limit that they coalesce. We investigate accuracy and performance of this method by simulating electronic relaxation in the spin-boson model and the Holstein model. Although this method involves only pairs of semiclassical trajectories, it shows a substantial improvement over mean field theory, capturing quantum coherence of nuclear dynamics as well as electron-nuclear correlations. This improvement is particularly evident in nonadiabatic systems, where the accuracy of this coupled trajectory method extends well beyond the perturbative electron-phonon coupling regime. This approach thus provides an attractive route forward to the ab-initio description of relaxation processes, such as thermalization, in condensed phase systems.
Automatic Differentiation (AD) is a powerful tool that allows calculating derivatives of implemented algorithms with respect to all of their parameters up to machine precision, without the need to explicitly add any additional functions. Thus, AD has great potential in quantum chemistry, where gradients are omnipresent but also difficult to obtain, and researchers typically spend a considerable amount of time finding suitable analytical forms when implementing derivatives. Here, we demonstrate that automatic differentiation can be used to compute gradients with respect to any parameter throughout a complete quantum chemistry method. We implement DiffiQult, a fully autodifferentiable Hartree-Fock (HF) algorithm, which serves as a proof-of-concept that illustrates the capabilities of AD for quantum chemistry. We leverage the obtained gradients to optimize the parameters of one-particle basis sets in the context of the floating Gaussian framework.
Ordinary differential equation models have become a standard tool for the mechanistic description of biochemical processes. If parameters are inferred from experimental data, such mechanistic models can provide accurate predictions about the behavior of latent variables or the process under new experimental conditions. Complementarily, inference of model structure can be used to identify the most plausible model structure from a set of candidates, and thus gain novel biological insight. Several toolboxes can infer model parameters and structure for small- to medium-scale mechanistic models out of the box. However, models for highly multiplexed datasets can require hundreds to thousands of state variables and parameters. For the analysis of such large-scale models, most algorithms require intractably high computation times. This chapter provides an overview of state-of-the-art methods for parameter and model inference, with an emphasis on scalability.
The magnetic tuning of the low rotational levels in the Xsigma+(0,0,0), Api(0,0,0), and Bsigma+ (0,0,0) electronic states of strontium hydroxide, SrOH, have been experimentally investigated using high resolution optical Zeeman spectroscopy of a cold molecular beam sample. The observed Zeeman shifts and splittings are successfully modeled using a traditional effective Hamiltonian approach to account for the interaction between the Api and Bsigma+ states. The determined magnetic g-factors for the Xsigma+, Api , and Bsigma+ states are compared to those predicted by perturbation theory. The dispersed fluorescence resulting from laser excitation of rotationally resolved branch features of the Bsigma+ (0,0,0)<Xsigma+(0,0,0) Api(0,0,0)< Xsigma+(0,0,0) transitions have been recorded and analyzed. The measured fluorescence branching ratios are compared with Franck-Condon calculations. The required bending motion wave functions are derived using a discrete variable representation (DVR) method. Implications for laser slowing and magneto-optical trapping experiments for SrOH are described.