Dendrimer building toolkit: Model building and characterization of various dendrimer architectures

We have developed a graphical user interface based dendrimer builder toolkit (DBT) which can be used to generate the dendrimer configuration of desired generation for various dendrimer architectures. The validation of structures generated by this tool was carried out by studying the structural properties of two well known classes of dendrimers: ethylenediamine cored poly(amidoamine) (PAMAM) dendrimer, diaminobutyl cored poly(propylene imine) (PPI) dendrimer. Using full atomistic molecular dynamics (MD) simulation we have calculated the radius of gyration, shape tensor and monomer density distribution for PAMAM and PPI dendrimer at neutral and high pH. A good agreement between the available simulation and experimental (small angle X-ray and neutron scattering; SAXS, SANS) results and calculated radius of gyration was observed. With this validation we have used DBT to build another new class of nitrogen cored poly(propyl ether imine) dendrimer and study it's structural features using all atomistic MD simulation. DBT is a versatile tool and can be easily used to generate other dendrimer structures with different chemistry and topology. The use of general amber force field to describe the intra-molecular interactions allows us to integrate this tool easily with the widely used molecular dynamics software AMBER. This makes our tool a very useful utility which can help to facilitate the study of dendrimer interaction with nucleic acids, protein and lipid bilayer for various biological applications. (c) 2012 Wiley Periodicals, Inc.

Contrasting Electronic Requirements for C-H Binding and C-H Activation in d(6) Half-Sandwich Complexes of Rhenium and Tungsten

A computational study of the interaction half-sandwich metal fragments (metal=Re/W, electron count=d(6)), containing linear nitrosyl (NO+), carbon monoxide (CO), trifluorophosphine (PF3), N-heterocyclic carbene (NHC) ligands with alkanes are conducted using density functional theory employing the hybrid meta-GGA functional (M06). Electron deficiency on the metal increases with the ligand in the order NHC < CO < PF3 < NO+. Electron-withdrawing ligands like NO+ lead to more stable alkane complexes than NHC, a strong electron donor. Energy decomposition analysis shows that stabilization is due to orbital interaction involving charge transfer from the alkane to the metal. Reactivity and dynamics of the alkane fragment are facilitated by electron donors on the metal. These results match most of the experimental results known for CO and PF3 complexes. The study suggests activation of alkane in metal complexes to be facile with strong donor ligands like NHC. (C) 2015 Wiley Periodicals, Inc.

Continuum in the X-Z-Y Weak Bonds: Z= Main Group Elements

The Continuum in the variation of the X-Z bond length change from blue-shifting to red-shifting through zero-shifting in the X-Z---Y complex is inevitable. This has been analyzed by ab-initio molecular orbital calculations using Z= Hydrogen, Halogens, Chalcogens, and Pnicogens as prototypical examples. Our analysis revealed that, the competition between negative hyperconjugation within the donor (X-Z) molecule and Charge Transfer (CT) from the acceptor (Y) molecule is the primary reason for the X-Z bond length change. Here, we report that, the proper tuning of X-and Y-group for a particular Z-can change the blue-shifting nature of X-Z bond to zero-shifting and further to red-shifting. This observation led to the proposal of a continuum in the variation of the X-Z bond length during the formation of X-Z---Y complex. The varying number of orbitals and electrons available around the Z-atom differentiates various classes of weak interactions and leads to interactions dramatically different from the H-Bond. Our explanations based on the model of anti-bonding orbitals can be transferred from one class of weak interactions to another. We further take the idea of continuum to the nature of chemical bonding in general. (C) 2015 Wiley Periodicals, Inc.

Advancing EDL technologies for future space missions : from ground testing facilities to ablative heatshields

Motivated by recent MSL results where the ablation rate of the PICA heatshield was over-predicted, and staying true to the objectives outlined in the NASA Space Technology Roadmaps and Priorities report, this work focuses on advancing EDL technologies for future space missions.

Due to the difficulties in performing flight tests in the hypervelocity regime, a new ground testing facility called the vertical expansion tunnel is proposed. The adverse effects from secondary diaphragm rupture in an expansion tunnel may be reduced or eliminated by orienting the tunnel vertically, matching the test gas pressure and the accelerator gas pressure, and initially separating the test gas from the accelerator gas by density stratification. If some sacrifice of the reservoir conditions can be made, the VET can be utilized in hypervelocity ground testing, without the problems associated with secondary diaphragm rupture.

The performance of different constraints for the Rate-Controlled Constrained-Equilibrium (RCCE) method is investigated in the context of modeling reacting flows characteristic to ground testing facilities, and re-entry conditions. The effectiveness of different constraints are isolated, and new constraints previously unmentioned in the literature are introduced. Three main benefits from the RCCE method were determined: 1) the reduction in number of equations that need to be solved to model a reacting flow; 2) the reduction in stiffness of the system of equations needed to be solved; and 3) the ability to tabulate chemical properties as a function of a constraint once, prior to running a simulation, along with the ability to use the same table for multiple simulations.

Finally, published physical properties of PICA are compiled, and the composition of the pyrolysis gases that form at high temperatures internal to a heatshield is investigated. A necessary link between the composition of the solid resin, and the composition of the pyrolysis gases created is provided. This link, combined with a detailed investigation into a reacting pyrolysis gas mixture, allows a much needed consistent, and thorough description of many of the physical phenomena occurring in a PICA heatshield, and their implications, to be presented.

Through the use of computational fluid mechanics and computational chemistry methods, significant contributions have been made to advancing ground testing facilities, computational methods for reacting flows, and ablation modeling.

First principles based multiparadigm modeling of electronic structures and dynamics

Electronic structures and dynamics are the key to linking the material composition and structure to functionality and performance.

An essential issue in developing semiconductor devices for photovoltaics is to design materials with optimal band gaps and relative positioning of band levels. Approximate DFT methods have been justified to predict band gaps from KS/GKS eigenvalues, but the accuracy is decisively dependent on the choice of XC functionals. We show here for CuInSe₂ and CuGaSe₂, the parent compounds of the promising CIGS solar cells, conventional LDA and GGA obtain gaps of 0.0-0.01 and 0.02-0.24 eV (versus experimental values of 1.04 and 1.67 eV), while the historically first global hybrid functional, B3PW91, is surprisingly the best, with band gaps of 1.07 and 1.58 eV. Furthermore, we show that for 27 related binary and ternary semiconductors, B3PW91 predicts gaps with a MAD of only 0.09 eV, which is substantially better than all modern hybrid functionals, including B3LYP (MAD of 0.19 eV) and screened hybrid functional HSE06 (MAD of 0.18 eV).

The laboratory performance of CIGS solar cells (> 20% efficiency) makes them promising candidate photovoltaic devices. However, there remains little understanding of how defects at the CIGS/CdS interface affect the band offsets and interfacial energies, and hence the performance of manufactured devices. To determine these relationships, we use the B3PW91 hybrid functional of DFT with the AEP method that we validate to provide very accurate descriptions of both band gaps and band offsets. This confirms the weak dependence of band offsets on surface orientation observed experimentally. We predict that the CBO of perfect CuInSe₂/CdS interface is large, 0.79 eV, which would dramatically degrade performance. Moreover we show that band gap widening induced by Ga adjusts only the VBO, and we find that Cd impurities do not significantly affect the CBO. Thus we show that Cu vacancies at the interface play the key role in enabling the tunability of CBO. We predict that Na further improves the CBO through electrostatically elevating the valence levels to decrease the CBO, explaining the observed essential role of Na for high performance. Moreover we find that K leads to a dramatic decrease in the CBO to 0.05 eV, much better than Na. We suggest that the efficiency of CIGS devices might be improved substantially by tuning the ratio of Na to K, with the improved phase stability of Na balancing phase instability from K. All these defects reduce interfacial stability slightly, but not significantly.

A number of exotic structures have been formed through high pressure chemistry, but applications have been hindered by difficulties in recovering the high pressure phase to ambient conditions (i.e., one atmosphere and room temperature). Here we use dispersion-corrected DFT (PBE-ulg flavor) to predict that above 60 GPa the most stable form of N₂O (the laughing gas in its molecular form) is a 1D polymer with an all-nitrogen backbone analogous to cis-polyacetylene in which alternate N are bonded (ionic covalent) to O. The analogous trans-polymer is only 0.03-0.10 eV/molecular unit less stable. Upon relaxation to ambient conditions both polymers relax below 14 GPa to the same stable non-planar trans-polymer, accompanied by possible electronic structure transitions. The predicted phonon spectrum and dissociation kinetics validate the stability of this trans-poly-NNO at ambient conditions, which has potential applications as a new type of conducting polymer with all-nitrogen chains and as a high-energy oxidizer for rocket propulsion. This work illustrates in silico materials discovery particularly in the realm of extreme conditions.

Modeling non-adiabatic electron dynamics has been a long-standing challenge for computational chemistry and materials science, and the eFF method presents a cost-efficient alternative. However, due to the deficiency of FSG representation, eFF is limited to low-Z elements with electrons of predominant s-character. To overcome this, we introduce a formal set of ECP extensions that enable accurate description of p-block elements. The extensions consist of a model representing the core electrons with the nucleus as a single pseudo particle represented by FSG, interacting with valence electrons through ECPs. We demonstrate and validate the ECP extensions for complex bonding structures, geometries, and energetics of systems with p-block character (C, O, Al, Si) and apply them to study materials under extreme mechanical loading conditions.

Despite its success, the eFF framework has some limitations, originated from both the design of Pauli potentials and the FSG representation. To overcome these, we develop a new framework of two-level hierarchy that is a more rigorous and accurate successor to the eFF method. The fundamental level, GHA-QM, is based on a new set of Pauli potentials that renders exact QM level of accuracy for any FSG represented electron systems. To achieve this, we start with using exactly derived energy expressions for the same spin electron pair, and fitting a simple functional form, inspired by DFT, against open singlet electron pair curves (H₂ systems). Symmetric and asymmetric scaling factors are then introduced at this level to recover the QM total energies of multiple electron pair systems from the sum of local interactions. To complement the imperfect FSG representation, the AMPERE extension is implemented, and aims at embedding the interactions associated with both the cusp condition and explicit nodal structures. The whole GHA-QM+AMPERE framework is tested on H element, and the preliminary results are promising.

Effects of phase-breaking on long-range charge transfer in DNA: Partially-coherent-tunneling model study

The mechanism of hole charge transfer in DNA of various lengths and sequences is investigated based on a partially coherent tunneling theory (Zhang et al., J Chem Phys 117:4578, 2002), where the effects of phase-breaking in adenine-thymine and guanine-cytosine base pairs are treated on equal foot. This work aims at providing a self-consistent microscopic interpretation for rate experiments on various DNA systems. We will also clarify the condition under which the simple superexchange-mediated-hopping picture is valid, and make some comments on the further development of present theory.

Density functional study on the reaction mechanism of palladium-catalyzed addition of cyanoboranes to Alkynes

The B3LYP hybrid density functional method has been carried Out to Study theoretically the mechanisin of Pd(0)-catalyzed alkyne cyanoboration reaction. Both the intermolecular and intramolecular alkyne cyanoboration reactions were studied. For each reaction, three paths were proposed. In path A of each reaction, the first step is B-CN bond oxidative addition to bisphosphine complex Pd(PH3)(2), in path B of each reaction, the first step is alkyne coordination to bisphosphine complex Pd(PH3)2, and in path C of each reaction, the first step is the PH3 dissociation front Pd(PH3)2 to form monophosphine complex Pd(PH3) For both reactions, path B is favored.

Investigation of Chemical Bond Characteristics, Thermal Expansion Coefficients and Bulk Moduli of alpha-R2MoO6 and R2Mo2O7 (R = rare earths) by Using a Dielectric Chemical Bond Method

Theoretical researches are performed on the alpha-R2MoO6 (R = Y, Gd, Tb Dy, Ho, Er, Tm and Yb) and pyrochlore-type R2Mo2O7 (R = Y, Nd, Sm, Gd, Tb and Dy) rare earth molybdates by using chemical bond theory of dielectric description. The chemical bonding characteristics and their relationship with thermal expansion property and compressibility are explored. The calculated values of linear thermal expansion coefficient (LTEC) and bulk modulus agree well with the available experimental values. The calculations reveal that the LTECs and the bulk moduli do have linear relationship with the ionic radii of the lanthanides: the LTEC decreases from 6.80 to 6.62 10(-6)/K and the bulk modulus increases from 141 to 154 GPa when R goes in the order Gd, Tb Dy, Ho, Er, Tm, and Yb in the alpha-R2MoO6 series; while in the R2Mo2O7 series, the LTEC ranges from 6.80 to 6.61 10(-6)/K and the bulk modulus ranges from 147 to 163 GPa when R varies in the order Nd, Sm, Gd, Tb and Dy.

Reaction Mechanism of Palladium-Catalyzed Silastannation of Allenes by Density Functional Theory

The reaction mechanism of Pd(O)-catalyzed allenes silastannation reaction is investigated by the density functional method B3LYP. The overall reaction mechanism is examined. For the allene insertion step, the Pd-Si bond is preferred over the Pd-Sn bond. The electronic mechanism of the allene insertion into Pd-Si bond to form sigma-vinylpalladium (terminal-insertion) and sigma-allylpalladium (internal-insertion) insertion products is discussed ill terms of the electron donation and back-donation. It is found that the electron back-donation is significant for both terminal- and internal-insertion. During allene insertion into Pd-Si bond, internal-insertion is preferred over terminal-insertion. By using methylallene, the regio-selectivity for the monosubstituted allene insertion into Pd-Si and Pd-Sn bond is analyzed.

First Principles Investigation on the Ultra-Incompressible and Hard TaN

The structural, electronic, and mechanical properties of TaN were investigated by use of the density functional theory (DFT). Eight structures were considered, i.e.. hexagonal WC TaN, NiAs, wurtzite, and CoSn structures. cubic NaCl. zinc-blende and CsCl structures. The results indicate that TaN in TaN-type structure is the most stable at ambient conditions among the considered structures. Above 5 GPa, TaN in WC-type structure becomes energetically the most stable phase. They are also stable both thermodynamically and mechanically. TaN in WC-type has the largest shear Modulus 243 GPa and large bulk modulus 337 GPa among the considered structures. The Volume compressibility is slightly larger than diamond, but smaller than c-BN at pressures from 0 to 100 GPa. The compressibility along the c axis is smaller than the linear compressibility of both diamond and c-BN.

SOLVENT AND SUBSTITUENT EFFECTS ON THE INTRAMOLECULAR AMIDE HYDROLYSIS OF N-METHYLMALEAMIC ACID

Intramolecular amide hydrolysis of N-methylmaleamic acid is revisited at the B3LYP/6-311G(2df,p)//B3LYP/6-31G(d,p)+ZVPE level, including solvent effects at the CPCM-B3LYP/6-311G(2df,p)//Onsager-B3LYP/6-31G(d,p)+ZPVE level. The concerted reaction mechanism is energetically favorable over stepwise reaction mechanisms in both the gas phase and solution. The calculated reaction barriers are significantly lower in solution than in the gas phase. In addition, it is concluded that the substituents of the four N-methylmaleamic acid derivatives considered herein have a significant effect on the gas-phase reaction barriers but a smaller, or little, effect on the barriers in solution.

Electronic structures and chemical bonding in diatomic ScX to ZnX (X = S, Se, Te)

Bond distances, vibrational frequencies, electron affinities, ionization potentials, dissociation energies, and dipole moments of the title molecules in neutral, positively, and negatively charged ions were studied using density functional method. Ground electronic state was assigned for each molecule. The bonding patterns were analyzed and compared with both the available data and across the series. It was found that, besides ionic component, covalent bonds are formed between the metal s, d orbitals, and the p orbital of S, Se, and Te. For neutral and cationic molecules, the covalent character increases from ScX to CrX and from FeX to CuX with an exception of decrease at MnX and ZnX, while for anionic molecules, the trend is not obvious. For both neutral and charged molecules, the sulfides have the shortest bond distance and largest vibrational frequency, while tellurides have the largest bond distance and smallest vibrational frequency. For neutral and anionic molecules, the dissociation energy of sulfides is the largest, that of tellurides is the smallest, while this only remains true for cationic molecules from ScX+ to FeX+.

Electronic structures and chemical bonding in 4d transition metal monohalides

Bond distances, vibrational frequencies, dipole moments, dissociation energies, electron affinities, and ionization potentials of NIX (XM = Y-Cd, X = F, Cl, Br, I) molecules in neutral, positively, and negatively charged ions were studied by density functional method, B3LYP. The bonding patterns were analyzed and compared with both the available data and across the series. It was found that besides ionic component, covalent bonds are formed between the 4d transition metal s, d orbitals, and the p orbital of halogen. For both neutral and charged molecules, the fluorides have the shortest bond distance, iodides the longest. Although the opposite situation is observed for vibrational frequency, that is, fluorides have the largest value, iodides the smallest. For neutral and anionic species, the dissociation energy tends to decrease with the increasing atomic number from Y to Cd, suggesting the decreasing or weakening of the bond strength. For cationic species, the trend is observed from Y to Ag.