Water molecule elimination from the protonated methanol dimer ion-An example of a size-selective intracluster reaction.

The abundance of extraterrestrial methanol makes the reaction between methanol molecules in a molecular cluster a possible key step in the search for mechanisms for the formation of more complex molecules under the conditions of the interstellar medium as well as circumstellar and planetary atmospheres. The reaction leading to the formation of the dimethyl ether ion from a methanol molecule interacting with a protonated methanol ion via the elimination of a water molecule is a basic mechanism for the formation of complex organic molecules. Here, we experimentally examine such reactions in the gas phase, analyzing the production and reactivity of protonated cluster ions formed by the ionization of a supersonic jet of methanol. Focusing especially on the post-collisional relaxation of the protonated methanol dimer and trimer ions after high-energy single collisions, the results indicate a strong size selectivity favoring the occurrence of this reaction only in the dimer ion. To elucidate this behavior, the velocity distribution of the eliminated water molecule was measured using an event-by-event coincidence analysis. These results are interpreted using quantum chemical calculations of the dissociation pathways. It turns out that in the dimer case, two transition states are able to contribute to this intracluster reaction. In the trimer case, methanol evaporation appears as the most energetically favorable relaxation pathway.

Coherent state switching using vibrational polaritons in an asymmetric double-well potential.

The quantum dynamics of vibrational polaritonic states arising from the interaction of a bistable molecule with the quantized mode of a Fabry-Perot microcavity is investigated using a generic asymmetric double-well potential as a simplified one-dimensional model of a reactive molecule. After discussing the role of the light-matter coupling strength in the emergence of avoided crossings between polaritonic states, we investigate the possibility of using these crossings to trigger a dynamical switching of these states from one potential well to the other. Two schemes are proposed to achieve this coherent state switching, either by preparing the molecule in an appropriate vibrational excited state before inserting it into the cavity, or by applying a short laser pulse inside the cavity to obtain a coherent superposition of polaritonic states. The respective influences of dipole moment amplitude and potential asymmetry on the coherent switching process are also discussed.

Vibrational Circular Dichroism Spectroscopy with a Classical Polarizable Force Field: Alanine in the Gas and Condensed Phases.

Polarizable force fields are an essential component for the chemically accurate modeling of complex molecular systems with a significant degree of fluxionality, beyond harmonic or perturbative approximations. In this contribution we examine the performance of such an approach for the vibrational spectroscopy of the alanine amino acid, in the gas and condensed phases, from the Fourier transform of appropriate time correlation functions generated along molecular dynamics (MD) trajectories. While the infrared (IR) spectrum only requires the electric dipole moment, the vibrational circular dichroism (VCD) spectrum further requires knowledge of the magnetic dipole moment, for which we provide relevant expressions to be used with polarizable force fields. The AMOEBA force field was employed here to model alanine in the neutral and zwitterionic isolated forms, solvated by water or nitrogen, and as a crystal. Within this framework, comparison of the electric and magnetic dipole moments to those obtained with nuclear velocity perturbation theory based on density-functional theory for the same MD trajectories are found to agree well with one another. The statistical convergence of the IR and VCD spectra is examined and found to be more demanding in the latter case. Comparisons with experimental frequencies are also provided for the condensed phases.

Modeling the dynamics of quantum systems coupled to large-dimensional baths using effective energy states.

The quantum dynamics of a low-dimensional system in contact with a large but finite harmonic bath is theoretically investigated by coarse-graining the bath into a reduced set of effective energy states. In this model, the couplings between the system and the bath are obtained from statistically averaging over the discrete, degenerate effective states. Our model is aimed at intermediate bath sizes in which non-Markovian processes and energy transfer between the bath and the main system are important. The method is applied to a model system of a Morse oscillator coupled to 40 harmonic modes. The results are found to be in excellent agreement with the direct quantum dynamics simulations presented in the work of Bouakline et al. [J. Phys. Chem. A 116, 11118-11127 (2012)], but at a much lower computational cost. Extension to larger baths is discussed in comparison to the time-convolutionless method. We also extend this study to the case of a microcanonical bath with finite initial internal energies. The computational efficiency and convergence properties of the effective bath states model with respect to relevant parameters are also discussed.

The interplay between size, shape, and surface segregation in high-entropy nanoalloys.

The respective influences of particle shape and size on the energetic stability of five-component multimetallic nanoparticles have been computationally investigated for AlCuFeCrNi and AuCuPdNiCo mixtures at equiconcentration. Using available embedded-atom model potentials, exchange Monte Carlo simulations possibly assisted with systematic quenching, we explore tools to approach ideal phase equilibrium in such high-entropy nanoalloys. In particular, we show how deviations to ideal solid solution behaviors can be characterized using percolation analyses, and how the contribution of alloying fluctuations at finite temperature can be inferred to evaluate the entropy of mixing in such nonideal cases. An approximation to the entropy of mixing based on pair correlations only is also found to capture the behavior of the thermodynamical mixing entropy quite well, and can be used as an order parameter of mixing. While the AlCuFeCrNi mixture appears to mix reasonably well in all cases considered, cobalt and nickel segregate significantly in AuCuPdNiCo nanoparticles, deviating strongly from ideal random mixtures. A simple Gaussian regression model applied to a coarse distribution of concentrations is found to correctly predict conditions for optimising the mixing thermodynamical properties of the miscible AlCuFeCrNi nanoparticle.

Influence of the environment on the infrared spectrum of alanine: An effective mode analysis.

The vibrational spectrum of the alanine amino acid was computationally determined in the infrared range 1000-2000 cm-1, under various environments encompassing the gas, hydrated, and crystalline phases, by means of classical molecular dynamics trajectories, carried out with the Atomic Multipole Optimized Energetics for Biomolecular Simulation polarizable force field. An effective mode analysis was performed, in which the spectra are optimally decomposed into different absorption bands arising from well-defined internal modes. In the gas phase, this analysis allows us to unravel the significant differences between the spectra obtained for the neutral and zwitterionic forms of alanine. In condensed phases, the method provides invaluable insight into the molecular origins of the vibrational bands and further shows that peaks with similar positions can be traced to rather different molecular motions.

Glycine Peptide Chain Formation in the Gas Phase via Unimolecular Reactions.

Peptide chain formation from amino acids such as glycine is a key step in the emergence of life. Unlike their synthesis by living systems, how peptide chains grow under abiotic conditions is an open question given the variety of organic compounds discovered in various astrophysical environments, comets and meteorites. We propose a new abiotic route in the presence of protonated molecular dimers of glycine in a cold gaseous atmosphere without further need for a solid catalytic substrate. The results provide evidence for the preferential formation of mixed protonated dimers of glycine consisting of a dipeptide and a glycine molecule instead of pure protonated glycine dimers. Additional measurements mimicking a cosmic-ray impact in terms of internal excitation show that a single gas-phase collision induces polymerization via dehydration in both the mixed and pure dimer ions. Peptide chain growth is thus demonstrated to occur via a unimolecular gas-phase reaction in an excited cluster ion.

Concluding remarks.

In 2007 in Birmingham, Riccardo Ferrando, Sarah Horswell, Claude Henry, Brian Johnson, Peter Lievens and the late Roy Johnston, organised a first Faraday Discussions on the topic of nanoparticles made of multiple metals, setting the stage for what would then become known as "nanoalloys". A second edition of these discussions, held in London in 2022, shows how significantly the field of nanoalloys has evolved in fifteen years.

Finite-temperature stability of hydrocarbons: Fullerenes vs flakes.

The effects of a finite temperature on the equilibrium structures of hydrocarbon molecules are computationally explored as a function of size and relative chemical composition in hydrogen and carbon. Using parallel tempering Monte Carlo simulations employing a reactive force field, we find that in addition to the phases already known for pure carbon, namely, cages, flakes, rings, and branched structures, strong changes due to temperature and the addition of little amounts of hydrogen are reported. Both entropy and the addition of moderate amounts of hydrogen favor planar structures such as nanoribbons over fullerenes. Accurate phase diagrams are proposed, highlighting the possible presence of multiple phase changes at finite size and composition. Astrophysical implications are also discussed.

Solvation of Large Polycyclic Aromatic Hydrocarbons in Helium: Cationic and Anionic Hexabenzocoronene.

The adsorption of helium on charged hexabenzocoronene (Hbc, C42H18), a planar polycyclic aromatic hydrocarbon (PAH) molecule of D6h symmetry, was investigated by a combination of high-resolution mass spectrometry and classical and quantum computational methods. The ion abundance of HenHbc+ complexes versus size n features prominent local anomalies at n = 14, 38, 68, 82, and a weak one at 26, indicating that for these "magic" sizes, the helium evaporation energies are relatively large. Surprisingly, the mass spectra of anionic HenHbc- complexes feature a different set of anomalies, namely at n = 14, 26, 60, and 62, suggesting that the preferred arrangement of the adsorbate atoms depends on the charge of the substrate. The results of our quantum calculations show that the adsorbate layer grows by successive filling of concentric rings that surround the central benzene ring, which is occupied by one helium atom each on either side of the substrate. The helium atoms are fairly localized in filled rings and they approximately preserve the D6h symmetry of the substrate, but helium atoms in partially filled rings are rather delocalized. The first three rings contain six atoms each; they account for magic numbers at n = 14, 26, and 38. The size of the first ring shrinks as atoms are filled into the second ring, and the position of atoms in the second ring changes from hollow sites to bridge sites as atoms are filled into the third ring. Beyond n = 38, however, the arrangement of helium atoms in the first three rings remains essentially frozen. Presumably, another ring is filled at n = 68 for cations and n = 62 for anions. The calculated structures and energies do not account for the difference between charge states, although they agree with the measurements for the cations and show that the first solvation shell of Hbc± is complete at n = 68. Beyond that size, the adsorbate layer becomes three-dimensional, and the circular arrangement of helium changes to hexagonal.

Adsorption of Helium and Hydrogen on Triphenylene and 1,3,5-Triphenylbenzene.

The adsorption of helium or hydrogen on cationic triphenylene (TPL, C18H12), a planar polycyclic aromatic hydrocarbon (PAH) molecule, and of helium on cationic 1,3,5-triphenylbenzene (TPB, C24H18), a propeller-shaped PAH, is studied by a combination of high-resolution mass spectrometry and classical and quantum computational methods. Mass spectra indicate that HenTPL+ complexes are particularly stable if n = 2 or 6, in good agreement with the quantum calculations that show that for these sizes, the helium atoms are strongly localized on either side of the central carbon ring for n = 2 and on either side of the three outer rings for n = 6. Theory suggests that He14TPL+ is also particularly stable, with the helium atoms strongly localized on either side of the central and outer rings plus the vacancies between the outer rings. For HenTPB+, the mass spectra hint at enhanced stability for n = 2, 4 and, possibly, 11. Here, the agreement with theory is less satisfactory, probably because TPB+ is a highly fluxional molecule. In the global energy minimum, the phenyl groups are rotated in the same direction, but when the zero-point harmonic correction is included, a structure with one phenyl group being rotated opposite to the other two becomes lower in energy. The energy barrier between the two isomers is very small, and TPB+ could be in a mixture of symmetric and antisymmetric states, or possibly even vibrationally delocalized.

Clustering, collision, and relaxation dynamics in pure and doped helium nanoclusters: Density- vs particle-based approaches.

The clustering, collision, and relaxation dynamics of pristine and doped helium nanodroplets is theoretically investigated in cases of pickup and clustering of heliophilic argon, collision of heliophobic cesium atoms, and coalescence of two droplets brought into contact by their mutual long-range van der Waals interaction. Three approaches are used and compared with each other. The He time-dependent density functional theory method considers the droplet as a continuous medium and accounts for its superfluid character. The ring-polymer molecular dynamics method uses a path-integral description of nuclear motion and incorporates zero-point delocalization while bosonic exchange effects are ignored. Finally, the zero-point averaged dynamics approach is a mixed quantum-classical method in which quantum delocalization is described by attaching a frozen wavefunction to each He atom, equivalent to classical dynamics with effective interaction potentials. All three methods predict that the growth of argon clusters is significantly hindered by the helium host droplet due to the impeding shell structure around the dopants and kinematic effects freezing the growing cluster in metastable configurations. The effects of superfluidity are qualitatively manifested by different collision dynamics of the heliophilic atom at high velocities, as well as quadrupole oscillations that are not seen with particle-based methods, for droplets experiencing a collision with cesium atoms or merging with each other.

Dimerization dynamics of carboxylic acids in helium nanodroplets.

The dimerization of molecules in helium nanodroplets is known to preferentially yield structures of higher energy than the global energy minimum structure for a number of quite different monomers. Here, we explore dimerization in this environment using an atomistic model within statistically converged molecular dynamics (MD) trajectories, treating the solvent implicitly through the use of a thermostat, or more explicitly by embedding one monomer in a He100 cluster. The focus is on the two simplest carboxylic acids, formic and acetic, both of which have been studied experimentally. While the global minimum structure, which comprises two CO⋯HO hydrogen bonds, is predicted to be the most abundant dimer in the absence of the helium solvent, this is no longer the case once helium atoms are included. The simulations confirm the importance of kinetic trapping effects and also shed light on the occurrence of specific dynamical effects, leading to the occasional formation of high-energy structures away from minima, such as saddle configurations. Theoretically predicted infrared spectra, based on the MD statistics, are in good agreement with the experimental spectra.

Energy Dispersion in Pyridinium-Water Nanodroplets upon Irradiation.

Postirradiation dissociation of molecular clusters has been mainly studied assuming energy redistribution in the entire cluster prior to the dissociation. Here, the evaporation of water molecules from out-of-equilibrium pyridinium-water cluster ions was investigated using the recently developed correlated ion and neutral time-of-flight (COINTOF) mass spectrometry technique in combination with a velocity-map imaging (VMI) device. This special setup enables the measurement of velocity distributions of the evaporated molecules upon high-velocity collisions with an argon atom. The distributions measured for pyridinium-water cluster ions are found to have two distinct components. Besides a low-velocity contribution, which corresponds to the statistical evaporation of water molecules after nearly complete redistribution of the excitation energy within the clusters, a high-velocity contribution is also found in which the molecules are evaporated before the energy redistribution is complete. These two different evaporation modes were previously observed and described for protonated water cluster ions. However, unlike in the case of pure water clusters, the low-velocity part of the distributions for pyridinium-doped water clusters is itself composed of two distinct Maxwell-Boltzmann distributions, indicating that evaporated molecules originate in this case from out-of-equilibrium processes. Statistical molecular dynamics simulations were performed to (i) understand the effects caused in the ensuing evaporation process by the various excitation modes at different initial cluster constituents and to (ii) simulate the distributions resulting from sequential evaporations. The presence of a hydrophobic impurity in water clusters is shown to impact water molecule evaporation due to the energy storage in the internal degrees of freedom of the impurity.

Adsorption of Helium on Small Cationic PAHs: Influence of Hydrocarbon Structure on the Microsolvation Pattern.

The adsorption of up to ∼100 helium atoms on cations of the planar polycyclic aromatic hydrocarbons (PAHs) anthracene, phenanthrene, fluoranthene, and pyrene was studied by combining helium nanodroplet mass spectrometry with classical and quantum computational methods. Recorded time-of-flight mass spectra reveal a unique set of structural features in the ion abundance as a function of the number of attached helium atoms for each of the investigated PAHs. Path-integral molecular dynamics simulations were used with a polarizable potential to determine the underlying adsorption patterns of helium around the studied PAH cations and in good general agreement with the experimental data. The calculated structures of the helium-PAH complexes indicate that the arrangement of adsorbed helium atoms is highly sensitive toward the structure of the solvated PAH cation. Closures of the first solvation shell around the studied PAH cations are suggested to lie between 29 and 37 adsorbed helium atoms depending on the specific PAH cation. Helium atoms are found to preferentially adsorb on these PAHs following the 3×3 commensurate pattern common for graphitic surfaces, in contrast to larger carbonaceous molecules like corannulene, coronene, and fullerenes that exhibit a 1 × 1 commensurate phase.

Infrared Spectroscopy of Chemically Diverse Carbon Clusters: A Data-Driven Approach.

Carbon clusters exhibit a broad diversity of topologies and shapes, encompassing fullerene-like cages, graphene-like flakes, and more disordered pretzel-like and branched structures. Here, we examine computationally their infrared spectra in relation with these structures from a statistical perspective. Individual spectra for broad samples of isomers were determined by means of the self-consistent charge density functional-based tight-binding method, and an interpolation scheme is designed to reproduce the spectral features by regression on a much smaller subset of the sample. This interpolation proceeds by encoding the structures using appropriate descriptors and selecting them through principal component analysis, Gaussian regression or inverse distance weighting providing the nonlinear weighting functions. Metric learning is employed to reduce the global error on a preselected testing set. The interpolated spectra satisfactorily reproduce the specific spectral features and their dependence on the size and shape, enabling quantitative prediction away from the testing set. Finally, the classification of structures within the four proposed families is critically discussed through a statistical analysis of the sample based on iterative label spreading.

Assessing cluster models of solvation for the description of vibrational circular dichroism spectra: synergy between static and dynamic approaches.

Solvation effects are essential for defining the shape of vibrational circular dichroism (VCD) spectra. Several approaches have been proposed to include them into computational models for calculating VCD signals, in particular those resting on the "cluster-in-a-liquid" model. Here we examine the capabilities of this ansatz on the example of flexible (1S,2S)-trans-1-amino-2-indanol solvated in dimethyl sulfoxide (DMSO). We compare cluster sets obtained from static calculations with results from explicit molecular dynamics (MD) trajectories based on either force field (FF) or first-principles (FP) methods. While the FFMD approach provides a broader sampling of configurational space, FPMD and time-correlation functions of dipole moments account for anharmonicity and entropy effects in the VCD calculation. They provide a means to evaluate the immediate effect of the solvent on the spectrum. This survey singles out several challenges associated with the use of clusters to describe solvation effects in systems showing shallow potential energy surfaces and non-covalent interactions. Static structures of clusters involving a limited number of solvent molecules satisfactorily capture the main effects of solvation in the bulk limit on the VCD spectra, if these structures are correctly weighted. The importance of taking into consideration their fluxionality, i.e. different solvent conformations sharing a same hydrogen bond pattern, and the limitations of small clusters for describing the solvent dynamics are discussed.

How interface properties control the equilibrium shape of core-shell Fe-Au and Fe-Ag nanoparticles.

While combining two metals in the same nanoparticle can lead to remarkable novel applications, the resulting structure in terms of crystallinity and shape remains difficult to predict. It is thus essential to provide a detailed atomistic picture of the underlying growth processes. In the present work we address the case of core-shell Fe-Au and Fe-Ag nanoparticles. Interface properties between Fe and the noble metals Au and Ag, computed using DFT, were used to parameterize Fe-Au and Fe-Ag pairwise interactions in combination with available many-body potentials for the pure elements. The growth of Au or Ag shells on nanometric Fe cores with prescribed shapes was then modelled by means of Monte Carlo simulations. The shape of the obtained Fe-Au nanoparticles is found to strongly evolve with the amount of metal deposited on the Fe core, a transition from the polyhedral Wulff shape of bare iron to a cubic shape taking place as the amount of deposited gold exceeds two monolayers. In striking contrast, the growth of silver proceeds in a much more anisotropic, Janus-like way and with a lesser dependence on the iron core shape. In both cases, the predicted morphologies are found to be in good agreement with experimental observations in which the nanoparticles are grown by physical deposition methods. Understanding the origin of these differences, which can be traced back to subtle variations in the electronic structure of the Au/Fe and Ag/Fe interfaces, should further contribute to the better design of core-shell bimetallic nanoparticles.

The structures of cationic gold clusters probed by far-infrared spectroscopy.

Determining the precise structures of small gold clusters is an essential step towards understanding their chemical and physical properties. Due to the relativistic nature of gold, its clusters remain planar (2D) up to appreciable sizes. Ion mobility experiments have suggested that positively charged gold clusters adopt three-dimensional (3D) structures from n = 8 onward. Computations predict, depending on the level of theory, 2D or 3D structures as putative energy-minimum for n = 8. In this work, far-infrared multiple photon dissociation spectroscopy, using Ar as tagging element, is combined with density-functional theory calculations to determine the structures of Aun+ (n≤ 9) clusters formed by laser ablation. While the Au frameworks in Au6Arm+ and Au7Arm+ complexes are confirmed to be planar and that in Au9Arm+ three-dimensional, we demonstrate the coexistence of 3D and planar Au8Arm+ (m = 1-3) isomers. Thus, it is revealed that at finite temperatures, the formal 2D to 3D transition takes place at n = 8 but is not sharp.