Bernard Kirtman

Bio: Bernard Kirtman is an academic researcher from University of California, Santa Barbara. The author has contributed to research in topics: Polarizability & Hyperpolarizability. The author has an hindex of 51, co-authored 228 publications receiving 9950 citations. Previous affiliations of Bernard Kirtman include University of California & Saarland University.

CRYSTAL14: A program for the ab initio investigation of crystalline solids

Abstract: The capabilities of the Crystal14 program are presented, and the improvements made with respect to the previous Crystal09 version discussed. Crystal14 is an ab initio code that uses a Gaussian-type basis set: both pseudopotential and all-electron strategies are permitted; the latter is not much more expensive than the former up to the first-second transition metal rows of the periodic table. A variety of density functionals is available, including as an extreme case Hartree–Fock; hybrids of various nature (global, range-separated, double) can be used. In particular, a very efficient implementation of global hybrids, such as popular B3LYP and PBE0 prescriptions, allows for such calculations to be performed at relatively low computational cost. The program can treat on the same grounds zero-dimensional (molecules), one-dimensional (polymers), two-dimensional (slabs), as well as three-dimensional (3D; crystals) systems. No spurious 3D periodicity is required for low-dimensional systems as happens when plane-waves are used as a basis set. Symmetry is fully exploited at all steps of the calculation; this permits, for example, to investigate nanotubes of increasing radius at a nearly constant cost (better than linear scaling!) or to perform self-consistent-field (SCF) calculations on fullerenes as large as (10,10), with 6000 atoms, 84,000 atomic orbitals, and 20 SCF cycles, on a single core in one day. Three versions of the code exist, serial, parallel, and massive-parallel. In the second one, the most relevant matrices are duplicated, whereas in the third one the matrices in reciprocal space are distributed for diagonalization. All the relevant vectors are now dynamically allocated and deallocated after use, making Crystal14 much more agile than the previous version, in which they were statically allocated. The program now fits more easily in low-memory machines (as many supercomputers nowadays are). Crystal14 can be used on parallel machines up to a high number of cores (benchmarks up to 10,240 cores are documented) with good scalability, the main limitation remaining the diagonalization step. Many tensorial properties can be evaluated in a fully automated way by using a single input keyword: elastic, piezoelectric, photoelastic, dielectric, as well as first and second hyperpolarizabilies, electric field gradients, Born tensors and so forth. Many tools permit a complete analysis of the vibrational properties of crystalline compounds. The infrared and Raman intensities are now computed analytically and related spectra can be generated. Isotopic shifts are easily evaluated, frequencies of only a fragment of a large system computed and nuclear contribution to the dielectric tensor determined. New algorithms have been devised for the investigation of solid solutions and disordered systems. The topological analysis of the electron charge density, according to the Quantum Theory of Atoms in Molecules, is now incorporated in the code via the integrated merge of the Topond package. Electron correlation can be evaluated at the Moller–Plesset second-order level (namely MP2) and a set of double-hybrids are presently available via the integrated merge with the Cryscor program. © 2014 Wiley Periodicals, Inc.

Quantum-mechanical condensed matter simulations with CRYSTAL

Abstract: The latest release of the Crystal program for solid-state quantum-mechanical ab initio simulations is presented. The program adopts atom-centered Gaussian-type functions as a basis set, which makes it possible to perform all-electron as well as pseudopotential calculations. Systems of any periodicity can be treated at the same level of accuracy (from 0D molecules, clusters and nanocrystals, to 1D polymers, helices, nanorods, and nanotubes, to 2D monolayers and slab models for surfaces, to actual 3D bulk crystals), without any artificial repetition along nonperiodic directions for 0–2D systems. Density functional theory calculations can be performed with a variety of functionals belonging to several classes: local-density (LDA), generalized-gradient (GGA), meta-GGA, global hybrid, range-separated hybrid, and self-consistent system-specific hybrid. In particular, hybrid functionals can be used at a modest computational cost, comparable to that of pure LDA and GGA formulations, because of the efficient implementation of exact nonlocal Fock exchange. Both translational and point-symmetry features are fully exploited at all steps of the calculation, thus drastically reducing the corresponding computational cost. The various properties computed encompass electronic structure (including magnetic spin-polarized open-shell systems, electron density analysis), geometry (including full or constrained optimization, transition-state search), vibrational properties (frequencies, infrared and Raman intensities, phonon density of states), thermal properties (quasi-harmonic approximation), linear and nonlinear optical properties (static and dynamic [hyper]polarizabilities), strain properties (elasticity, piezoelectricity, photoelasticity), electron transport properties (Boltzmann, transport across nanojunctions), as well as X-ray and inelastic neutron spectra. The program is distributed in serial, parallel, and massively parallel versions. In this paper, the original developments that have been devised and implemented in the last 4 years (since the distribution of the previous public version, Crystal14, occurred in December 2013) are described.

Assessment of Conventional Density Functional Schemes for Computing the Dipole Moment and (Hyper)polarizabilities of Push−Pull π-Conjugated Systems†

Abstract: DFT schemes based on conventional exchange-correlation (XC) functionals have been employed to determine the dipole moment (μ), polarizability (α), and first (β) and second (γ) hyperpolarizabilities of push−pull π-conjugated systems. In addition to the failures already pointed out for α and γ in a recent study on polyacetylene chains [J. Chem. Phys. 1998, 109, 10489; Phys. Rev. Lett. 1999, 83, 694], these functionals are also unsuitable for the evaluation of μ and β. In the case of β, in particular, an almost catastrophic behavior with respect to increasing chain length is found. We show that the C functional has a negligible effect on the calculated μ, α, β, and γ whereas the X-part is responsible for the large property overestimations when the size of the system increases. The overly large μ values are associated with an overestimation of the charge transfer between the donor and the acceptor whereas for α, β, and γ, incomplete screening of the external electric field is responsible for the large discrep...

Assessment of Conventional Density Functional Schemes for Computing the Polarizabilities and Hyperpolarizabilities of Conjugated Oligomers: An Ab Initio Investigation of Polyacetylene Chains

Abstract: DFT schemes based on conventional and less conventional exchange-correlation (XC) functionals have been employed to determine the polarizability and second hyperpolarizability of π-conjugated polyacetylene chains. These functionals fail in one or more of several ways: (i) the correlation correction to α is either much too small or in the wrong direction, leading to an overestimate; (ii) γ is significantly overestimated; (iii) the chain length dependence is excessively large, particularly for γ and for the more alternant system; and (iv) the bond length alternation effects on γ are either underestimated or qualitatively incorrect. The poor results with the asymptotically correct van Leeuwen–Baerends XC potential show that the overestimations are not related to the asymptotic behavior of the potential. These failures are described in terms of the separate effects of the exchange and the correlation parts of the XC functionals. They are related to the short-sightedness of the XC potentials which are relatively insensitive to the polarization charge induced by the external electric field at the chain ends.

Electric field dependence of the exchange-correlation potential in molecular chains

Abstract: Density functional calculations on the (non)linear optical properties of conjugated molecular chains using currently popular exchange-correlation (xc) potentials give overestimations of several orders of magnitude. By analyzing ``exact'' and Krieger-Li-Iafrate xc potentials, the error is traced back to an incorrect electric field dependence of the ``response part'' of the xc potential in local and gradient-corrected density approximations, which lack a linear term counteracting the applied electric field.

The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals

Abstract: We present two new hybrid meta exchange- correlation functionals, called M06 and M06-2X. The M06 functional is parametrized including both transition metals and nonmetals, whereas the M06-2X functional is a high-nonlocality functional with double the amount of nonlocal exchange (2X), and it is parametrized only for nonmetals.The functionals, along with the previously published M06-L local functional and the M06-HF full-Hartree–Fock functionals, constitute the M06 suite of complementary functionals. We assess these four functionals by comparing their performance to that of 12 other functionals and Hartree–Fock theory for 403 energetic data in 29 diverse databases, including ten databases for thermochemistry, four databases for kinetics, eight databases for noncovalent interactions, three databases for transition metal bonding, one database for metal atom excitation energies, and three databases for molecular excitation energies. We also illustrate the performance of these 17 methods for three databases containing 40 bond lengths and for databases containing 38 vibrational frequencies and 15 vibrational zero point energies. We recommend the M06-2X functional for applications involving main-group thermochemistry, kinetics, noncovalent interactions, and electronic excitation energies to valence and Rydberg states. We recommend the M06 functional for application in organometallic and inorganometallic chemistry and for noncovalent interactions.

Quantum mechanical continuum solvation models.

Abstract: 6.2.2. Definition of Effective Properties 3064 6.3. Response Properties to Magnetic Fields 3066 6.3.1. Nuclear Shielding 3066 6.3.2. Indirect Spin−Spin Coupling 3067 6.3.3. EPR Parameters 3068 6.4. Properties of Chiral Systems 3069 6.4.1. Electronic Circular Dichroism (ECD) 3069 6.4.2. Optical Rotation (OR) 3069 6.4.3. VCD and VROA 3070 7. Continuum and Discrete Models 3071 7.1. Continuum Methods within MD and MC Simulations 3072

Accurate molecular van der Waals interactions from ground-state electron density and free-atom reference data

Abstract: We present a parameter-free method for an accurate determination of long-range van der Waals interactions from mean-field electronic structure calculations. Our method relies on the summation of interatomic C6 coefficients, derived from the electron density of a molecule or solid and accurate reference data for the free atoms. The mean absolute error in the C6 coefficients is 5.5% when compared to accurate experimental values for 1225 intermolecular pairs, irrespective of the employed exchangecorrelation functional. We show that the effective atomic C6 coefficients depend strongly on the bonding environment of an atom in a molecule. Finally, we analyze the van der Waals radii and the damping function in the C6R � 6 correction method for density-functional theory calculations.

A new local density functional for main-group thermochemistry, transition metal bonding, thermochemical kinetics, and noncovalent interactions.

Abstract: We present a new local density functional, called M06-L, for main-group and transition element thermochemistry, thermochemical kinetics, and noncovalent interactions. The functional is designed to capture the main dependence of the exchange-correlation energy on local spin density, spin density gradient, and spin kinetic energy density, and it is parametrized to satisfy the uniform-electron-gas limit and to have good performance for both main-group chemistry and transition metal chemistry. The M06-L functional and 14 other functionals have been comparatively assessed against 22 energetic databases. Among the tested functionals, which include the popular B3LYP, BLYP, and BP86 functionals as well as our previous M05 functional, the M06-L functional gives the best overall performance for a combination of main-group thermochemistry, thermochemical kinetics, and organometallic, inorganometallic, biological, and noncovalent interactions. It also does very well for predicting geometries and vibrational frequencies. Because of the computational advantages of local functionals, the present functional should be very useful for many applications in chemistry, especially for simulations on moderate-sized and large systems and when long time scales must be addressed. © 2006 American Institute of Physics. DOI: 10.1063/1.2370993

Photoinduced electron transfer from a conducting polymer to buckminsterfullerene.

Abstract: Evidence for photoinduced electron transfer from the excited state of a conducting polymer onto buckminsterfullerene, C(60), is reported. After photo-excitation of the conjugated polymer with light of energy greater than the pi-pi* gap, an electron transfer to the C(60) molecule is initiated. Photoinduced optical absorption studies demonstrate a different excitation spectrum for the composite as compared to the separate components, consistent with photo-excited charge transfer. A photoinduced electron spin resonance signal exhibits signatures of both the conducting polymer cation and the C(60) anion. Because the photoluminescence in the conducting polymer is quenched by interaction with C(60), the data imply that charge transfer from the excited state occurs on a picosecond time scale. The charge-separated state in composite films is metastable at low temperatures.