Computational Physics on Graphics Processing Units

The use of graphics processing units for scientific computations is an emerging strategy that can significantly speed up various different algorithms. In this review, we discuss advances made in the field of computational physics, focusing on classical molecular dynamics, and on quantum simulations for electronic structure calculations using the density functional theory, wave function techniques, and quantum field theory.Comment: Proceedings of the 11th International Conference, PARA 2012, Helsinki, Finland, June 10-13, 201

Evolution of Linear Absorption and Nonlinear Optical Properties in V-Shaped Ruthenium(II)-Based Chromophores

In this article, we describe a series of complexes with electron-rich cis-{Ru^(II)(NH_3)_4}^(2+) centers coordinated to two pyridyl ligands bearing N-methyl/arylpyridinium electron-acceptor groups. These V-shaped dipolar species are new, extended members of a class of chromophores first reported by us (Coe, B. J. et al. J. Am. Chem. Soc. 2005, 127, 4845−4859). They have been isolated as their PF_6− salts and characterized by using various techniques including ^1H NMR and electronic absorption spectroscopies and cyclic voltammetry. Reversible Ru^(III/II) waves show that the new complexes are potentially redox-switchable chromophores. Single crystal X-ray structures have been obtained for four complex salts; three of these crystallize noncentrosymmetrically, but with the individual molecular dipoles aligned largely antiparallel. Very large molecular first hyperpolarizabilities β have been determined by using hyper-Rayleigh scattering (HRS) with an 800 nm laser and also via Stark (electroabsorption) spectroscopic studies on the intense, visible d → π^* metal-to-ligand charge-transfer (MLCT) and π → π^* intraligand charge-transfer (ILCT) bands. The latter measurements afford total nonresonant β_0 responses as high as ca. 600 × 10^(−30) esu. These pseudo-C_(2v) chromophores show two substantial components of the β tensor, β_(zzz) and β_(zyy), although the relative significance of these varies with the physical method applied. According to HRS, β_(zzz) dominates in all cases, whereas the Stark analyses indicate that β_(zyy) is dominant in the shorter chromophores, but β_(zzz) and β_(zyy) are similar for the extended species. In contrast, finite field calculations predict that β_(zyy) is always the major component. Time-dependent density functional theory calculations predict increasing ILCT character for the nominally MLCT transitions and accompanying blue-shifts of the visible absorptions, as the ligand π-systems are extended. Such unusual behavior has also been observed with related 1D complexes (Coe, B. J. et al. J. Am. Chem. Soc. 2004, 126, 3880−3891)

Diquat Derivatives: Highly Active, Two-Dimensional Nonlinear Optical Chromophores with Potential Redox Switchability

In this article, we present a detailed study of structure−activity relationships in diquaternized 2,2′-bipyridyl (diquat) derivatives. Sixteen new chromophores have been synthesized, with variations in the amino electron donor substituents, π-conjugated bridge, and alkyl diquaternizing unit. Our aim is to combine very large, two-dimensional (2D) quadratic nonlinear optical (NLO) responses with reversible redox chemistry. The chromophores have been characterized as their PF_6^− salts by using various techniques including electronic absorption spectroscopy and cyclic voltammetry. Their visible absorption spectra are dominated by intense π → π^* intramolecular charge-transfer (ICT) bands, and all show two reversible diquat-based reductions. First hyperpolarizabilities β have been measured by using hyper-Rayleigh scattering with an 800 nm laser, and Stark spectroscopy of the ICT bands affords estimated static first hyperpolarizabilities β_0. The directly and indirectly derived β values are large and increase with the extent of π-conjugation and electron donor strength. Extending the quaternizing alkyl linkage always increases the ICT energy and decreases the E_(1/2) values for diquat reduction, but a compensating increase in the ICT intensity prevents significant decreases in Stark-based β_0 responses. Nine single-crystal X-ray structures have also been obtained. Time-dependent density functional theory clarifies the molecular electronic/optical properties, and finite field calculations agree with polarized HRS data in that the NLO responses of the disubstituted species are dominated by ‘off-diagonal’ β_(zyy) components. The most significant findings of these studies are: (i) β_0 values as much as 6 times that of the chromophore in the technologically important material (E)-4′-(dimethylamino)-N-methyl-4-stilbazolium tosylate; (ii) reversible electrochemistry that offers potential for redox-switching of optical properties over multiple states; (iii) strongly 2D NLO responses that may be exploited for novel practical applications; (iv) a new polar material, suitable for bulk NLO behavior

Can impact excitation explain efficient carrier multiplication in carbon nanotube photodiodes?

We address recent experiments (Science 325, 1367 (2009)) reporting on highly efficient multiplication of electron-hole pairs in carbon nanotube photodiodes at photon energies near the carrier multiplication threshold (twice the quasi-particle band gap). This result is surprising in light of recent experimental and theoretical work on multiexciton generation in other confined materials, such as semiconducting nanocrystals. We propose a detailed mechanism based on carrier dynamics and impact excitation resulting in highly efficient multiplication of electron-hole pairs. We discuss the important time and energy scales of the problem and provide analysis of the role of temperature and the length of the diode

Ab Initio Screening Approach for the Discovery of Lignin Polymer Breaking Pathways

The directed depolymerization of lignin biopolymers is of utmost relevance for the valorization or commercialization of biomass fuels. We present a computational and theoretical screening approach to identify potential cleavage pathways and resulting fragments that are formed during depolymerization of lignin oligomers containing two to six monomers. We have developed a chemical discovery technique to identify the chemically relevant putative fragments in eight known polymeric linkage types of lignin. Obtaining these structures is a crucial precursor to the development of any further kinetic modeling. We have developed this approach by adapting steered molecular dynamics calculations under constant force and varying the points of applied force in the molecule to diversify the screening approach. Key observations include relationships between abundance and breaking frequency, the relative diversity of potential pathways for a given linkage, and the observation that readily cleaved bonds can destabilize adjacent bonds, causing subsequent automatic cleavage.Massachusetts Institute of Technology (Research Support Corporation, Reed Grant)United States. Dept. of Energy. Computational Science Graduate Fellowship Program (DOE-CSGF)Burroughs Wellcome Fund (Career Award at the Scientific Interface

Excited-State Electronic Structure with Configuration Interaction Singles and Tamm–Dancoff Time-Dependent Density Functional Theory on Graphical Processing Units

Excited-state calculations are implemented in a development version of the GPU-based TeraChem software package using the configuration interaction singles (CIS) and adiabatic linear response Tamm–Dancoff time-dependent density functional theory (TDA-TDDFT) methods. The speedup of the CIS and TDDFT methods using GPU-based electron repulsion integrals and density functional quadrature integration allows full ab initio excited-state calculations on molecules of unprecedented size. CIS/6-31G and TD-BLYP/6-31G benchmark timings are presented for a range of systems, including four generations of oligothiophene dendrimers, photoactive yellow protein (PYP), and the PYP chromophore solvated with 900 quantum mechanical water molecules. The effects of double and single precision integration are discussed, and mixed precision GPU integration is shown to give extremely good numerical accuracy for both CIS and TDDFT excitation energies (excitation energies within 0.0005 eV of extended double precision CPU results)