Interatomic potentials for fusion reactor material simulations

Fusion energy is a clean and safe solution for the intricate question of how to produce non-polluting and sustainable energy for the constantly growing population. The fusion process does not result in any harmful waste or green-house gases, since small amounts of helium is the only bi-product that is produced when using the hydrogen isotopes deuterium and tritium as fuel. Moreover, deuterium is abundant in seawater and tritium can be bred from lithium, a common metal in the Earth's crust, rendering the fuel reservoirs practically bottomless. Due to its enormous mass, the Sun has been able to utilize fusion as its main energy source ever since it was born. But here on Earth, we must find other means to achieve the same. Inertial fusion involving powerful lasers and thermonuclear fusion employing extreme temperatures are examples of successful methods. However, these have yet to produce more energy than they consume. In thermonuclear fusion, the fuel is held inside a tokamak, which is a doughnut-shaped chamber with strong magnets wrapped around it. Once the fuel is heated up, it is controlled with the help of these magnets, since the required temperatures (over 100 million degrees C) will separate the electrons from the nuclei, forming a plasma. Once the fusion reactions occur, excess binding energy is released as energetic neutrons, which are absorbed in water in order to produce steam that runs turbines. Keeping the power losses from the plasma low, thus allowing for a high number of reactions, is a challenge. Another challenge is related to the reactor materials, since the confinement of the plasma particles is not perfect, resulting in particle bombardment of the reactor walls and structures. Material erosion and activation as well as plasma contamination are expected. Adding to this, the high energy neutrons will cause radiation damage in the materials, causing, for instance, swelling and embrittlement. In this thesis, the behaviour of a material situated in a fusion reactor was studied using molecular dynamics simulations. Simulations of processes in the next generation fusion reactor ITER include the reactor materials beryllium, carbon and tungsten as well as the plasma hydrogen isotopes. This means that interaction models, {\it i.e. interatomic potentials}, for this complicated quaternary system are needed. The task of finding such potentials is nonetheless nearly at its end, since models for the beryllium-carbon-hydrogen interactions were constructed in this thesis and as a continuation of that work, a beryllium-tungsten model is under development. These potentials are combinable with the earlier tungsten-carbon-hydrogen ones. The potentials were used to explain the chemical sputtering of beryllium due to deuterium plasma exposure. During experiments, a large fraction of the sputtered beryllium atoms were observed to be released as BeD molecules, and the simulations identified the swift chemical sputtering mechanism, previously not believed to be important in metals, as the underlying mechanism. Radiation damage in the reactor structural materials vanadium, iron and iron chromium, as well as in the wall material tungsten and the mixed alloy tungsten carbide, was also studied in this thesis. Interatomic potentials for vanadium, tungsten and iron were modified to be better suited for simulating collision cascades that are formed during particle irradiation, and the potential features affecting the resulting primary damage were identified. Including the often neglected electronic effects in the simulations was also shown to have an impact on the damage. With proper tuning of the electron-phonon interaction strength, experimentally measured quantities related to ion-beam mixing in iron could be reproduced. The damage in tungsten carbide alloys showed elemental asymmetry, as the major part of the damage consisted of carbon defects. On the other hand, modelling the damage in the iron chromium alloy, essentially representing steel, showed that small additions of chromium do not noticeably affect the primary damage in iron. Since a complete assessment of the response of a material in a future full-scale fusion reactor is not achievable using only experimental techniques, molecular dynamics simulations are of vital help. This thesis has not only provided insight into complicated reactor processes and improved current methods, but also offered tools for further simulations. It is therefore an important step towards making fusion energy more than a future goal.

Simplified theory of carrier back-scattering in semiconducting carbon nanotubes: A Kane's model approach

We present a simplified yet analytical formulation of the carrier backscattering coefficient for zig-zag semiconducting single walled carbon nanotubes under diffusive regime. The electron-phonon scattering rate for longitudinal acoustic, optical, and zone-boundary phonon emissions for both inter- and intrasubband transition rates have been derived using Kane's nonparabolic energy subband model.The expressions for the mean free path and diffusive resistance have been formulated incorporating the aforementioned phonon scattering. Appropriate overlap function in Fermi's golden rule has been incorporated for a more general approach. The effect of energy subbands on low and high bias zones for the onset of longitudinal acoustic, optical, and zone-boundary phonon emissions and absorption have been analytically addressed. 90% transmission of the carriers from the source to the drain at 400 K for a 5 mu m long nanotube at 105 V m(-1) has been exhibited. The analytical results are in good agreement with the available experimental data. (c) 2010 American Institute of Physics.

Tunnel transport in polypyrrole at low temperature

The anomalous behaviour of conductivity below 4 K in polypyrrole can be attributed to the possibility of tunnel transport in disordered polaronic systems. The deviation from T-1/3 and T-1/4, depending on disorder, can be due to the onset of tunnel transport between localised states, apart from the hopping contribution to the conductivity. In intermediately and lightly doped polypyrrole films, the tunnel contribution to conductivity increases with decreasing temperature in a narrow temperature range, which is a feature of the presence of polarons taking part in the conduction mechanisms of disordered systems with strong electron-phonon coupling. The transition from hopping to tunneling dominated process can be observed either by the increase in conductivity in some cases or by the saturation of conductivity, depending crucially on the extent of disorder in the sample. In both cases the transition temperature is seen to increase with the reduction in the number of localised states.

Correlations between frequency of infra-red active vibrational modes and copper-oxygen distance in copper oxides. Application to superconductors

The infra-red spectra of a large number of ternary Cu(II) oxides with at least a quasi square-planar coordination of oxygen around the copper ions have been studied. The frequency of the bands with the highest frequency,v max, is found to correlate extremely well with the shortest Cu–O distance.v max increases at an impressive rate of sim20 cm–1 per 0.01 Å when the Cu–O distance becomes less than 1.97 Å, which is the Cu2+–O2– distance in square-planar CuO4 complexes as obtained from empirical ionic radii considerations. The marked sensitivity may be used as a ldquotitrationrdquo procedure not only to assign bands but also to obtain diagnostic information about local coordination in compounds derived, for example, from the YBa2Cu3O7–d structure such as LaCaBaCu3O7–d . The only example where this correlation fails is in the two-layer non-superconducting oxides derived from La2(Ca, Sr)Cu2O6. The significance of this result is discussed. The marked dependence of frequency on the bond-distance is qualitatively examined in terms of an increased electron-phonon coupling to account for the observed tendency of the superconducting transition temperature to go through a maximum as the average basal plane Cu–O distance is decreased.

Large-amplitude chirped coherent phonons in tellurium mediated by ultrafast photoexcited carrier diffusion

We report femtosecond time-resolved reflectivity measurements of coherent phonons in tellurium performed over a wide range of temperatures (3-296 K) and pump-laser intensities. A totally symmetric A(1) coherent phonon at 3.6 THz responsible for the oscillations in the reflectivity data is observed to be strongly positively chirped (i.e., phonon time period decreases at longer pump-probe delay times) with increasing photoexcited carrier density, more so at lower temperatures. We show that the temperature dependence of the coherent phonon frequency is anomalous (i.e, increasing with increasing temperature) at high photoexcited carrier density due to electron-phonon interaction. At the highest photoexcited carrier density of (1.4 x 10(21) cm(-3) and the sample temperature of 3 K, the lattice displacement of the coherent phonon mode is estimated to be as high as similar to 0.24 angstrom. Numerical simulations based on coupled effects of optical absorption and carrier diffusion reveal that the diffusion of carriers dominates the nonoscillatory electronic part of the time-resolved reflectivity. Finally, using the pump-probe experiments at low carrier density of 6 x 10(18) cm(-3), we separate the phonon anharmonicity to obtain the electron-phonon coupling contribution to the phonon frequency and linewidth.

Electronic excitations and alternation of conjugated polymers

The electronic excitations and fluorescence of conjugated polymers are related to large or small alternation ? of the transfer integrals t(1 ± ?) along the backbone. The fluorescence of polysilanes (PSs) and poly (para-phenylenevinylene (PPV) is linked to large ?, which places the one-photon gap Eg below the lowest two-photon gap Ea and reduces distortions due to electron-phonon (e-p) coupling. In contrast to small ? not, vert, similar 0.1 in ?-conjugated polymers, such as polyacetylene (PA), para-conjugated phenyls lead to an extended ?-system with increased alternation, to states localized on each ring and to charge-transfer excitations between them. Surprisingly good agreement is found between semiempirical parametric method 3 (PM3) bond lengths and exact Pariser-Parr-Pople (PPP) ?-bond orders for trans-stilbene, where the PPV bipolarons are confined to two phenyls. Stilbene spectra are consistent with increased alternation and small e-p distortions.

Probing Single and Bilayer Graphene Field Effect Transistors By Raman Spectroscopy

This article is a review of our work related to Raman studies of single layer and bilayer graphenes as a function Fermi level shift achieved by electrochemically top gating a field effect transistor. Combining the transport and in situ Raman studies of the field effect devices, a quantitative understanding is obtained of the phonon renormalization due to doping of graphene. Results are discussed in the light of time dependent perturbation theory, with electron phonon coupling parameter as an input from the density functional theory. It is seen that phonons near and Gamma and K points of the Brillouin zone are renormalized very differently by doping. Further, Gamma-phonon renormalization is different in bilayer graphene as compared to single layer, originating from their different electronic band structures near the zone boundary K-point. Thus Raman spectroscopy is not only a powerful probe to characterize the number of layers and their quality in a graphene sample, but also to quantitatively evaluate electron phonon coupling required to understand the performance of graphene devices.

Degradation and Failure of Field Emitting Carbon Nanotube Arrays

It has been observed experimentally that the collective field emission from an array of Carbon Nanotubes (CNTs) exhibits fluctuation and degradation, and produces thermal spikes, resulting in electro-mechanical fatigue and failure of CNTs. Based on a new coupled multiphysics model incorporating the electron-phonon transport and thermo-electrically activated breakdown, a novel method for estimating accurately the lifetime of CNT arrays has been developed in this paper. The main results are discussed for CNT arrays during the field emission process. It is shown that the time-to-failure of CNT arrays increases with the decrease in the angle of tip orientation. This observation has important ramifications for such areas as biomedical X-ray devices using patterned films of CNTs.

Field emission properties of carbon nanotube arrays with defects and impurities

Field emission from carbon nanotubes (CNTs) in the form of arrays or thin films give rise to several strongly correlated process of electromechanical interaction and degradation. Such processes are mainly due to (1) electron-phonon interaction (2) electromechanical force field leading to stretching of CNTs (3) ballistic transport induced thermal spikes, coupled with high dynamic stress, leading to degradation of emission performance at the device scale. Fairly detailed physics based models of CNTs considering the aspects (1) and (2) above have already been developed by these authors, and numerical results indicate good agreement with experimental results. What is missing in such a system level modeling approach is the incorporation of structural defects and vacancies or charge impurities. This is a practical and important problem due to the fact that degradation of field emission performance is indeed observed in experimental I-V curves. What is not clear from these experiments is whether such degradation in the I-V response is due to dynamic reorientation of the CNTs or due to the defects or due to both of these effects combined. Non-equilibrium Green’s function based simulations using a tight-binding Hamiltonian for single CNT segment show up the localization of carrier density at various locations of the CNTs. About 11% decrease in the drive current with steady difference in the drain current in the range of 0.2-0.4V of the gate voltage was reported in literature when negative charge impurity was introduced at various locations of the CNT over a length of ~20nm. In the context of field emission from CNT tips, a simplistic estimate of defects have been introduced by a correction factor in the Fowler-Nordheim formulae. However, a more detailed physics based treatment is required, while at the same time the device-scale simulation is necessary. The novelty of our present approach is the following. We employ a concept of effective stiffness degradation for segments of CNTs, which is due to structural defects, and subsequently, we incorporate the vacancy defects and charge impurity effects in the Green’s function based approach. Field emission induced current-voltage characteristics of a vertically aligned CNT array on a Cu-Cr substrate is then simulated using a detailed nonlinear mechanistic model of CNTs coupled with quantum hydrodynamics. An array of 10 vertically aligned and each 12 m long CNTs is considered for the device scale analysis. Defect regions are introduced randomly over the CNT length. The result shows the decrease in the longitudinal strain due to defects. Contrary to the expected influence of purely mechanical degradation, this result indicates that the charge impurity and hence weaker transport can lead to a different electromechanical force field, which ultimately can reduce the strain. However, there could be significant fluctuation in such strain field due to electron-phonon coupling. The effect of such fluctuations (with defects) is clearly evident in the field emission current history. The average current also decreases significantly due to such defects.

Modeling the Field Emission Current Fluctuation in Carbon Nanotube Thin Films

Owing to their distinct properties, carbon nanotubes (CNTs) have emerged as promising candidate for field emission devices. It has been found experimentally that the results related to the field emission performance show variability. The design of an efficient field emitting device requires the analysis of the variabilities with a systematic and multiphysics based modeling approach. In this paper, we develop a model of randomly oriented CNTs in a thin film by coupling the field emission phenomena, the electron-phonon transport and the mechanics of single isolated CNT. A computational scheme is developed by which the states of CNTs are updated in time incremental manner. The device current is calculated by using Fowler-Nordheim equation for field emission to study the performance at the device scale.

Effect of dimerization on dynamics of spin-charge separation in Pariser-Parr-Pople model: A time-dependent density matrix renormalization group study

We investigate the effect of static electron-phonon coupling on real-time dynamics of spin and charge transport in pi-conjugated polyene chains. The polyene chain is modeled by the Pariser-Parr-Pople Hamiltonian with dimerized nearest-neighbor parameter t(0)(1 + delta) for short bonds and t(0)(1 - delta) for long bonds, and long-range electron-electron interactions. We follow the time evolution of the spin and charge using time-dependent density matrix renormalization group technique when a hole is injected at one end of the chain in its ground state. We find that spin and charge dynamics followed through spin and charge velocities depend both on chain length and extent of dimerization delta. Analysis of the results requires focusing on physical quantities such as average spin and charge polarizations, particularly in the large dimerization limit. In the dimerization range 0.0 <= delta <= 0.15, spin-charge dynamics is found to have a well-defined behavior, with spin-charge separation (measured as the ratio of charge velocity to spin velocity) as well as the total amount of charge and spin transported in a given time along the chain decreasing as dimerization increases. However, in the range 0.3 <= delta <= 0.5, it is observed that the dynamics of spin and charge transport becomes complicated. It is observed that, for large delta values, spin-charge separation is suppressed and the injected hole fails to travel the entire length of the chain.

Pressure induced phase transition in metallic LaB6: A Raman study

We report high pressure Raman studies on single crystals of metallic LaB6 upto a pressure of 16.$ GPa. Raman spectra shows three lines at 680 cm(-1) (T-2g), 1120 cm(-1) (E-g) and 1258 cm(-1) (A(1g)), associated with the internal modes of B-6 molecule. The T-2g mode shows an asymmetric Fano line shape, arising from the interference between the phonon line and the electronic continuum. The line is fitted with I(omega) = I-0(q + epsilon(2))/(I + epsilon(2)) where epsilon = (omega - omega(0))/Gamma, omega(0) is the phonon frequency renormalised due to electron-phonon self energy corrections, Tis the width parameter proportional to the square of the matrix element of the electron-phonon interaction potential. The parameter a signifies the strength of interference. Most interestingly our pressure data for the T-2g mode shows a significant change in the slope of the mode frequency with pressure d omega(0)/dP and Gamma at 9.5 GPa. This clearly indicates that LaB6 undergoes a subtle phase transition at 9.5 GPa within the metallic phase.

Molecular charge-transfer interaction with single-layer graphene

While the effect of electrochemical doping on single-layer graphene (SG) with holes and electrons has been investigated, the effect of charge-transfer doping on SG has not been examined hitherto. Effects of varying the concentration of electron donor and acceptor molecules such as tetrathiafulvalene (TTF) and tetracyanoethylene (TCNE) on SG produced by mechanical exfoliation as well as by the reduction of single-layer graphene oxide have been investigated. TTF softens the G-band in the Raman spectrum, whereas TCNE stiffens the G-band. The full-width-at-half-maximum of the G-band increases on interaction with both TTF and TCNE. These effects are similar to those found with few-layer graphene, but in contrast to those found with electrochemical doping. A common feature between the two types of doping is found in the case of the 2-D band, which shows softening and stiffening on electron and hole doping, respectively. The experimental results are explained on the basis of the frequency shifts, electron-phonon coupling and structural inhomogeneities that are relevant to molecule-graphene interaction.

Gap-Dependent Quasiparticle Dynamics and Coherent Acoustic Phonons in CaFe2As2 across Spin Density Wave Phase Transition

We report ultrafast quasiparticle (QP) dynamics and coherent acoustic phonons in undoped CaFe2As2 iron pnictide single crystals exhibiting spin-density wave (SDW) and concurrent structural phase transition at temperature T-SDW similar to 165K using femtosecond time-resolved pump-probe spectroscopy. The contributions in transient differential reflectivity arising from exponentially decaying QP relaxation and oscillatory coherent acoustic phonon mode show large variations in the vicinity of T-SDW. From the temperature-dependence of the QP recombination dynamics in the SDW phase, we evaluate a BCS-like temperature dependent charge gap with its zero-temperature value of similar to(1.6 perpendicular to 0.2)k(B)T(SDW), whereas, much above T-SDW, an electron-phonon coupling constant of similar to 0.13 has been estimated from the linear temperature-dependence of the QP relaxation time. The long-wavelength coherent acoustic phonons with typical time-period of similar to 100 ps have been analyzed in the light of propagating strain pulse model providing important results for the optical constants, sounds velocity and the elastic modulus of the crystal in the whole temperature range of 3 to 300 K.

Surface ferromagnetism and superconducting properties of nanocrystalline niobium nitride

Nanocrystalline delta-NbNx samples have been synthesized by reacting NbCl5 and urea at three different temperatures. A comparison of their structural, magnetic, transport and thermal properties is reported in the present study. The size of the particles and their agglomeration extent increase with increasing reaction temperature. The sample prepared at 900 degrees C showed the highest superconducting transition temperature (T-c) of 16.2 K with a transition width, similar to 1.8 K, as obtained from the resistivity measurement on cold-pressed bars. Above T-c, magnetization measurements revealed the presence of surface ferromagnetism which coexists with superconductivity below T-c. Heat capacity measurements confirm superconductivity with strong electron-phonon coupling constant. The sample prepared at 800 degrees C shows a lower T-c (10 K) while that prepared at 700 degrees C exhibit no superconductivity down to the lowest temperature (3 K) measured.