Source author record

Michele Invernizzi

Michele Invernizzi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

5works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2022arXiv

Exploration vs Convergence Speed in Adaptive-bias Enhanced Sampling

In adaptive-bias enhanced sampling methods, a bias potential is added to the system to drive transitions between metastable states. The bias potential is a function of a few collective variables and is gradually modified according to the underlying free energy surface. We show that when the collective variables are suboptimal, there is an exploration-convergence tradeoff, and one must choose between a quickly converging bias that will lead to fewer transitions, or a slower to converge bias that can explore the phase space more efficiently but might require a much longer time to produce an accurate free energy estimate. The recently proposed On-the-fly Probability Enhanced Sampling (OPES) method focuses on fast convergence, but there are cases where fast exploration is preferred instead. For this reason, we introduce a new variant of the OPES method that focuses on quickly escaping metastable states, at the expense of convergence speed. We illustrate the benefits of this approach on prototypical systems and show that it outperforms the popular metadynamics method.

preprint2020arXiv

Attenuating the fermion sign problem in path integral Monte Carlo simulations using the Bogoliubov inequality and thermodynamic integration

Accurate thermodynamic simulations of correlated fermions using path integral Monte Carlo (PIMC) methods are of paramount importance for many applications such as the description of ultracold atoms, electrons in quantum dots, and warm-dense matter. The main obstacle is the fermion sign problem (FSP), which leads to an exponential increase in computation time both with increasing the system-size and with decreasing temperature. Very recently, Hirshberg et al. [J. Chem. Phys. 152, 171102 (2020)] have proposed to alleviate the FSP based on the Bogoliubov inequality. In the present work, we extend this approach by adding a parameter that controls the perturbation, allowing for an extrapolation to the exact result. In this way, we can also use thermodynamic integration to obtain an improved estimate of the fermionic energy. As a test system, we choose electrons in 2D and 3D quantum dots and find in some cases a speed-up exceeding 10^6 , as compared to standard PIMC, while retaining a relative accuracy of $\sim0.1\%$. Our approach is quite general and can readily be adapted to other simulation methods.

preprint2020arXiv

Path Integral Molecular Dynamics for Fermions: Alleviating the Sign Problem with the Bogoliubov Inequality

We present a method for performing path integral molecular dynamics (PIMD) simulations for fermions and address its sign problem. PIMD simulations are widely used for studying many-body quantum systems at thermal equilibrium. However, they assume that the particles are distinguishable and neglect bosonic and fermionic exchange effects. Interacting fermions play a key role in many chemical and physical systems, such as electrons in quantum dots and ultracold trapped atoms. A direct sampling of the fermionic partition function is impossible using PIMD since its integrand is not positive definite. We show that PIMD simulations for fermions are feasible by employing our recently developed method for bosonic PIMD and reweighting the results to obtain fermionic expectation values. The approach is tested against path integral Monte Carlo (PIMC) simulations for up to 7 electrons in a two-dimensional quantum dot for a range of interaction strengths. However, like PIMC, the method suffers from the sign problem at low temperatures. We propose a simple approach for alleviating it by simulating an auxiliary system with a larger average sign and obtaining an upper bound to the energy of the original system using the Bogoliubov inequality. This allows fermions to be studied at temperatures lower than would otherwise have been feasible using PIMD, as demonstrated in the case of a three-electron quantum dot. Our results extend the boundaries of PIMD simulations of fermions and will hopefully stimulate the development of new approaches for tackling the sign problem.

preprint2020arXiv

Rethinking Metadynamics: from bias potentials to probability distributions

Metadynamics is an enhanced sampling method of great popularity, based on the on-the-fly construction of a bias potential that is function of a selected number of collective variables. We propose here a change in perspective that shifts the focus from the bias to the probability distribution reconstruction, while keeping some of the key characteristics of metadynamics, such as the flexible on-the-fly adjustments to the free energy estimate. The result is an enhanced sampling method that presents a drastic improvement in convergence speed, especially when dealing with suboptimal and/or multidimensional sets of collective variables. The method is especially robust and easy to use, in fact it requires only few simple parameters to be set, and it has a straightforward reweighting scheme to recover the statistics of the unbiased ensemble. Furthermore it gives more control on the desired exploration of the phase space, since the deposited bias is not allowed to grow indefinitely and it does not push the simulation to uninteresting high free energy regions. We demonstrate the performance of the method in a number of representative examples.

preprint2015arXiv

Superlubric-Pinned Transition in Sliding Incommensurate Colloidal Monolayers

Two-dimensional (2D) crystalline colloidal monolayers sliding over a laser-induced optical lattice recently emerged as a new tool for the study of friction between ideal crystal surfaces. Here we focus in particular on static friction, the minimal sliding force necessary to depin one lattice from the other. If the colloid and the optical lattices are mutually commensurate, the colloid sliding is always pinned by static friction; but when they are incommensurate the presence or absence of pinning can be expected to depend upon the system parameters. If a 2D analogy to the mathematically established Aubry transition of one-dimensional systems were to hold, an increasing periodic corrugation strength $U_0$ should turn an initially free-sliding monolayer into a pinned state through a well-defined dynamical phase transition. We address this problem by the simulated sliding of a realistic model 2D colloidal lattice, confirming the existence of a clear and sharp superlubric-pinned transition for increasing corrugation strength. Unlike the 1D Aubry transition which is continuous, the 2D transition exhibits a definite first-order character. With no change of symmetry, the transition entails a structural character, with a sudden increase of the colloid-colloid interaction energy, accompanied by a compensating downward jump of the colloid-corrugation energy. The transition value for the corrugation amplitude $U_0$ depends upon the misalignment angle $θ$ between the optical and the colloidal lattices, superlubricity surviving until larger corrugations for angles away from the energetically favored orientation, which is itself generally slightly misaligned, as shown in recent work. The observability of the superlubric-pinned colloid transition is proposed and discussed.