First-Principle Calculations of Large Fullerenes

Magnetic response properties of carbon nano-onions

The magnetic response of a number of double- and triple-layer carbon nano-onions (CNOs) is analyzed by calculating the magnetically induced current density and the induced magnetic field using the pseudo-π model. Qualitatively the same magnetic response was obtained in calculations at the all-electron level. The calculations show that the CNOs exhibit strong net diatropic (paratropic) ring currents when the external magnetic field points perpendicularly to one of the six-membered (five-membered) rings. They are deshielded inside and shielded outside the CNO; the latter dominates for larger CNOs. The magnetic response originates from a combination of spherical paratropic current densities on the inside of each carbon layer and diatropic current densities on the outside of them. The quantitative differences in the aromaticity of the CNOs as compared to single fullerenes are discussed in terms of ring-current strengths. The magnetic response of some of the CNOs is approximately the sum of the magnetic response of the individual layers, whereas deviations are significant for CNOs containing C₈₀. For the largest CNOs, the deviation from the sum of the fullerene contributions is larger, especially when the external magnetic field is perpendicular to a six-membered ring.

Self-Consistent Auxiliary Density Perturbation Theory

The working equations for the extension of auxiliary density perturbation theory (ADPT) to hybrid functionals, employing the variational fitting of the Fock potential, are derived. The response equations in the resulting self-consistent ADPT (SC-ADPT) are solved iteratively with an adapted Eirola-Nevanlinna algorithm. As a result, a memory and CPU time efficient implementation of perturbation theory free of four-center electron repulsion integrals (ERIs) is obtained. Our validation calculations of SC-ADPT static and dynamic polarizabilities show quantitative agreement with corresponding coupled perturbed Hartree-Fock and Kohn-Sham results employing four-center ERIs. The comparison of SC-ADPT hybrid functional polarizabilities with coupled cluster reference calculations yield semiquantitative agreement. The presented systematic study of the dynamic polarizabilities of oligothiophenes shows that hybrid functionals can overcome the pathological misplacement of excitation poles by the local density and generalized gradient approximations. Good agreement with experimental dynamic polarizabilities for all studied oligothiophenes is achieved with range-separated hybrid functionals in the framework of SC-ADPT.

Buckling transitions and soft-phase invasion of two-component icosahedral shells

What is the optimal distribution of two types of crystalline phases on the surface of icosahedral shells, such as of many viral capsids? We here investigate the distribution of a thin layer of soft material on a crystalline convex icosahedral shell. We demonstrate how the shapes of spherical viruses can be understood from the perspective of elasticity theory of thin two-component shells. We develop a theory of shape transformations of an icosahedral shell upon addition of a softer, but still crystalline, material onto its surface. We show how the soft component "invades" the regions with the highest elastic energy and stress imposed by the 12 topological defects on the surface. We explore the phase diagram as a function of the surface fraction of the soft material, the shell size, and the incommensurability of the elastic moduli of the rigid and soft phases. We find that, as expected, progressive filling of the rigid shell by the soft phase starts from the most deformed regions of the icosahedron. With a progressively increasing soft-phase coverage, the spherical segments of domes are filled first (12 vertices of the shell), then the cylindrical segments connecting the domes (30 edges) are invaded, and, ultimately, the 20 flat faces of the icosahedral shell tend to be occupied by the soft material. We present a detailed theoretical investigation of the first two stages of this invasion process and develop a model of morphological changes of the cone structure that permits noncircular cross sections. In conclusion, we discuss the biological relevance of some structures predicted from our calculations, in particular for the shape of viral capsids.

Bootstrap Embedding For Large Molecular Systems

Recent developments in quantum embedding theories have provided attractive approaches to correlated calculations for large systems. In this work, we extend our previous work [J. Chem. Theory Comput. 2019, 15, 4497-4506; J. Phys. Chem. Lett. 2019, 10, 6368-6374] on bootstrap embedding (BE) to enable correlated ab initio calculations at the coupled cluster with singles and doubles (CCSD) level for large molecules. We introduce several new algorithmic developments that significantly reduce the computational cost of BE, while maintaining its accuracy. The resulting implementation scales as O(N³) for the integral transform and O(N) for the CCSD calculation. Numerical results on a series of conjugated molecules suggest that BE with reasonably sized fragments can recover more than 99.5% of the total correlation energy of a full CCSD calculation, while the required computational resources (time and storage) compare favorably to one popular local correlation scheme: domain localized pair natural orbital (DLPNO). The largest BE calculation in this work involves ∼2900 basis functions and can be performed on a single node with 16 CPU cores and 64 GB of memory in a few days. We anticipate that these developments represent an important step toward the application of BE to solve practical problems.

Hybrid ADFT Study of the C104 and C106 IPR Isomers

This work presents a hybrid auxiliary density functional theory (ADFT) study of the neutral and hexaanionic C₁₀₄ and C₁₀₆ fullerenes with the aim to determine their ground state structures. To this end, all C₁₀₄ and C₁₀₆ fullerene structures that obey the isolated pentagon rule (IPR) were optimized with the Perdew-Burke-Ernzerhof generalized gradient approximation followed by a single-point energy calculation with the PBE0 hybrid functional. Our studies show that this composite approach yields relative energies of giant fullerenes that are accurate to around 1 kcal/mol. As a result, the ground states of C₁₀₄, C₁₀₄^6-, and C₁₀₆^6- can be assigned to the isomers 234:C_s, 821:D₂, and 891:C_s, respectively. On the other hand, the energetically lowest lying IPR isomers of C₁₀₆, 331:C_s, 1194:C₂, 534:C₁ are separated by less than 1 kcal/mol which makes an unequivocal ground state assignment by hybrid DFT methods impossible. To guide future experiments, we also report the simulated IR and Raman spectra of the most stable neutral and hexaanionic C₁₀₄ and C₁₀₆ fullerenes.

Generalized structural motif model for studying the thermodynamic stability of fullerenes: from C60to graphene passing through giant fullerenes

A generalized motif model to describe the stability of neutral fullerenes, covering the full range of cage sizes, starting from C₆₀, going through giant fullerenes, and ultimately leading to graphene.

Efficient calculation of nuclear spin-rotation constants from auxiliary density functional theory

The computation of the spin-rotation tensor within the framework of auxiliary density functional theory (ADFT) in combination with the gauge including atomic orbital (GIAO) scheme, to treat the gauge origin problem, is presented. For the spin-rotation tensor, the calculation of the magnetic shielding tensor represents the most demanding computational task. Employing the ADFT-GIAO methodology, the central processing unit time for the magnetic shielding tensor calculation can be dramatically reduced. In this work, the quality of spin-rotation constants obtained with the ADFT-GIAO methodology is compared with available experimental data as well as with other theoretical results at the Hartree-Fock and coupled-cluster level of theory. It is found that the agreement between the ADFT-GIAO results and the experiment is good and very similar to the ones obtained by the coupled-cluster single-doubles-perturbative triples-GIAO methodology. With the improved computational performance achieved, the computation of the spin-rotation tensors of large systems or along Born-Oppenheimer molecular dynamics trajectories becomes feasible in reasonable times. Three models of carbon fullerenes containing hundreds of atoms and thousands of basis functions are used for benchmarking the performance. Furthermore, a theoretical study of temperature effects on the structure and spin-rotation tensor of the H(12)C-(12)CH-DF complex is presented. Here, the temperature dependency of the spin-rotation tensor of the fluorine nucleus can be used to identify experimentally the so far unknown bent isomer of this complex. To the best of our knowledge this is the first time that temperature effects on the spin-rotation tensor are investigated.

Efficient calculation of the rotational g tensor from auxiliary density functional theory

The computation of the rotational g tensor with the recently developed auxiliary density functional theory (ADFT) gauge including atomic orbital (GIAO) methodology is presented. For the rotational g tensor, the calculation of the magnetizability tensor represents the most demanding computational task. With the ADFT-GIAO methodology, the CPU time for the magnetizability tensor calculation can be dramatically reduced. Therefore, it seems most desirable to employ the ADFT-GIAO methodology also for the computation of the rotational g tensor. In this work, the quality of rotational g tensors obtained with the ADFT-GIAO methodology is compared with available experimental data as well as with other theoretical results at the Hartree-Fock and coupled-cluster level of theory. It is found that the agreement between the ADFT-GIAO results and the experiment is good. Furthermore, we also show that the ADFT-GIAO g tensor calculation is applicable to large systems like carbon nanotube models containing hundreds of atom and thousands of basis functions.

Structural, electronic and energetic properties of giant icosahedral fullerenes up to C6000: insights from an ab initio hybrid DFT study

An accurate and detailed electronic structure study of giant carbon fullerenes, benefitting from improved symmetry exploitation in the CRYSTAL14 code.