Computational Methods for Studies of Quantum Dots and Boson Condensates

This thesis presents ab initio studies of two kinds of physical systems, quantum dots and bosons, using two program packages of which the bosonic one has mainly been developed by the author. The implemented models, \emph{i.e.}, configuration interaction (CI) and coupled cluster (CC) take the correlated motion of the particles into account, and provide a hierarchy of computational schemes, on top of which the exact solution, within the limit of the single-particle basis set, is obtained. The theory underlying the models is presented in some detail, in order to provide insight into the approximations made and the circumstances under which they hold. Some of the computational methods are also highlighted. In the final sections the results are summarized. The CI and CC calculations on multiexciton complexes in self-assembled semiconductor quantum dots are presented and compared, along with radiative and non-radiative transition rates. Full CI calculations on quantum rings and double quantum rings are also presented. In the latter case, experimental and theoretical results from the literature are re-examined and an alternative explanation for the reported photoluminescence spectra is found. The boson program is first applied on a fictitious model system consisting of bosonic electrons in a central Coulomb field for which CI at the singles and doubles level is found to account for almost all of the correlation energy. Finally, the boson program is employed to study Bose-Einstein condensates confined in different anisotropic trap potentials. The effects of the anisotropy on the relative correlation energy is examined, as well as the effect of varying the interaction potential.}

The influence of lone-pair repulsions on C–C bond lengths: a critical evaluation of the experimental and theoretical evidence

X-Ray structural data, as well as semiempirical and ab initio molecular orbital calculations, reveal no systematic and substantial difference between the C–C bond lengths of cis and trans 1,2-diketones. Additional results on various conformations of 1,2-diimines and 1,2-dithiones follow the same pattern. Therefore, lone-pair repulsions cannot be implicated in the observed lengthening of C–C bonds in isatin and several related molecules. Conjugation in these systems occurs peripherally avoiding the participation of the central C–C bond. Negative hyperconjugative interaction between the oxygen lone pairs and the adjacent C–C σ* orbital is suggested to be the principal reason for the relatively long C–C bond in diketones. This effect is found in both the cis and trans conformations.

Patterning Nanoroads and Quantum Dots on Fluorinated Graphene

Using ab initio methods we have investigated the fluorination of graphene and find that different stoichiometric phases can be formed without a nucleation barrier, with the complete “2D-Teflon” CF phase being thermodynamically most stable. The fluorinated graphene is an insulator and turns out to be a perfect matrix-host for patterning nanoroads and quantum dots of pristine graphene. The electronic and magnetic properties of the nanoroads can be tuned by varying the edge orientation and width. The energy gaps between the highest occupied and lowest unoccupied molecular orbitals (HOMO-LUMO) of quantum dots are size-dependent and show a confinement typical of Dirac fermions. Furthermore, we study the effect of different basic coverage of F on graphene (with stoichiometries CF and C4F) on the band gaps, and show the suitability of these materials to host quantum dots of graphene with unique electronic properties.

A parametric method 3 (PM3) study of trans‐stilbene

We report a comparative modified neglect of diatomic overlap (MNDO), Austin method one (AM1), and parametric method 3 (PM3) study of trans‐stilbene (tS) in its ground, excited (singlet and triplet), and ionic (positive and negative polarons and bipolarons) states. We have also calculated the barrier for ring rotation about the backbone single bond. Our results show that PM3 geometries are superior to MNDO and AM1, at least for tS. PM3 predicts, in contrast with MNDO, AM1 and even ab initio 3‐21G, a coplanar structure for tS, in accordance with recent experimental data. Singlet and triplet energies obtained from heats of formation are in surprisingly good agreement with experimental data.

Conformational Studies of Some Alpha-substituted Acetones CH3COCH2XCH3 Where X=CH2, O Or S And Phenacyl Sulfide

Both semiempirical and ab initio calculations are reported for conformational studies of a series of alpha-substituted acetones CH3COCH2XCH3 where X = CH2, O or S and of phenacyl sulfide PhCOCH2SCH3. For conformational studies in the lowest triplet state of these molecules, the MINDO/3 method was employed in the unrestricted Hartree-Fock frame. Results reveal that rotation around a bond alpha to the carbonyl group is more favourable than that around the beta bond. The preferred conformations in the lowest triplet state are nearly the same as in the ground state.

Evidence for O-H.cntdot..cntdot..cntdot.C and N-H.cntdot..cntdot..cntdot.C hydrogen bonding in crystalline alkynes, alkenes, and aromatics

Attempts have been made to understand the nature and significance of hydrogen bonds of the type X-H-C (X = 0, N) . These unusual interactions have been discussed recently. Crystallographic studies on 17a-ethynylandrosta- 2,4-dieno[2,3-d]dihydroxazol1-7 8-01 (donazole) provide direct evidence of such an 0-H.0-C interaction. Ab initio computations, IR spectroscopy, and database studies show that these hydrogen bonds, while uncommon, are energetically and structurally significant.

Is a CImage C quadruple bond possible?

Ab initio computations Predict a minimum for a carbon-carbon quadruple bond in C2. While the computed C-C separation, 1.128 angstrom, fits Pauling's bond order-bond length relationship well, this result is an artifact of the restricted Hartree Fock theory.

Syntheses, spectroscopy, structures, and conformations of .lambda.3-cyclotriphosphazanes: role of negative hyperconjugation

The reactions of As-chlorocyclotriphosphazane [EtNPCl], with phenols or trifluoroethanol yield the respective aryloxy- or trifluoroethoxy-containingX 3-cyclotriphosphazanes [EtNP(OR)]3 (R = C6H4Br-4 (2),C 6H5 (3C,6 H3-Mez-3,5 (4), C6H3Mez-2,6 (5), CH2CF3 (6)) as their cis-transisomericmixtures. The products have beencharacterized by IRand NMRspectroscopy. Thecrystalstructuresofboth thecis (2a) and trans(2b) isomer_softhep-bromophenoxy derivative have been determined by X-ray diffraction. Crystal data for 2a: triclinic, P1, a = 9.872(4) A, b = 13.438(6) A, c = 13.548(8) A, CY = 117.02(5)', 0 = 96.00(6)', y = 105.38(4)O, Z = 2, final R = 0.080. Crystal data for 2b: monoclinic, P21/n, a = 12.721(6) A, b = 13.468(7) A, c = 17.882(5) A, /3 = 101.62(3)O, Z = 4, final R = 0.066. The cis isomer exhibits a chair-triaxial conformation and the trans isomer a boat-triaxial conformation. Conformational preferences of X3-cyclotriphosphazanes have been probed by both MNDO and ab initio calculations on model systems [HNPXIp (X = H, F). In addition to vicinal lone pair repulsions, negative hyperconjugative interactions involving the nitrogen lone pairs and adjacent P-X Q* orbitals are found to be important (especially when X is an electronegative substituent) in determining the conformational preferences of X3-cyclotriphosphazanes. The calculations also show that the axial - equatorial conversion at phosphorus has a large activation barrier in these systems

Binding of active site directed ligands to phospholipase A2: Implications on the molecular constraints and catalytic mechanism

Molecular constraints for the localization of active site directed ligands (competitive inhibitors and substrates) in the active site of phospholipase A2 (PLA2) are characterized. Structure activity relationships with known inhibitors suggest that the head : group interactions dominate the selectivity as well as a substantial part of the affinity. The ab initio fitting of the amide ligands in the active site was carried out to characterize the head group interactions. Based on a systematic coordinate space search, formamide is docked with known experimental constraints such as coordination of the carbonyl group to Ca2+ and hydrogen bond between amide nitrogen and ND1 of His48. An optimal position for a bound water molecule is identified and its significance for the catalytic mechanism is postulated. Unlike the traditional ''pseudo-triad'' mechanism, the ''Ca-coordinatedoxyanion'' mechanism proposed here invokes activation of the catalytic water to form the oxyanion in the coordination sphere of calcium. As it attacks the carbonyl carbon of the ester, a near-tetrahedral intermediate is formed. As the second proton of the catalytic water is abstracted by the ester oxygen, its reorientation and simultaneous cleavage form hydrogen bond with ND1 of His48. In this mechanism of esterolysis, a catalytic role for the water co-ordinated to Ca2+ is recognised.

Factors contributing to phantom bonds in inorganic molecules: Interpretations and predictions

The validity of various qualitative proposals for interpreting and predicting the existence of short contacts between formally non-bonded atoms, as in cyclodisiloxane and related inorganic ring systems, is critically evaluated. The models range from simple considerations of geometric constraints, lone pair repulsions and pi-complex formation to proposals such as the unsupported pi-bond model and the sigma-bridged-pi bond concept. It is pointed out that a unified description based on a combination of closed and open 3-centre 2-electron bonds is possible. The role of hybridisation is emphasized in the short phantom bond computed in an earlier model system. These insights are used to predict structures with exceptionally short Si..Si and B..B phantom bonds. The proposals are confirmed by ab initio calculations.

Reactions of the hexachlorocyclodiphosphazane [MeNPCl3]2 with primary aromatic amines: formation of highly basic bisphosphinimines

Reactions of hexachlorocyclodiphosphazane [MeNPCl3]2 with primary aromatic amines afforded the bisphosphinimine hydrochlorides [(RNH)2(RN)PN(Me)P(NHMe)(NHR)2]+Cl- (R = Ph 1, C6H4Me-4 2 or C6H4OMe-4 3). Dehydrochlorination of 2 and 3 by methanolic KOH yielded highly basic bisphosphinimines [(RNH)2(RN)PN(Me)P(NMe)(NHR)2] (R = C6H4Me-4 4 or C6H4OMe-4 5). Compounds 1-5 have been characterised by elemental analysis and IR and NMR (H-1, C-13, P-31) spectroscopy. The structure of 2 has been confirmed by single-crystal X-ray diffraction. The short P-N bond lengths and the conformations of the PN, units can be explained on the basis of cumulative negative hyperconjugative interactions between nitrogen lone pairs and adjacent P-N sigma* orbitals. Ab initio calculations on the model phosphinimine (H2N)3P=NH and its protonated form suggest that (amino)phosphinimines would be stronger bases compared to many organic bases such as guanidine.

Conformational preferences of α-functionalised keten-S,N-acetals: Potential role of SO and SS interactions in solution

PMR spectra of carbonyl compounds 2a-k reveal significant variations in the population of E and Z isomers on changing the solvent from CDCl3 to DMSO-d6. In non-polar media, the intramolecular N-H…. O hydrogen bonded form is exclusively observed. In DMSO-d6, the alternative Z form is also populated. A similar conformational switch is also noted in the corresponding thiones. Different interpretations are critically analysed. The most consistent explanation is suggested to involvean interplay of N-H….X hydrogen bonding and S…X attractive interaction (X = O,S) in these systems. Ab initio calculations support this interpretation.

Role of free valences and spin densities on the migrating atom in location of saddle points in radical-exchange reactions

The principle of the conservation of bond orders during radical-exchange reactions is examined using Mayer's definition of bond orders. This simple intuitive approximation is not valid in a quantitative sense. Ab initio results reveal that free valences (or spin densities) develop on the migrating atom during reactions. For several examples of hydrogen-transfer reactions, the sum of the reaction coordinate bond orders in the transition state was found to be 0.92 +/- 0.04 instead of the theoretical 1.00 because free valences (or spin densities) develop on the migrating atom during reactions. It is shown that free valence is almost equal to the square of the spin density on the migrating hydrogen atom and the maxima in the free valence (or spin density) profiles coincide (or nearly coincide) with the saddle points in the corresponding energy profiles.

First-Order Hyperpolarizabilities of Sulfophthalein Dyes

In this paper we report the first hyperpolarizabilities (beta) of 12, sulfophthalein dyes. Since these dyes are ionic in nature, their second-order nonlinearities were measured by the hyper-Rayleigh scattering technique in solution. The measured beta values are large and highly solvent dependent. Inclusion of solvent polarity in ab initio estimates of static second-order polarizability does not fully account for the experimental beta values. Contributions from the dissociated forms of the dye in different solvents seem to play an important role in enhancing beta in these systems.

Reactions of cyanogen radical with alkanes and an explanation for negative temperature dependence of rate constants

Reactions of cyanide radicals with alkanes have been investigated by ab initio methods. It is found that the potential energy surface for reaction of CN with a primary C-H bond in methane has a small positive barrier while reactions of CN with a secondary and a tertiary C-H bond in alkanes are barrierless at the correlated level. A simple explanation for the obtained negative temperature dependence of rate constants for reactions of CN with a secondary and a tertiary C-H bond in alkanes are given in terms of the collision theory of bimolecular reactions. It is shown that for barrierless reactions the negative temperature dependence of the rate constants is attributed to the variation of the pre-exponential factor with temperature.