Catalytic role of boron atoms in self-interstitial clustering in Si

Using density functional theory (DFT) calculations and kinetic simulations, we have investigated the influence of boron atoms on self-interstitial clustering in Si. From DFT calculations of neutral interstitial clusters with a single B atom (BsIn, nIn–1 + BsI) becomes substantially weaker than that of an interstitial (BsIn-->BsIn–1 + I) when n>=4. This implies boron can be liberated while leaving an interstitial cluster behind. Our kinetic simulations including the boron liberation explain well experimental observations reported by J. L. Benton et al., J. Appl. Phys. 82, 120 (1997)

Shouldering in B diffusion profiles in Si: Role of di-boron diffusion

The role of di-boron diffusion in evolution of B diffusion profiles has been investigated. We find that boron pair (B-s-B-i) diffusion can become as important as boron-interstitial pair (B-s-Si-i) diffusion when both boron concentration and annealing temperature are very high, leading to concentration-dependent B diffusion. Our simulated B diffusion profiles with dramatic shouldering are in excellent agreement with experimental ones reported by Schroer [Appl. Phys. Lett. 74, 3996 (1999)] for high-temperature (approximate to 1200 degrees C) postimplantion annealing of ultralow-energy (approximate to500 eV) implanted high-concentration (>10(19) cm(-3)) boron in silicon

Computational Study of Copper(II) Complexation and Hydrolysis in Aqueous Solutions Using Mixed Cluster/Continuum Models

We use density functional theory (B3LYP) and the COSMO continuum solvent model to characterize the structure and stability of the hydrated Cu(II) complexes [Cu(MeNH_2)(H_2O)_(n−1)]^(2+) and [Cu(OH)_x(H_2O)_(n−x)]^(2−x) (x = 1−3) as a function of metal coordination number (4−6) and cluster size (n = 4−8, 18). The small clusters with n = 4−8 are found to be the most stable in the nearly square-planar four-coordinate configuration, except for [Cu(OH)_3(H_2O)]^−, which is three-coordinate. In the presence of the two full hydration shells (n = 18), however, the five-coordinate square-pyramidal geometry is the most favorable for Cu(MeNH_2)^(2+) (5, 6) and Cu(OH)^+ (5, 4, 6), and the four-coordinate geometry is the most stable for Cu(OH)_2 (4, 5) and Cu(OH)_3^− (4). (Other possible coordination numbers for these complexes in the aqueous phase are shown in parentheses.) A small energetic difference between these structures (0.23−2.65 kcal/mol) suggests that complexes with different coordination numbers may coexist in solution. Using two full hydration shells around the Cu^(2+) ion (18 ligands) gives Gibbs free energies of aqueous reactions that are in excellent agreement with experiment. The mean unsigned error is 0.7 kcal/mol for the three consecutive hydrolysis steps of Cu^(2+) and the complexation of Cu^(2+) with methylamine. Conversely, calculations for the complexes with only one coordination shell (four equatorial ligands) lead to a mean unsigned error that is >6.0 kcal/mol. Thus, the explicit treatment of the first and the second shells is critical for the accurate prediction of structural and thermodynamic properties of Cu(II) species in aqueous solution

PAMAM Dendrimers Undergo pH Responsive Conformational Changes without Swelling

Atomistic molecular dynamics (MD) simulations of a G4-NH2 PAMAM dendrimer were carried out in aqueous solution using explicit water molecules and counterions (with the Dreiding III force field optimized using quantum mechanics). Our simulations predict that the radius of gyration (R_g) of the dendrimer changes little with pH from 21.1 Å at pH 10 (uncharged PAMAM) to 22.1 Å at pH 5 (charged with 126 protons), which agrees quantitatively with recent small angle neutron scattering (SANS) experiments (from 21.4 Å at pH 10 to 21.5 Å at pH 5). Even so we predict a dramatic change in the conformation. The ion pairing in the low pH form leads to a locally compact dense shell with an internal surface area only 37% of the high pH form with a dense core. This transformation from “dense core” at high pH to “dense shell” at low pH could facilitate the encapsulation and release of guest molecules (e.g., drugs) using pH as the trigger, making dendrimers a unique drug delivery vehicle

Shock Induced Decomposition and Sensitivity of Energetic Materials by ReaxFF Molecular Dynamics

We develop strain-driven compression-expansion technique using molecular dynamics (MD) with reactive force fields (ReaxFF) to study the impact sensitivity of energetic materials. It has been applied to simulation of 1,3,5-trinitrohexahydro-s-triazine (RDX) crystal subjected to high-rate compression typical at the detonation front. The obtained results show that at lower compression ratio x = 1-V/V040%) all molecules decompose very quickly. We have observed both primary and secondary reactions during the decomposition process as well as production of various intermediates (NO2, NO, HONO, OH) and final products (H2O, N2, CO, CO2). The results of strain-driven compression-expansion modeling are in a good agreement with previous ReaxFF-MD shock simulations in RDX. Proposed approach might be useful for a quick test of sensitivity of energetic materials under conditions of high strain rate loading

Predicting Solid-State Heats of Formation of Newly Synthesized Polynitrogen Materials by Using Quantum Mechanical Calculations

We present density functional theory level predictions and analysis of the basic properties of newly synthesized high-nitrogen compounds together with 3,6-bis(2H-tetrazol-5-yl)-1,2,4,5-tetrazine (BTT) and 3,3′-azobis(6-amino-1,2,4,5-tetrazine) (DAAT), for which experimental data are available. The newly synthesized high-nitrogen compounds are based on tricycle fused 1,2,4-triazine and 1,2,4,5-tetrazine heterocycles. In this work, the molecules BTT and DAAT have been studied in order to validate the theoretical approach and to facilitate further progress developments for the molecules of interest. Molecular structural properties are clarified, and IR spectra predictions are provided to help detection of those compounds in the experiment. The energy content of the molecules in the gas phase is evaluated by calculating standard enthalpies of formation, by using a special selection of isodesmic reaction paths. We also include estimates of the condensed-phase heats of formation and heats of sublimation in the framework of the Politzer approach. The obtained properties are consistent with those new high-nitrogen compounds being a promising set of advanced energetic materials

Evaluation of B3LYP, X3LYP, and M06-class density functionals for predicting the binding energies of neutral, protonated, and deprotonated water clusters

In this paper we assess the accuracy of the B3LYP, X3LYP, and newly developed M06-L, M06-2X, and M06 functionals to predict the binding energies of neutral and charged water clusters including (H_2O)_n, n = 2−8, 20), H_3O+(H_2O_)n, n = 1−6, and OH−(H_2O)_n, n = 1−6. We also compare the predicted energies of two ion hydration and neutralization reactions on the basis of the calculated binding energies. In all cases, we use as benchmarks calculated binding energies of water clusters extrapolated to the complete basis set limit of the second-order Møller−Plesset perturbation theory with the effects of higher order correlation estimated at the coupled-cluster theory with single, double, and perturbative triple excitations in the aug-cc-pVDZ basis set. We rank the accuracy of the functionals on the basis of the mean unsigned error (MUE) between calculated benchmark and density functional theory energies. The corresponding MUE (kcal/mol) for each functional is listed in parentheses. We find that M06-L (0.73) and M06 (0.84) give the most accurate binding energies using very extended basis sets such as aug-cc-pV5Z. For more affordable basis sets, the best methods for predicting the binding energies of water clusters are M06-L/aug-cc-pVTZ (1.24), B3LYP/6-311++G(2d,2p) (1.29), and M06/aug-cc-PVTZ (1.33). M06-L/aug-cc-pVTZ also gives more accurate energies for the neutralization reactions (1.38), whereas B3LYP/6-311++G(2d,2p) gives more accurate energies for the ion hydration reactions (1.69)

Foliation of the Kottler-Schwarzschild-De Sitter Spacetime by Flat Spacelike Hypersurfaces

There exist Kruskal like coordinates for the Reissner-Nordstrom (RN) black hole spacetime which are regular at coordinate singularities. Non existence of such coordinates for the extreme RN black hole spacetime has already been shown. Also the Carter coordinates available for the extreme case are not manifestly regular at the coordinate singularity, therefore, a numerical procedure was developed to obtain free fall geodesics and flat foliation for the extreme RN black hole spacetime. The Kottler-Schwarzschild-de Sitter (KSSdS) spacetime geometry is similar to the RN geometry in the sense that, like the RN case, there exist non-singular coordinates when there are two distinct coordinate singularities. There are no manifestly regular coordinates for the extreme KSSdS case. In this paper foliation of all the cases of the KSSdS spacetime by flat spacelike hypersurfaces is obtained by introducing a non-singular time coordinate.Comment: 12 pages, 4 figure

Iridium complexes bearing a PNP ligand, favoring facile C(sp^3)–H bond cleavage

Hydrogen iodide is lost upon reaction of PNP with IrI_3, where PNP = 2,6-bis-(di-t-butylphosphinomethyl)pyridine to give crystallographically characterized Ir(PNP)*(I)_2, which reacts with H_2 to give Ir(PNP)(H)(I)_2. Ir(PNP)(Cl)_3 is relatively inert towards the intramolecular C–H activation of the tert-butyl's of the PNP ligand