Related papers: Computationally Efficient Prediction of Area per L…
As one of the most robust global optimization methods, simulated annealing has received considerable attention, with many variations that attempt to improve the cooling schedule. This paper introduces a variant of simulated annealing that…
Approximate Bayes Computations (ABC) are used for parameter inference when the likelihood function of the model is expensive to evaluate but relatively cheap to sample from. In particle ABC, an ensemble of particles in the product space of…
A topology optimization formulation including a model of the layer-by-layer additive manufacturing (AM) process is considered. Defined as a multi-objective minimization problem, the formulation accounts for the performance and cost of both…
Continuum elastic models that account for membrane thickness variations are especially useful in the description of nanoscale deformations due to the presence of membrane proteins with hydrophobic mismatch. We show that terms involving the…
The equilibrium and fluctuation methods for determining the surface tension, $\sigma$, and bending modulus, $\kappa$, of a bilayer membrane with a fixed projected area are discussed. In the fluctuation method the elastic coefficients…
The process of numerical thermalization in particle-in-cell (PIC) simulations has been studied extensively. It is analogous to Coulomb collisions in real plasmas, causing particle velocity distributions (VDFs) to evolve towards a Maxwellian…
We have adapted a set of classification algorithms, also known as Machine Learning, to the identification of fluid and gel domains close to the main transition of dipalmitoyl-phosphatidylcholine (DPPC) bilayers. Using atomistic molecular…
A formal linear vector field representation for scientific equations is developed to rationalize the intuitive methods that are constantly employed. It is shown that unlike mechanical units that appear in the basis of the space, the reduced…
Machine Learning Interatomic Potentials (MLIPs) are becoming a central tool in simulation-based chemistry. However, like most deep learning models, MLIPs struggle to make accurate predictions on out-of-distribution data or when trained in a…
The use of explicit particle-in-cell (PIC) method for relativistic plasma simulations is restricted by numerical heating and instabilities that may significantly constrain the choice of time and space steps. To partially eliminate these…
Lipid membranes are abundant in living organisms, where they constitute a surrounding shell for cells and their organelles. There are many circumstances in which the deformations of lipid membranes are involved in living cells: fusion and…
Differentiable programming, enabled by automatic differentiation (AD), provides a robust framework for gradient-based optimization in computational plasma physics. While optimization is often only used towards design, we demonstrate that it…
Using molecular dynamics, this study investigates the local elastic properties of transverse shear deformation of lipid membranes. The analysis demonstrates that transverse shear deformation induces anisotropy in the local stress profile of…
Numerical simulations of turbulent Rayleigh-B\'enard convection in an ideal gas, using either the anelastic approximation or the fully compressible equations, are compared. Theoretically, the anelastic approximation is expected to hold in…
Simulation models often have parameters as input and return outputs to understand the behavior of complex systems. Calibration is the process of estimating the values of the parameters in a simulation model in light of observed data from…
Population annealing (PA) is a population-based algorithm that is designed for equilibrium simulations of thermodynamic systems with a rough free energy landscape. It is known to be more efficient in doing so than standard Markov chain…
Machine learning interatomic potentials (MLIPs) have become increasingly effective at approximating quantum mechanical calculations at a fraction of the computational cost. However, lower errors on held out test sets do not always translate…
High-entropy alloys (HEAs) have attracted increasing attention due to their unique structural and functional properties. In the study of HEAs, thermodynamic properties and phase stability play a crucial role, making phase diagram…
Nonlinear regression is a useful statistical tool, relating observed data and a nonlinear function of unknown parameters. When the parameter-dependent nonlinear function is computationally intensive, a straightforward regression analysis by…
We investigate the low-temperature phase of QCD and the crossover region with two light flavors of quarks. The chiral expansion around the point $(T,m=0)$ in the temperature vs. quark-mass plane indicates that a sharp real-time excitation…