Aggregation patterns of curcumin and piperine mixtures in different polar media.

This work reports a thorough molecular dynamics investigation on the aggregation patterns of curcumin and piperine in water, ethanol and a mixture of both solvents. The low solubility of curcumin in water results in a rapid formation of very stable dimers for both keto and enol tautomers. In agreement with a higher solubility, piperine molecules move closer and farther apart several times during the simulation, which indicates the formation of a less stable dimer in water. In contrast, both curcumin and piperine are soluble in ethanol and, thus, dimers can hardly be formed in this media. In comparison with a pure-water solvent, a 30 : 70 mixture of ethanol and water significantly reduces the probability of formation of most dimers of curcumin and piperine molecules. The simulations show that larger clusters may be complex structures, but the formation of stacks (in the case of piperine and enol tautomer of curcumin) and cages (when the keto tautomer of curcumin is involved) are not rare. Furthermore, it is shown that each single molecule presents a certain degree of mobility in the cluster, especially on the surface, but without leading to dissociation.

A global optimization perspective on molecular clusters.

Although there is a long history behind the idea of chemical structure, this is a key concept that continues to challenge chemists. Chemical structure is fundamental to understanding most of the properties of matter and its knowledge for complex systems requires the use of state-of-the-art techniques, either experimental or theoretical. From the theoretical view point, one needs to establish the interaction potential among the atoms or molecules of the system, which contains all the information regarding the energy landscape, and employ optimization algorithms to discover the relevant stationary points. In particular, global optimization methods are of major importance to search for the low-energy structures of molecular aggregates. We review the application of global optimization techniques to several molecular clusters; some new results are also reported. Emphasis is given to evolutionary algorithms and their application in the study of the microsolvation of alkali-metal and Ca2+ ions with various types of solvents.This article is part of the themed issue 'Theoretical and computational studies of non-equilibrium and non-statistical dynamics in the gas phase, in the condensed phase and at interfaces'.

On the use of big-bang method to generate low-energy structures of atomic clusters modeled with pair potentials of different ranges.

The efficiency of the so-called big-bang method for the optimization of atomic clusters is analysed in detail for Morse pair potentials with different ranges; here, we have used Morse potentials with four different ranges, from long- ρ = 3) to short-ranged ρ = 14) interactions. Specifically, we study the efficacy of the method in discovering low-energy structures, including the putative global minimum, as a function of the potential range and the cluster size. A new global minimum structure for long-ranged ρ = 3) Morse potential at the cluster size of n= 240 is reported. The present results are useful to assess the maximum cluster size for each type of interaction where the global minimum can be discovered with a limited number of big-bang trials.

How different are two chemical structures?

We extend the scope of a recent method for superimposing two molecules ( J. Chem. Phys. 2009, 131, 124126-1-124126-10 ) to include the identification of chiral structures. This methodology is tested by applying it to several organic molecules and water clusters that were subjected to geometry optimization. The accuracy of four simpler, non-superimposing approaches is then analyzed by comparing pairs of structures for argon and water clusters. The structures considered in this work were obtained by a Markovian walk in the coordinate space. First, a random geometry is generated, and then, the iterative application of a mutation operator ensures the creation of increasingly dissimilar structures. The discriminating power of the non-superimposing approaches is tested by comparing the corresponding dissimilarity measures with the root-mean-square distance obtained from the superimposing method. Finally, we showcase the application of those methods to characterize the diversity of solutions in global geometry optimization by evolutionary algorithms.

Generation and characterization of low-energy structures in atomic clusters.

Factors relevant for controlling the structures determined in the local optimization of argon clusters are investigated. In particular, the role of volume and shape for the box where initial structures are generated is assessed. A thorough characterization of the optimization is also presented, based on a nearest-neighbor analysis, in clusters ranging from 30 to 55 atoms. This includes the assessment of the degree of preservation of aspects of the initial randomly generated structure in the final optimized counterpart, and the correlation between optimized energy and the number of nearest neighbors and average departure from the diatomic reference distance. The usefulness of this analysis to explore the energy landscape of atomic clusters is also highlighted.

Structure Dependence of Hyperpolarizability in Octopolar Molecules.

Recent Hyper Rayleigh Scattering measurements report a significant increase of second-order hyperpolarizability upon introduction of positive charges at the pyridyl end groups in trispyridyl octopolar chromophores. We calculated the geometries, linear response, and first-order hyperpolarizabilities of a series of six trispyridyl molecules both in the neutral and protonated forms. The calculations were performed with ab initio and semiempirical methods. The results are in good agreement with the experimental values and a correlation between the first hyperpolarizability and two structural properties, the N-C bond elongation and the C-C bond length alternation, ⟨Δr⟩ was established. To test these effects we computed the hyperpolarizability for several constrained geometries and confirmed the importance of planarity on the hyperpolarizability values. However the ⟨Δr⟩ values alone seem to have little influence both on the hyperpolarizability and on the gap values. Replacing the triple C≡C bond by a double C=C bond in the conjugation bridge has no significant effect due to the strong hyperpolarizability dependence on the pyridyl-benzene dihedral angle.