Claudio Perego

Molecular-dynamics simulations of urea nucleation from aqueous solution

Significance Nucleation from solution is a ubiquitous process that plays important roles in physics, chemistry, engineering, and material science. Despite its importance, nucleation is far from being completely understood. In this work, we combine advanced molecular-dynamics simulation techniques and theory to provide a description of urea nucleation from aqueous solution. In particular, our analysis shows that a two-step nucleation mechanism is favorable and that two polymorphs are seen to compete in the early stages of the nucleation process. In our analysis, we have derived and validated a theoretical correction to finite-size effects to compute free-energy profiles in the limit of a macroscopic system at constant supersaturation. Despite its ubiquitous character and relevance in many branches of science and engineering, nucleation from solution remains elusive. In this framework, molecular simulations represent a powerful tool to provide insight into nucleation at the molecular scale. In this work, we combine theory and molecular simulations to describe urea nucleation from aqueous solution. Taking advantage of well-tempered metadynamics, we compute the free-energy change associated to the phase transition. We find that such a free-energy profile is characterized by significant finite-size effects that can, however, be accounted for. The description of the nucleation process emerging from our analysis differs from classical nucleation theory. Nucleation of crystal-like clusters is in fact preceded by large concentration fluctuations, indicating a predominant two-step process, whereby embryonic crystal nuclei emerge from dense, disordered urea clusters. Furthermore, in the early stages of nucleation, two different polymorphs are seen to compete.

Molecular dynamics simulations of solutions at constant chemical potential

Molecular dynamics studies of chemical processes in solution are of great value in a wide spectrum of applications, which range from nano-technology to pharmaceutical chemistry. However, these calculations are affected by severe finite-size effects, such as the solution being depleted as the chemical process proceeds, which influence the outcome of the simulations. To overcome these limitations, one must allow the system to exchange molecules with a macroscopic reservoir, thus sampling a grand-canonical ensemble. Despite the fact that different remedies have been proposed, this still represents a key challenge in molecular simulations. In the present work, we propose the Constant Chemical Potential Molecular Dynamics (CμMD) method, which introduces an external force that controls the environment of the chemical process of interest. This external force, drawing molecules from a finite reservoir, maintains the chemical potential constant in the region where the process takes place. We have applied the CμMD method to the paradigmatic case of urea crystallization in aqueous solution. As a result, we have been able to study crystal growth dynamics under constant supersaturation conditions and to extract growth rates and free-energy barriers.

Chemical potential calculations in dense liquids using metadynamics

The calculation of chemical potential has traditionally been a challenge in atomistic simulations. One of the most used approaches is Widoms insertion method in which the chemical potential is calculated by periodically attempting to insert an extra particle in the system. In dense systems this method fails since the insertion probability is very low. In this paper we show that in a homogeneous fluid the insertion probability can be increased using metadynamics. We test our method on a supercooled high density binary Lennard-Jones fluid. We find that we can obtain efficiently converged results even when Widoms method fails.

Advances in target normal sheath acceleration theory

A theoretical model of the Target Normal Sheath Acceleration (TNSA) process, able to go beyond the limits of available descriptions, is developed. It allows to achieve a more satisfactory interpretation of TNSA. The theory, also supported by two dimensional particle-in-cell simulations, elucidates the role played by the main laser and target parameters. Comparison between model predictions and experimental data related to the target thickness dependence of the maximum ion energy is discussed, showing satisfactory agreement. The model can be used as a simple but effective tool to guide the design of future experiments.

Target normal sheath acceleration analytical modeling, comparative study and developments

Ultra-intense laser interaction with solid targets appears to be an extremely promising technique to accelerate ions up to several MeV, producing beams that exhibit interesting properties for many foreseen applications. Nowadays, most of all the published experimental results can be theoretically explained in the framework of the target normal sheath acceleration (TNSA) mechanism proposed by Wilks et al. [Phys. Plasmas 8(2), 542 (2001)]. As an alternative to numerical simulation various analytical or semi-analytical TNSA models have been published in the latest years, each of them trying to provide predictions for some of the ion beam features, given the initial laser and target parameters. However, the problem of developing a reliable model for the TNSA process is still open, which is why the purpose of this work is to enlighten the present situation of TNSA modeling and experimental results, by means of a quantitative comparison between measurements and theoretical predictions of the maximum ion energy. Moreover, in the light of such an analysis, some indications for the future development of the model proposed by Passoni and Lontano [Phys. Plasmas 13(4), 042102 (2006)] are then presented.

Chemical potential calculations in non-homogeneous liquids

The numerical computation of chemical potential in dense non-homogeneous fluids is a key problem in the study of confined fluid thermodynamics. To this day, several methods have been proposed; however, there is still need for a robust technique, capable of obtaining accurate estimates at large average densities. A widely established technique is the Widom insertion method, which computes the chemical potential by sampling the energy of insertion of a test particle. Non-homogeneity is accounted for by assigning a density dependent weight to the insertion points. However, in dense systems, the poor sampling of the insertion energy is a source of inefficiency, hampering a reliable convergence. We have recently presented a new technique for the chemical potential calculation in homogeneous fluids. This novel method enhances the sampling of the insertion energy via well-tempered metadynamics, reaching accurate estimates at very large densities. In this paper, we extend the technique to the case of non-homogeneous fluids. The method is successfully tested on a confined Lennard-Jones fluid. In particular, we show that, thanks to the improved sampling, our technique does not suffer from a systematic error that affects the classic Widom method for non-homogeneous fluids, providing a precise and accurate result.

Preparation of the high power laser system PETAL for experimental studies of inertial confinement fusion and high energy density states of matter

The paper describes the preparation of the short-pulse high-energy laser PETAL that will be coupled to the French megajoule laser (LMJ) of CEA. The LMJ/PETAL facility will be opened to academic access for the international research community. In parallel diagnostics are being developed within the PETAL project and many physical problems are being addressed ranging from the study of the problems of radiation generation and activation issues to the problem of generation of large amplitude electromagnetic pulses.

SPECTRAL WIDTH VARIATION OF ULTRASHORT LASER PULSES IN MONOMODE OPTICAL FIBERS

Femtosecond laser pulse propagation in monomode optical fibers is demonstrated and investigated numerically (by simulations) and experimentally in this paper. A passively mode locked Nd:glass laser giving a pulse duration of about 200 fsec at 1053 nm wavelength and 120 mW average optical power with 100 MHz repetition rate is used in the experimental work. Numerical simulations are done by solving the nonlinear Schrodinger equation with the aid of Matlab program. The results show that self phase modulation (SPM) leads to compression of the spectral width from 5 nm to 2.1 nm after propagation of different optical powers (34, 43, 86 and 120 mW) in fibers of different length (5, 15, 35 m). The varying optical powers produced a varying phase shift. The output spectral width also changed with the fiber length at a given peak power.