Theoretical study of finite temperature spectroscopy in van der Waals clusters. II Time-dependent absorption spectra

Using approximate partition functions and a master equation approach, we investigate the statistical relaxation toward equilibrium in selected CaAr$_n$ clusters. The Gaussian theory of absorption (previous article) is employed to calculate the average photoabsorption intensity associated with the 4s^2-> 4s^14p^1 transition of calcium as a function of time during relaxation. In CaAr_6 and CaAr_10 simple relaxation is observed with a single time scale. CaAr_13 exhibits much slower dynamics and the relaxation occurs over two distinct time scales. CaAr_37 shows much slower relaxation with multiple transients, reminiscent of glassy behavior due to competition between different low-energy structures. We interpret these results in terms of the underlying potential energy surfaces for these clusters.Comment: 10 pages, 9 figure

A Poset Connected to Artin Monoids of Simply Laced Type

Let W be a Weyl group whose type is a simply laced Dynkin diagram. On several W-orbits of sets of mutually commuting reflections, a poset is described which plays a role in linear representatons of the corresponding Artin group A. The poset generalizes many properties of the usual order on positive roots of W given by height. In this paper, a linear representation of the positive monoid of A is defined by use of the poset

Understanding fragility in supercooled Lennard-Jones mixtures. II. Potential energy surface

We numerically investigated the connection between isobaric fragility and the properties of high-order stationary points of the potential energy surface in different supercooled Lennard-Jones mixtures. The increase of effective activation energies upon supercooling appears to be driven by the increase of average potential energy barriers measured by the energy dependence of the fraction of unstable modes. Such an increase is sharper, the more fragile is the mixture. Correlations between fragility and other properties of high-order stationary points, including the vibrational density of states and the localization features of unstable modes, are also discussed.Comment: 13 pages, 13 figures, minor revisions, one figure adde

Communication: optimal parameters for basin-hopping global optimization based on Tsallis statistics.

A fundamental problem associated with global optimization is the large free energy barrier for the corresponding solid-solid phase transitions for systems with multi-funnel energy landscapes. To address this issue we consider the Tsallis weight instead of the Boltzmann weight to define the acceptance ratio for basin-hopping global optimization. Benchmarks for atomic clusters show that using the optimal Tsallis weight can improve the efficiency by roughly a factor of two. We present a theory that connects the optimal parameters for the Tsallis weighting, and demonstrate that the predictions are verified for each of the test cases.This work was supported by the ERC and the EPSRC.This is the accepted manuscript version of the article. Copyright 2014 American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics. The following article appeared in J. Chem. Phys. 141, 071101 (2014) and may be found at http://dx.doi.org/10.1063/1.4893344

Non existence of a phase transition for the Penetrable Square Wells in one dimension

Penetrable Square Wells in one dimension were introduced for the first time in [A. Santos et. al., Phys. Rev. E, 77, 051206 (2008)] as a paradigm for ultra-soft colloids. Using the Kastner, Schreiber, and Schnetz theorem [M. Kastner, Rev. Mod. Phys., 80, 167 (2008)] we give strong evidence for the absence of any phase transition for this model. The argument can be generalized to a large class of model fluids and complements the van Hove's theorem.Comment: 14 pages, 7 figures, 1 tabl

Quantum annealing of the Traveling Salesman Problem

We propose a path-integral Monte Carlo quantum annealing scheme for the symmetric Traveling Salesman Problem, based on a highly constrained Ising-like representation, and we compare its performance against standard thermal Simulated Annealing. The Monte Carlo moves implemented are standard, and consist in restructuring a tour by exchanging two links (2-opt moves). The quantum annealing scheme, even with a drastically simple form of kinetic energy, appears definitely superior to the classical one, when tested on a 1002 city instance of the standard TSPLIB.Comment: 5 pages, 2 figure

Protein Structure Prediction Using Basin-Hopping

Associative memory Hamiltonian structure prediction potentials are not overly rugged, thereby suggesting their landscapes are like those of actual proteins. In the present contribution we show how basin-hopping global optimization can identify low-lying minima for the corresponding mildly frustrated energy landscapes. For small systems the basin-hopping algorithm succeeds in locating both lower minima and conformations closer to the experimental structure than does molecular dynamics with simulated annealing. For large systems the efficiency of basin-hopping decreases for our initial implementation, where the steps consist of random perturbations to the Cartesian coordinates. We implemented umbrella sampling using basin-hopping to further confirm when the global minima are reached. We have also improved the energy surface by employing bioinformatic techniques for reducing the roughness or variance of the energy surface. Finally, the basin-hopping calculations have guided improvements in the excluded volume of the Hamiltonian, producing better structures. These results suggest a novel and transferable optimization scheme for future energy function development

Stacked clusters of polycyclic aromatic hydrocarbon molecules

Clusters of polycyclic aromatic hydrocarbon (PAH) molecules are modelled using explicit all-atom potentials using a rigid body approximation. The PAH's considered range from pyrene (C10H8) to circumcoronene (C54H18), and clusters containing between 2 and 32 molecules are investigated. In addition to the usual repulsion-dispersion interactions, electrostatic point-charge interactions are incorporated, as obtained from density functional theory calculations. The general electrostatic distribution in neutral or singly charged PAH's is reproduced well using a fluctuating charges analysis, which provides an adequate description of the multipolar distribution. Global optimization is performed using a variety of methods, including basin-hopping and parallel tempering Monte Carlo. We find evidence that stacking the PAH molecules generally yields the most stable motif. A structural transition between one-dimensional stacks and three-dimensional shapes built from mutiple stacks is observed at larger sizes, and the threshold for this transition increases with the size of the monomer. Larger aggregates seem to evolve toward the packing observed for benzene in bulk.Difficulties met in optimizing these clusters are analysed in terms of the strong anisotropy of the molecules. We also discuss segregation in heterogeneous clusters and vibrational properties in the context of astrophysical observations.Comment: 12 pages, 7 figure