The Emergence of the Hexagonal Lattice in Two-Dimensional Wigner Fragments.

At very low density, the electrons in a uniform electron gas spontaneously break symmetry and form a crystalline lattice called a Wigner crystal. But which type of crystal will the electrons form? We report a numerical study of the density profiles of fragments of Wigner crystals from first principles. To simulate Wigner fragments, we use Clifford periodic boundary conditions and a renormalized distance in the Coulomb potential. Moreover, we show that high-spin restricted open-shell Hartree-Fock theory becomes exact in the low-density limit. We are thus able to accurately capture the localization in two-dimensional Wigner fragments with many electrons. No assumptions about the positions where the electrons will localize are made. The density profiles we obtain emerge naturally when we minimize the total energy of the system. We clearly observe the emergence of the hexagonal crystal structure, which has been predicted to be the ground-state structure of the two-dimensional Wigner crystal.

Development of accurate potentials for the physisorption of water on graphene.

From coupled-cluster singles and doubles model including connected triples corrections [CCSD(T)] calculations on the water dimer and B97D/CC on the water-circumcoronene complex at a large number of randomly generated conformations, interaction potentials for the physisorption of water on graphene are built, accomplishing almost sub-chemical accuracy. The force fields were constructed by decomposing the interaction into electrostatic and van der Waals contributions, the latter represented through improved Lennard-Jones potentials. Besides, a Chemistry at Harvard Macromolecular Mechanics (CHARMM)-like term was included in the water-water potential to improve the description of hydrogen bonds, and an induction term was added to model the polarization effects in the interaction between water and polyaromatic hydrocarbons (PAHs) or graphene. Two schemes with three and six point charges were considered for the interactions water-water and water-PAH, as Coulomb contributions are zero in the water-graphene system. The proposed fitted potentials reproduce the ab initio data used to build them in the whole range of distances and conformations and provide results for selected points very close to CCSD(T) benchmarks. When applied to the water-graphene system, the obtained results are in excellent agreement with p-CCSD(T), revised symmetry-adapted perturbation theory based on density functional theory monomer properties (DFT-SAPT), and diffusion Monte Carlo reference values. Furthermore, the stability of the various conformers water-PAH and water-graphene, as well as the different trends observed between these systems are rationalized in terms of the modifications of the electrostatic contribution.

Grand Canonical Monte Carlo Simulations to Determine the Optimal Interlayer Distance of a Graphene Slit-Shaped Pore for Adsorption of Methane, Hydrogen and their Equimolar Mixture.

The adsorption-for separation, storage and transportation-of methane, hydrogen and their mixture is important for a sustainable energy consumption in present-day society. Graphene derivatives have proven to be very promising for such an application, yet for a good design a better understanding of the optimal pore size is needed. In this work, grand canonical Monte Carlo simulations, employing Improved Lennard-Jones potentials, are performed to determine the ideal interlayer distance for a slit-shaped graphene pore in a large pressure range. A detailed study of the adsorption behavior of methane, hydrogen and their equimolar mixture in different sizes of graphene pores is obtained through calculation of absolute and excess adsorption isotherms, isosteric heats and the selectivity. Moreover, a molecular picture is provided through z-density profiles at low and high pressure. It is found that an interlayer distance of about twice the van der Waals distance of the adsorbate is recommended to enhance the adsorbing ability. Furthermore, the graphene structures with slit-shaped pores were found to be very capable of adsorbing methane and separating methane from hydrogen in a mixture at reasonable working conditions (300 K and well below 15 atm).

Molecular Dynamics of CH4/N2 Mixtures on a Flexible Graphene Layer: Adsorption and Selectivity Case Study.

We theoretically investigate graphene layers, proposing them as membranes of subnanometer size suitable for CH4/N2 separation and gas uptake. The observed potential energy surfaces, representing the intermolecular interactions within the CH4/N2 gaseous mixtures and between these and the graphene layers, have been formulated by adopting the so-called Improved Lennard-Jones (ILJ) potential, which is far more accurate than the traditional Lennard-Jones potential. Previously derived ILJ force fields are used to perform extensive molecular dynamics simulations on graphene's ability to separate and adsorb the CH4/N2 mixture. Furthermore, the intramolecular interactions within graphene were explicitly considered since they are responsible for its flexibility and the consequent out-of-plane movements of the constituting carbon atoms. The effects on the adsorption capacity of graphene caused by introducing its flexibility in the simulations are assessed via comparison of different intramolecular force fields giving account of flexibility against a simplified less realistic model that considers graphene to be rigid. The accuracy of the potentials guarantees a quantitative description of the interactions and trustable results for the dynamics, as long as the appropriate set of intramolecular and intermolecular force fields is chosen. In particular it is shown that only if the flexibility of graphene is explicitly taken into account, a simple united-atom interaction potential can provide correct predictions. Conversely, when using an atomistic model, neglecting in the simulations the intrinsic flexibility of the graphene sheet has a minor effect. From a practical point of view, the global analysis of the whole set of results proves that these nanostructures are versatile materials competitive with other carbon-based adsorbing membranes suitable to cope with CH4 and N2 adsorption.

Modeling the Interaction of Carbon Monoxide with Flexible Graphene: From Coupled Cluster Calculations to Molecular-Dynamics Simulations.

The interaction of CO with graphene was studied at different theoretical levels. Quantum-mechanical calculations on finite graphene models with the use of coronene for coupled cluster calculations and circumcoronene for B97D calculations showed that there was no preferential site for adsorption and that the most important factor was the orientation of CO relative to graphene. The parallel orientation was preferred, with binding energies around 9 kJ mol-1 at the CCSD(T) and B97D levels, which was in good agreement with experimental findings. From a large number of CO-circumcoronene and CO-CO interactions, computed at different distances and randomly generated orientations, parameters were fit to the improved Lennard-Jones potential. Such potentials, together with others describing the intramolecular dynamics of graphene, were subsequently employed in classical molecular-dynamics simulations of the adsorption of CO on graphene by using the canonical ensemble. The obtained results showed that the introduction of flexibility in graphene, which simulated the effects associated to curvature of the surface, diminished the adsorption level and that, as expected, adsorption also diminished with temperature.

π-Ring currents in doped coronenes with nitrogen and boron: diatropic-paratropic duality.

The change in the electronic structure of coronene upon doping with nitrogen or boron has been theoretically studied by means of its magnetic properties and magnetic field induced current density maps. The addition of two atoms of nitrogen or boron to the central ring of coronene causes a drastic variation in the delocalization of π-electrons, which does not depend on its nature but instead on its position. Then, doping in the para position makes coronene more aromatic while doping in the meta position makes it to become antiaromatic. The magnetic behavior of the pristine molecule is characterized by two concentric currents flowing in opposite senses that are converted into hemi-perimetric currents in the ortho and meta isomers, so dividing the molecule into aromatic and antiaromatic regions. The paratropic and diatropic ring currents of the coronene moiety may, therefore, be modulated through the position of the heteroatom and, consequently, also the localized/delocalized behavior.

Multi-level coupled cluster theory.

We present a general formalism where different levels of coupled cluster theory can be applied to different parts of the molecular system. The system is partitioned into subsystems by Cholesky decomposition of the one-electron Hartree-Fock density matrix. In this way the system can be divided across chemical bonds without discontinuities arising. The coupled cluster wave function is defined in terms of cluster operators for each part and these are determined from a set of coupled equations. The total wave function fulfills the Pauli-principle across all borders and levels of electron correlation. We develop the associated response theory for this multi-level coupled cluster theory and present proof of principle applications. The formalism is an essential tool in order to obtain size-intensive complexity in the calculation of local molecular properties.

Calculation of excitation energies from the CC2 linear response theory using Cholesky decomposition.

A new implementation of the approximate coupled cluster singles and doubles CC2 linear response model is reported. It employs a Cholesky decomposition of the two-electron integrals that significantly reduces the computational cost and the storage requirements of the method compared to standard implementations. Our algorithm also exploits a partitioning form of the CC2 equations which reduces the dimension of the problem and avoids the storage of doubles amplitudes. We present calculation of excitation energies of benzene using a hierarchy of basis sets and compare the results with conventional CC2 calculations. The reduction of the scaling is evaluated as well as the effect of the Cholesky decomposition parameter on the quality of the results. The new algorithm is used to perform an extrapolation to complete basis set investigation on the spectroscopically interesting benzylallene conformers. A set of calculations on medium-sized molecules is carried out to check the dependence of the accuracy of the results on the decomposition thresholds. Moreover, CC2 singlet excitation energies of the free base porphin are also presented.

Energy interactions in amyloid-like fibrils from NNQQNY.

We use large-scale MP2 calculations to analyze the interactions appearing in amyloid fibers, which are difficult to determine experimentally. To this end, dimers and trimers of the hexapeptide NNQQNY from the yeast prion-like protein Sup35 were considered as model systems. We studied the energy interactions present in the three levels of organization in which the formation of amyloid fibrils is structured. The structural changes in the hydrogen bonds were studied too. It was found that the most energetic process is the formation of the ß-sheet, which is equally due to both hydrogen bonds and van der Waals interactions. The aromatic rings help stabilize these aggregates through stacking of the aromatic rings of tyrosine, the stability produced by the aromatics residues increasing with their aromaticity. The formation of the basic unit of the assembled proto-fiber, the steric zipper, is less energetic and is associated to both dispersion forces and hydrogen bonds. The interactions between pair of ß-sheets across the peptide-to-peptide contact through the tyrosine rings are cooperative and due to dispersion effects. Moreover, the strength of this interaction can rationalize the variation of mobility of the aromatic ring in the tyrosine units found in solid NMR experiments.

Cholesky decomposition-based definition of atomic subsystems in electronic structure calculations.

Decomposing the Hartree-Fock one-electron density matrix and a virtual pseudodensity matrix, we obtain an orthogonal set of normalized molecular orbitals with local character to be used in post-Hartree-Fock calculations. The applicability of the procedure is illustrated by calculating CCSD(T) energies and CCSD molecular properties in reduced active spaces.

How nitrogen modifies the nuclear magnetic shielding in tetraazanaphthalenes.

Although for planar conjugated hydrocarbons the out-of-plane component of proton magnetic shielding is an unquestionable quantitative aromaticity indicator, the same is not true for tetraazanaphthalenes. As in these compounds the (core + sigma)-currents associated to the nitrogen nuclei diminish the perpendicular component of shielding, abnormal values of (1)H NMR sigma(zz) are obtained. Therefore, a consistent aromaticity measure must be based only on the pi-contribution to the out-of-plane component of proton magnetic shielding. Otherwise, the behavior of these compounds in presence of an external magnetic field parallels that of naphthalene, with the nitrogen nuclei contributing to the ring current in a comparable amount to carbon nuclei. The pi-current contribution to magnetic shielding represents 6-8% of the out-of-plane shielding for nitrogen and 9-12% for carbon.

Method specific Cholesky decomposition: coulomb and exchange energies.

We present a novel approach to the calculation of the Coulomb and exchange contributions to the total electronic energy in self consistent field and density functional theory. The numerical procedure is based on the Cholesky decomposition and involves decomposition of specific Hadamard product matrices that enter the energy expression. In this way, we determine an auxiliary basis and obtain a dramatic reduction in size as compared to the resolution of identity (RI) method. Although the auxiliary basis is determined from the energy expression, we have complete control of the errors in the gradient or Fock matrix. Another important advantage of this method specific Cholesky decomposition is that the exchange energy and Fock matrix can be evaluated with a linear scaling effort contrary to the RI method or standard Cholesky decomposition of the two-electron integral matrix. The methods presented show the same scaling properties as the so-called local density fitting methods, but with full error control.

Accurate ab initio density fitting for multiconfigurational self-consistent field methods.

Using Cholesky decomposition and density fitting to approximate the electron repulsion integrals, an implementation of the complete active space self-consistent field (CASSCF) method suitable for large-scale applications is presented. Sample calculations on benzene, diaquo-tetra-mu-acetato-dicopper(II), and diuraniumendofullerene demonstrate that the Cholesky and density fitting approximations allow larger basis sets and larger systems to be treated at the CASSCF level of theory with controllable accuracy. While strict error control is an inherent property of the Cholesky approximation, errors arising from the density fitting approach are managed by using a recently proposed class of auxiliary basis sets constructed from Cholesky decomposition of the atomic electron repulsion integrals.

Electronic structure of the ground and excited states of beta-carboline.

Coupled-cluster calculations are used to compute the energy of conversion between the neutral and the zwitterionic forms of beta-carboline. The stability of the different species is discussed in terms of charge separation and aromatic character, which is related to magnetic criteria. By means of a linear response formalism the vertical excitation energies and oscillator strengths of the lowest singlet states of both structures as well as of the cationic species are determined. General agreement of the relative position and intensity of the different peaks with experimental data is achieved, but the overall spectra are slightly displaced because of solvent effects.

Variation of polarizability in the [4n+2] annulene series: from [22]- to [66]-annulene.

Using correlated ab initio methods, the polarizability of large [4n+2]-annulenes is determined, showing that there exists an almost linear relation between the exaltation of magnetic susceptibility (a measure of aromaticity) and an equivalent enlargement of polarizability.

Fast noniterative orbital localization for large molecules.

We use Cholesky decomposition of the density matrix in atomic orbital basis to define a new set of occupied molecular orbital coefficients. Analysis of the resulting orbitals ("Cholesky molecular orbitals") demonstrates their localized character inherited from the sparsity of the density matrix. Comparison with the results of traditional iterative localization schemes shows minor differences with respect to a number of suitable measures of locality, particularly the scaling with system size of orbital pair domains used in local correlation methods. The Cholesky procedure for generating orthonormal localized orbitals is noniterative and may be made linear scaling. Although our present implementation scales cubically, the algorithm is significantly faster than any of the conventional localization schemes. In addition, since this approach does not require starting orbitals, it will be useful in local correlation treatments on top of diagonalization-free Hartree-Fock optimization algorithms.

Lower Rydberg series of methane: a combined coupled cluster linear response and molecular quantum defect orbital calculation.

Vertical excitation energies as well as related absolute photoabsorption oscillator strength data are very scarce in the literature for methane. In this study, we have characterized the three existing series of low-lying Rydberg states of CH4 by computing coupled cluster linear response (CCLR) vertical excitation energies together with oscillator strengths in the molecular-adapted quantum defect orbital formalism from a distorted Cs geometry selected on the basis of outer valence green function calculations. The present work provides a wide range of data of excitation energies and absolute oscillator strengths which correspond to the Rydberg series converging to the three lower ionization potential values of the distorted methane molecule, in energy regions for which experimentally measured data appear to be unavailable.

Theoretical study on the structures and electronic spectra of TCNE2-.

Investigations into the charge-separated states and electron-transfer transitions in tetracyanoethylene (TCNE) complexes have recently generated much interest. In this work we present theoretical calculations showing that the most stable structure of the dianion TCNE2- has D2d symmetry in vacuum as well as in the solvents dichloromethane and acetonitrile. By means of the coupled cluster linear response, we compute the vertical electronic spectrum in both the gas phase and solution. The theoretical results are compared to the experimental data and good agreement is achieved.

Structure, magnetizability, and nuclear magnetic shielding tensors of bis-heteropentalenes. IV. Dihydrophospholophosphole isomers.

The geometry of the heteropentalenes formed by two phosphole units has been determined at the DFT level. The magnetic susceptibility and the nuclear magnetic shielding at the nuclei of these systems have also been calculated using gauge-including atomic orbitals and a large Gaussian basis set to achieve near Hartree-Fock estimates. A comparative study of the various isomers, of their flattened analogs, and of the parent phosphole molecule, shows that the [3,4-c] isomer is the most aromatic system in the set considered, assuming diatropicity and degree of planarity as indicators, even if it is the less stable in terms of total molecular energy. Plots of magnetic field-induced current densities confirm diatropicity of P-containing bis-heteropentalenenes, showing, however, significant differences from the analogous systems with distinct heteroatoms. The maps give evidence of spiral flow nearby C-C bonds, compatible with prevalent distortive behavior of pi electrons exalted by pyramidalization at P, and competing against the sigma electron compression, which would favor planar structure.

From pentalene to dicyclopenta[b,g]naphthalene, or the change towards delocalized structures.

Using triples-corrected coupled-cluster methods as well as other high-level theoretical approximations, the optimized parameters and isomerization barriers of the family of compounds cyclopentadiene-(benzene)x-cyclopentadiene (x = 0, 1, 2) are computed. In contrast to previous studies, s-indacene presents a localized C(2h) geometry. Also, the localized structure of pentalene is found to be the most stable, but when two benzene rings are intercalated between the five-member rings of pentalene, the resulting molecule preferably adopts a delocalized D(2h) conformation.