Graphene-like two-dimensional ionic boron with double dirac cones at ambient condition

Recently, partially ionic boron (γ-B28) has been predicted and observed in pure boron, in bulk phase and controlled by pressure [Nature, 457 (2009) 863]. By using ab initio evolutionary structure search, we report the prediction of ionic boron at a reduced dimension and ambient pressure, namely, the two-dimensional (2D) ionic boron. This 2D boron structure consists of graphene-like plane and B2 atom pairs, with the P6/mmm space group and 6 atoms in the unit cell, and has lower energy than the previously reported α-sheet structure and its analogues. Its dynamical and thermal stability are confirmed by the phonon-spectrum and ab initio molecular dynamics simulation. In addition, this phase exhibits double Dirac cones with massless Dirac fermions due to the significant charge transfer between the graphene-like plane and B2 pair that enhances the energetic stability of the P6/mmm boron. A Fermi velocity (vf) as high as 2.3 x 106 m/s, which is even higher than that of graphene (0.82 x 106 m/s), is predicted for the P6/mmm boron. The present work is the first report of the 2D ionic boron at atmospheric pressure. The unique electronic structure renders the 2D ionic boron a promising 2D material for applications in nanoelectronics.

Multiscale study on hydrogen mobility in metallic fusion divertor material

For achieving efficient fusion energy production, the plasma-facing wall materials of the fusion reactor should ensure long time operation. In the next step fusion device, ITER, the first wall region facing the highest heat and particle load, i.e. the divertor area, will mainly consist of tiles based on tungsten. During the reactor operation, the tungsten material is slowly but inevitably saturated with tritium. Tritium is the relatively short-lived hydrogen isotope used in the fusion reaction. The amount of tritium retained in the wall materials should be minimized and its recycling back to the plasma must be unrestrained, otherwise it cannot be used for fueling the plasma. A very expensive and thus economically not viable solution is to replace the first walls quite often. A better solution is to heat the walls to temperatures where tritium is released. Unfortunately, the exact mechanisms of hydrogen release in tungsten are not known. In this thesis both experimental and computational methods have been used for studying the release and retention of hydrogen in tungsten. The experimental work consists of hydrogen implantations into pure polycrystalline tungsten, the determination of the hydrogen concentrations using ion beam analyses (IBA) and monitoring the out-diffused hydrogen gas with thermodesorption spectrometry (TDS) as the tungsten samples are heated at elevated temperatures. Combining IBA methods with TDS, the retained amount of hydrogen is obtained as well as the temperatures needed for the hydrogen release. With computational methods the hydrogen-defect interactions and implantation-induced irradiation damage can be examined at the atomic level. The method of multiscale modelling combines the results obtained from computational methodologies applicable at different length and time scales. Electron density functional theory calculations were used for determining the energetics of the elementary processes of hydrogen in tungsten, such as diffusivity and trapping to vacancies and surfaces. Results from the energetics of pure tungsten defects were used in the development of an classical bond-order potential for describing the tungsten defects to be used in molecular dynamics simulations. The developed potential was utilized in determination of the defect clustering and annihilation properties. These results were further employed in binary collision and rate theory calculations to determine the evolution of large defect clusters that trap hydrogen in the course of implantation. The computational results for the defect and trapped hydrogen concentrations were successfully compared with the experimental results. With the aforedescribed multiscale analysis the experimental results within this thesis and found in the literature were explained both quantitatively and qualitatively.

Computer simulation of H and He effects in fusion reactor materials

Fusion power is an appealing source of clean and abundant energy. The radiation resistance of reactor materials is one of the greatest obstacles on the path towards commercial fusion power. These materials are subject to a harsh radiation environment, and cannot fail mechanically or contaminate the fusion plasma. Moreover, for a power plant to be economically viable, the reactor materials must withstand long operation times, with little maintenance. The fusion reactor materials will contain hydrogen and helium, due to deposition from the plasma and nuclear reactions because of energetic neutron irradiation. The first wall divertor materials, carbon and tungsten in existing and planned test reactors, will be subject to intense bombardment of low energy deuterium and helium, which erodes and modifies the surface. All reactor materials, including the structural steel, will suffer irradiation of high energy neutrons, causing displacement cascade damage. Molecular dynamics simulation is a valuable tool for studying irradiation phenomena, such as surface bombardment and the onset of primary damage due to displacement cascades. The governing mechanisms are on the atomic level, and hence not easily studied experimentally. In order to model materials, interatomic potentials are needed to describe the interaction between the atoms. In this thesis, new interatomic potentials were developed for the tungsten-carbon-hydrogen system and for iron-helium and chromium-helium. Thus, the study of previously inaccessible systems was made possible, in particular the effect of H and He on radiation damage. The potentials were based on experimental and ab initio data from the literature, as well as density-functional theory calculations performed in this work. As a model for ferritic steel, iron-chromium with 10% Cr was studied. The difference between Fe and FeCr was shown to be negligible for threshold displacement energies. The properties of small He and He-vacancy clusters in Fe and FeCr were also investigated. The clusters were found to be more mobile and dissociate more rapidly than previously assumed, and the effect of Cr was small. The primary damage formed by displacement cascades was found to be heavily influenced by the presence of He, both in FeCr and W. Many important issues with fusion reactor materials remain poorly understood, and will require a huge effort by the international community. The development of potential models for new materials and the simulations performed in this thesis reveal many interesting features, but also serve as a platform for further studies.

Synthesis and Characterization of a Class of Donor-Acceptor Conjugated Molecules: Experiments and Theoretical Calculations

A class of conjugated molecules containing donor (thiophene) and acceptor (malononitrile) is synthesized by Knoevenagel condensation reaction between 2-(2,6-dimethy1-4H-pyran-4-ylidene) malononitrile and thiophene carbaldehyde containing two and three thiophene units. The resulting molecules are characterized by H-1 and C-13 NMR. We have performed UV-vis absorption, fluorescence, and cyclic voltammetry measurements on these materials. The spectroscopic and electrochemical measurements proved beyond doubt that these materials possess lowexcitation gap and are suitable for being an active material in various electronic devices. We have also performed electronic structure calculations using density functional theory (DFT) and INDO/SCI methods to characterize the ground and excited states of this class of molecules. These donor-acceptor molecules show a strong charge transfercharacter that increases with the increase in the number of thiophene rings coupled to the malononitrile acceptor moiety. We have also calculated the pi-coherence length, Stoke's shift, and effect of solvents on excited states for this class of molecules, Our theoretical values agree well with experimental results.

A Unique Nickel System having Versatile Catalytic Activity of Biological Significance

A new dinuclear nickel(II) complex, [Ni-2(LH2)(H2O)(2)(OH)(NO3)](NO3)(3) (1), of an ``end-off'' compartmental ligand 2,6-bis(N-ethylpiperazine-iminomethyl)-4-methyl-phenolato, has been synthesized and structurally characterized. The X-ray single crystal structure analysis shows that the piperazine moieties assume the expected chair conformation and are protonated. The complex 1 exhibits versatile catalytic activities of biological significance, viz. catecholase, phosphatase, and DNA cleavage activities, etc. The catecholase activity of the complex observed is very dependent on the nature of the solvent. In acetonitrile medium, the complex is inactive to exhibit catecholase activity. On the other hand, in methanol, it catalyzes not only the oxidation of 3,5-di-tert-butylcatechol (3,5-DTBC) but also tetrachlorocatechol (TCC), a catechol which is very difficult to oxidize, under aerobic conditions. UV vis spectroscopic investigation shows that TCC oxidation proceeds through the formation of an intermediate. The intermediate has been characterized by an electron spray ionizaton-mass spectrometry study, which suggests a bidentate rather than a monodentate mode of TCC coordination in that intermediate, and this proposition have been verified by density functional theory calculation. The complex also exhibits phosphatase (with substrate p-nitrophenylphosphate) and DNA cleavage activities. The DNA cleavage activity exhibited by complex 1 most probably proceeds through a hydroxyl radical pathway. The bioactivity study suggests the possible applications of complex 1 as a site specific recognition of DNA and/or as an anticancer agent.

A Unique Nickel System having Versatile Catalytic Activity of Biological Significance

A new dinuclear nickel(II) complex, [Ni-2(LH2)(H2O)(2)(OH)(NO3)](NO3)(3) (1), of an "end-off" compartmental ligand 2,6-bis(N-ethylpiperazine-iminomethyl)-4-methyl-phenolato, has been synthesized and structurally characterized. The X-ray single crystal structure analysis shows that the piperazine moieties assume the expected chair conformation and are protonated. The complex 1 exhibits versatile catalytic activities of biological significance, viz. catecholase, phosphatase, and DNA cleavage activities, etc. The catecholase activity of the complex observed is very dependent on the nature of the solvent. In acetonitrile medium, the complex is inactive to exhibit catecholase activity. On the other hand, in methanol, it catalyzes not only the oxidation of 3,5-di-tert-butylcatechol (3,5-DTBC) but also tetrachlorocatechol (TCC), a catechol which is very difficult to oxidize, under aerobic conditions. UV vis spectroscopic investigation shows that TCC oxidation proceeds through the formation of an intermediate. The intermediate has been characterized by an electron spray ionizaton-mass spectrometry study, which suggests a bidentate rather than a monodentate mode of TCC coordination in that intermediate, and this proposition have been verified by density functional theory calculation. The complex also exhibits phosphatase (with substrate p-nitrophenylphosphate) and DNA cleavage activities. The DNA cleavage activity exhibited by complex 1 most probably proceeds through a hydroxyl radical pathway. The bioactivity study suggests the possible applications of complex 1 as a site specific recognition of DNA and/or as an anticancer agent.

Electronic structure and bonding of beta-rhombohedral boron using cluster fragment approach

Theoretical studies using density functional theory are carried out to understand the electronic structure and bonding and electronic properties of elemental beta-rhombohedral boron. The calculated band structure of ideal beta-rhombohedral boron (B-105) shows valence electron deficiency and depicts metallic behavior. This is in contrast to the experimental result that it is a semiconductor. To understand this ambiguity we discuss the electronic structure and bonding of this allotrope with cluster fragment approach using our recently proposed mno rule. This helps us to comprehend in greater detail the structure of B-105 and materials which are closely related to beta-rhombohedral boron. The molecular structures B12H12-2, B28H21+1, BeB27H21, LiB27H21-1, CB27H21+2, B57H36+3, Be3B54H36, and Li2CB54H36, and corresponding solids Li8Be3B102 and Li10CB102 are arrived at using these ideas and studied using first principles density functional theory calculations.

A Comparative Study of Aromaticity in Substituted Tetracyclic and Hexacyclic Thiophenes

We have studied the nature of aromaticity in expanded porphyrinic analogues of thiophenes formed by four and six thiophenes. Using density functional theory (DFT) we have analyzed the aromaticity of both the macrocycle and individual molecular fragments. We find paramagnetic annulenic ring currents in the case of tetracyclic molecules and diamagnetic annulenic ring currents for hexacyclic molecules. We have also studied the effect of substitution of benzene rings within the macrocycle. We find that as the number of benzene rings is increased the aromaticity increases for tetracyclic systems and decreases for hexacyclic systems. All the results have been analyzed with various microscopic parameters, including geometry, excitation gap, and NMR criteria.

Molecular Dynamics Simulations of Homogeneous Nucleation

Nucleation is the first step in a phase transition where small nuclei of the new phase start appearing in the metastable old phase, such as the appearance of small liquid clusters in a supersaturated vapor. Nucleation is important in various industrial and natural processes, including atmospheric new particle formation: between 20 % to 80 % of atmospheric particle concentration is due to nucleation. These atmospheric aerosol particles have a significant effect both on climate and human health. Different simulation methods are often applied when studying things that are difficult or even impossible to measure, or when trying to distinguish between the merits of various theoretical approaches. Such simulation methods include, among others, molecular dynamics and Monte Carlo simulations. In this work molecular dynamics simulations of the homogeneous nucleation of Lennard-Jones argon have been performed. Homogeneous means that the nucleation does not occur on a pre-existing surface. The simulations include runs where the starting configuration is a supersaturated vapor and the nucleation event is observed during the simulation (direct simulations), as well as simulations of a cluster in equilibrium with a surrounding vapor (indirect simulations). The latter type are a necessity when the conditions prevent the occurrence of a nucleation event in a reasonable timeframe in the direct simulations. The effect of various temperature control schemes on the nucleation rate (the rate of appearance of clusters that are equally able to grow to macroscopic sizes and to evaporate) was studied and found to be relatively small. The method to extract the nucleation rate was also found to be of minor importance. The cluster sizes from direct and indirect simulations were used in conjunction with the nucleation theorem to calculate formation free energies for the clusters in the indirect simulations. The results agreed with density functional theory, but were higher than values from Monte Carlo simulations. The formation energies were also used to calculate surface tension for the clusters. The sizes of the clusters in the direct and indirect simulations were compared, showing that the direct simulation clusters have more atoms between the liquid-like core of the cluster and the surrounding vapor. Finally, the performance of various nucleation theories in predicting simulated nucleation rates was investigated, and the results among other things highlighted once again the inadequacy of the classical nucleation theory that is commonly employed in nucleation studies.

InMnO3: A biferroic

InMnO3 which has a hexagonal structure similar to that of YMnO3 is found to show a canted antiferromagnetic behavior below 50 K (T-N) and a ferroelectric (FE) transition at 500 K accompanied by hystersis. We have determined the structure, polarization, and energetics of the FE and paraelectric (PE) phases of InMnO3 using first-principles density functional theory calculations based on pseudopotentials and a plane-wave basis, and find the polarization of the PE phase to be a half-integer quantum. The difference in polarization of the FE and PE phases calculated along a simple path is different from the absolute value of polarization of the FE phase. A weak piezoelectric response is exhibited by InMnO3 to uniaxial strain.

A series of transition metal-azido extended complexes with various anionic and neutral co-ligands: synthesis, structure and their distinct magnetic behavior

The crystal structures and magnetic properties of five new transition metal-azido complexes with two anionic [pyrazine-2-carboxylate (pyzc) and p-aminobenzoate (paba)] and two neutral [pyrazine (pyz) and pyridine (py)] coligands are reported All five complexes were synthesized bysolvothermal methods The complex [Co-2(pyzc)(2)(N-3)(2)(H2O)(2)](n) (1) is 1D and exhibit canted antiferromagnetism, while the 3D complex [MnNa(pyzc)(N-3)(2)(H2O)(2)](n) (2) has a complicated structure and is weakly ferromagnetic in nature [Mn-2(paba)(2)(N-3)(2)(H2O)(2)](n) (3). is a 2D sheet and the Mn-II ions are found to be antiferromagnetically coupled The isostructural 2D complexes [Cu-3(pyz)(2)(N-3)(6)](n) (4) and [Cu-3(py)(2)(N-3)(6)](n) (5) resemble remarkably in their magnetic properties exhibiting moderately strong ferromagnetism. Density functional theory calculations (B3LYP functional) have been performed to provide a qualitative theoretical interpietation of the overall magnetic behavior shown by these complexes.

Freezing of a colloidal liquid subject to shear flow

A nonequilibrium generalization of the density-functional theory of freezing is proposed to investigate the shear-induced first-order phase transition in colloidal suspensions. It is assumed that the main effect of a steady shear is to break the symmetry of the structure factor of the liquid and that for small shear rate, the phenomenon of a shear-induced order-disorder transition may be viewed as an equilibrium phase transition. The theory predicts that the effective density at which freezing takes place increases with shear rate. The solid (which is assumed to be a bcc lattice) formed upon freezing is distorted and specifically there is less order in one plane compared with the order in the other two perpendicular planes. It is shown that there exists a critical shear rate above which the colloidal liquid does not undergo a transition to an ordered (or partially ordered) state no matter how large the density is. Conversely, above the critical shear rate an initially formed bcc solid always melts into an amorphous or liquidlike state. Several of these predictions are in qualitative agreement with the light-scattering experiments of Ackerson and Clark. The limitations as well as possible extensions of the theory are also discussed.

Completeness of the Coulomb wave functions in quantum mechanics

We give an explicit, direct, and fairly elementary proof that the radial energy eigenfunctions for the hydrogen atom in quantum mechanics, bound and scattering states included, form a complete set. The proof uses only some properties of the confluent hypergeometric functions and the Cauchy residue theorem from analytic function theory; therefore it would form useful supplementary reading for a graduate course on quantum mechanics.

Study of the random vibration of nonlinear systems by the Gaussian closure technique

A technique is developed to study random vibration of nonlinear systems. The method is based on the assumption that the joint probability density function of the response variables and input variables is Gaussian. It is shown that this method is more general than the statistical linearization technique in that it can handle non-Gaussian excitations and amplitude-limited responses. As an example a bilinear hysteretic system under white noise excitation is analyzed. The prediction of various response statistics by this technique is in good agreement with other available results.