Advanced density functional theory methods for materials science

In this work we chiefly deal with two broad classes of problems in computational materials science, determining the doping mechanism in a semiconductor and developing an extreme condition equation of state. While solving certain aspects of these questions is well-trodden ground, both require extending the reach of existing methods to fully answer them. Here we choose to build upon the framework of density functional theory (DFT) which provides an efficient means to investigate a system from a quantum mechanics description.

Zinc Phosphide (Zn₃P₂) could be the basis for cheap and highly efficient solar cells. Its use in this regard is limited by the difficulty in n-type doping the material. In an effort to understand the mechanism behind this, the energetics and electronic structure of intrinsic point defects in zinc phosphide are studied using generalized Kohn-Sham theory and utilizing the Heyd, Scuseria, and Ernzerhof (HSE) hybrid functional for exchange and correlation. Novel 'perturbation extrapolation' is utilized to extend the use of the computationally expensive HSE functional to this large-scale defect system. According to calculations, the formation energy of charged phosphorus interstitial defects are very low in n-type Zn₃P₂ and act as 'electron sinks', nullifying the desired doping and lowering the fermi-level back towards the p-type regime. Going forward, this insight provides clues to fabricating useful zinc phosphide based devices. In addition, the methodology developed for this work can be applied to further doping studies in other systems.

Accurate determination of high pressure and temperature equations of state is fundamental in a variety of fields. However, it is often very difficult to cover a wide range of temperatures and pressures in an laboratory setting. Here we develop methods to determine a multi-phase equation of state for Ta through computation. The typical means of investigating thermodynamic properties is via ’classical’ molecular dynamics where the atomic motion is calculated from Newtonian mechanics with the electronic effects abstracted away into an interatomic potential function. For our purposes, a ’first principles’ approach such as DFT is useful as a classical potential is typically valid for only a portion of the phase diagram (i.e. whatever part it has been fit to). Furthermore, for extremes of temperature and pressure quantum effects become critical to accurately capture an equation of state and are very hard to capture in even complex model potentials. This requires extending the inherently zero temperature DFT to predict the finite temperature response of the system. Statistical modelling and thermodynamic integration is used to extend our results over all phases, as well as phase-coexistence regions which are at the limits of typical DFT validity. We deliver the most comprehensive and accurate equation of state that has been done for Ta. This work also lends insights that can be applied to further equation of state work in many other materials.

Theory of s-d exchange scattering in dilute magnetic alloys

The problem of s-d exchange scattering of conduction electrons off localized magnetic moments in dilute magnetic alloys is considered employing formal methods of quantum field theoretical scattering. It is shown that such a treatment not only allows for the first time, the inclusion of multiparticle intermediate states in single particle scattering equations but also results in extremely simple and straight forward mathematical analysis. These equations are proved to be exact in the thermodynamic limit. A self-consistent integral equation for electron self energy is derived and approximately solved. The ground state and physical parameters of dilute magnetic alloys are discussed in terms of the theoretical results. Within the approximation of single particle intermediate states our results reduce to earlier versions. The following additional features are found as a consequence of the inclusion of multiparticle intermediate states;

(i) A non analytic binding energy is pre sent for both, antiferromagnetic (J < o) and ferromagnetic (J > o) couplings of the electron plus impurity system.

(ii) The correct behavior of the energy difference of the conduction electron plus impurity system and the free electron system is found which is free of unphysical singularities present in earlier versions of the theories.

(iii) The ground state of the conduction electron plus impurity system is shown to be a many-body condensate state for J < o and J > o, both. However, a distinction is made between the usual terminology of "Singlet" and "Triplet" ground states and nature of our ground state.

(iv) It is shown that a long range ordering, leading to an ordering of the magnetic moments can result from a contact interaction such as the s-d exchange interaction.

(v) The explicit dependence of the excess specific heat of the Kondo systems is obtained and found to be linear in temperatures as T→ o and T ℓnT for 0.3 T_K ≤ T ≤ 0.6 T_K. A rise in (ΔC/T) for temperatures in the region 0 < T ≤ 0.1 T_K is predicted. These results are found to be in excellent agreement with experiments.

(vi) The existence of a critical temperature for Ferromagnetic coupling (J > o) is shown. On the basis of this the apparent contradiction of the simultaneous existence of giant moments and Kondo effect is resolved.

Radio frequency studies of surface resistance and critical magnetic field of type I and type II superconductors

The surface resistance and the critical magnetic field of lead electroplated on copper were studied at 205 MHz in a half-wave coaxial resonator. The observed surface resistance at a low field level below 4.2°K could be well described by the BCS surface resistance with the addition of a temperature independent residual resistance. The available experimental data suggest that the major fraction of the residual resistance in the present experiment was due to the presence of an oxide layer on the surface. At higher magnetic field levels the surface resistance was found to be enhanced due to surface imperfections.

The attainable rf critical magnetic field between 2.2°K and T_c of lead was found to be limited not by the thermodynamic critical field but rather by the superheating field predicted by the one-dimensional Ginzburg-Landau theory. The observed rf critical field was very close to the expected superheating field, particularly in the higher reduced temperature range, but showed somewhat stronger temperature dependence than the expected superheating field in the lower reduced temperature range.

The rf critical magnetic field was also studied at 90 MHz for pure tin and indium, and for a series of SnIn and InBi alloys spanning both type I and type II superconductivity. The samples were spherical with typical diameters of 1-2 mm and a helical resonator was used to generate the rf magnetic field in the measurement. The results of pure samples of tin and indium showed that a vortex-like nucleation of the normal phase was responsible for the superconducting-to-normal phase transition in the rf field at temperatures up to about 0.98-0.99 T_c' where the ideal superheating limit was being reached. The results of the alloy samples showed that the attainable rf critical fields near T_c were well described by the superheating field predicted by the one-dimensional GL theory in both the type I and type II regimes. The measurement was also made at 300 MHz resulting in no significant change in the rf critical field. Thus it was inferred that the nucleation time of the normal phase, once the critical field was reached, was small compared with the rf period in this frequency range.

Microwave dynamics of quasiparticles and critical fields in superconducting films

The microwave response of the superconducting state in equilibrium and non-equilibrium configurations was examined experimentally and analytically. Thin film superconductors were mostly studied in order to explore spatial effects. The response parameter measured was the surface impedance.

For small microwave intensity the surface impedance at 10 GHz was measured for a variety of samples (mostly Sn) over a wide range of sample thickness and temperature. A detailed analysis based on the BCS theory was developed for calculating the surface impedance for general thickness and other experimental parameters. Experiment and theory agreed with each other to within the experimental accuracy. Thus it was established that the samples, thin films as well as bulk, were well characterised at low microwave powers (near equilibrium).

Thin films were perturbed by a small dc supercurrent and the effect on the superconducting order parameter and the quasiparticle response determined by measuring changes in the surface resistance (still at low microwave intensity and independent of it) due to the induced current. The use of fully superconducting resonators enabled the measurement of very small changes in the surface resistance (< 10^-9 Ω/sq.). These experiments yield information regarding the dynamics of the order parameter and quasiparticle systems. For all the films studied the results could be described at temperatures near T_c by the thermodynamic depression of the order parameter due to the static current leading to a quadratic increase of the surface resistance with current.

For the thinnest films the low temperature results were surprising in that the surface resistance decreased with increasing current. An explanation is proposed according to which this decrease occurs due to an additional high frequency quasiparticle current caused by the combined presence of both static and high frequency fields. For frequencies larger than the inverse of the quasiparticle relaxation time this additional current is out of phase (by π) with the microwave electric field and is observed as a decrease of surface resistance. Calculations agree quantitatively with experimental results. This is the first observation and explanation of this non-equilibrium quasiparticle effect.

For thicker films of Sn, the low temperature surface resistance was found to increase with applied static current. It is proposed that due to the spatial non-uniformity of the induced current distribution across the thicker films, the above purely temporal analysis of the local quasiparticle response needs to be generalised to include space and time non-equilibrium effects.

The nonlinear interaction of microwaves arid superconducting films was also examined in a third set of experiments. The surface impedance of thin films was measured as a function of the incident microwave magnetic field. The experiments exploit the ability to measure the absorbed microwave power and applied microwave magnetic field absolutely. It was found that the applied surface microwave field could not be raised above a certain threshold level at which the absorption increased abruptly. This critical field level represents a dynamic critical field and was found to be associated with the penetration of the app1ied field into the film at values well below the thermodynamic critical field for the configuration of a field applied to one side of the film. The penetration occurs despite the thermal stability of the film which was unequivocally demonstrated by experiment. A new mechanism for such penetration via the formation of a vortex-antivortex pair is proposed. The experimental results for the thinnest films agreed with the calculated values of this pair generation field. The observations of increased transmission at the critical field level and suppression of the process by a metallic ground plane further support the proposed model.

High frequency wave propagation in the earth: theory and observation

The wave-theoretical analysis of acoustic and elastic waves refracted by a spherical boundary across which both velocity and density increase abruptly and thence either increase or decrease continuously with depth is formulated in terms of the general problem of waves generated at a steady point source and scattered by a radially heterogeneous spherical body. A displacement potential representation is used for the elastic problem that results in high frequency decoupling of P-SV motion in a spherically symmetric, radially heterogeneous medium. Through the application of an earth-flattening transformation on the radial solution and the Watson transform on the sum over eigenfunctions, the solution to the spherical problem for high frequencies is expressed as a Weyl integral for the corresponding half-space problem in which the effect of boundary curvature maps into an effective positive velocity gradient. The results of both analytical and numerical evaluation of this integral can be summarized as follows for body waves in the crust and upper mantle:

1) In the special case of a critical velocity gradient (a gradient equal and opposite to the effective curvature gradient), the critically refracted wave reduces to the classical head wave for flat, homogeneous layers.

2) For gradients more negative than critical, the amplitude of the critically refracted wave decays more rapidly with distance than the classical head wave.

3) For positive, null, and gradients less negative than critical, the amplitude of the critically refracted wave decays less rapidly with distance than the classical head wave, and at sufficiently large distances, the refracted wave can be adequately described in terms of ray-theoretical diving waves. At intermediate distances from the critical point, the spectral amplitude of the refracted wave is scalloped due to multiple diving wave interference.

These theoretical results applied to published amplitude data for P-waves refracted by the major crustal and upper mantle horizons (the Pg, P*, and Pn travel-time branches) suggest that the 'granitic' upper crust, the 'basaltic' lower crust, and the mantle lid all have negative or near-critical velocity gradients in the tectonically active western United States. On the other hand, the corresponding horizons in the stable eastern United States appear to have null or slightly positive velocity gradients. The distribution of negative and positive velocity gradients correlates closely with high heat flow in tectonic regions and normal heat flow in stable regions. The velocity gradients inferred from the amplitude data are generally consistent with those inferred from ultrasonic measurements of the effects of temperature and pressure on crustal and mantle rocks and probable geothermal gradients. A notable exception is the strong positive velocity gradient in the mantle lid beneath the eastern United States (2 x 10^-3 sec^-1), which appears to require a compositional gradient to counter the effect of even a small geothermal gradient.

New seismic-refraction data were recorded along a 800 km profile extending due south from the Canadian border across the Columbia Plateau into eastern Oregon. The source for the seismic waves was a series of 20 high-energy chemical explosions detonated by the Canadian government in Greenbush Lake, British Columbia. The first arrivals recorded along this profile are on the Pn travel-time branch. In northern Washington and central Oregon their travel time is described by T = Δ/8.0 + 7.7 sec, but in the Columbia Plateau the Pn arrivals are as much as 0.9 sec early with respect to this line. An interpretation of these Pn arrivals together with later crustal arrivals suggest that the crust under the Columbia Plateau is thinner by about 10 km and has a higher average P-wave velocity than the 35-km-thick, 62-km/sec crust under the granitic-metamorphic terrain of northern Washington. A tentative interpretation of later arrivals recorded beyond 500 km from the shots suggests that a thin 8.4-km/sec horizon may be present in the upper mantle beneath the Columbia Plateau and that this horizon may form the lid to a pronounced low-velocity zone extending to a depth of about 140 km.

Study on the critical submergence of surface vortices and the design of anti-vortex intakes

Surface vortex behavior in front of the tunnel intake was investigated in this paper. The critical submergence of vortex was discussed based on the concept of 'critical spherical sink surface' (CSSS). The vortex formation and evolution at the tunnel intake were analyzed based on the theory of CSSS considering the effect of circulation. A theory was proposed to explain the surface vortex. The theoretical development was verified by the physical model experiments of Xiluodu hydropower station. The radial velocity and vortex circulation were considered as the main factors that influence the formation and evolution of surface vortex. Finally, an anti-vortex intake configuration was proposed to weaken the air-core vortex in front of the tunnel intakes of the hydraulic structures. © 2011 Science China Press and Springer-Verlag Berlin Heidelberg.

Application of level II reliability theory

Level II reliability theory provides an approximate method whereby the reliability of a complex engineering structure which has multiple strength and loading variables may be estimated. This technique has been applied previously to both civil and offshore structures with considerable success. The aim of the present work is to assess the applicability of the method for aircraft structures, and to this end landing gear design is considered in detail. It is found that the technique yields useful information regarding the structural reliability, and further it enables the critical design parameters to be identified.

Know thy neighbour: A normative theory of synaptic depression

Synapses exhibit an extraordinary degree of short-term malleability, with release probabilities and effective synaptic strengths changing markedly over multiple timescales. From the perspective of a fixed computational operation in a network, this seems like a most unacceptable degree of added variability. We suggest an alternative theory according to which short-term synaptic plasticity plays a normatively-justifiable role. This theory starts from the commonplace observation that the spiking of a neuron is an incomplete, digital, report of the analog quantity that contains all the critical information, namely its membrane potential. We suggest that a synapse solves the inverse problem of estimating the pre-synaptic membrane potential from the spikes it receives, acting as a recursive filter. We show that the dynamics of short-term synaptic depression closely resemble those required for optimal filtering, and that they indeed support high quality estimation. Under this account, the local postsynaptic potential and the level of synaptic resources track the (scaled) mean and variance of the estimated presynaptic membrane potential. We make experimentally testable predictions for how the statistics of subthreshold membrane potential fluctuations and the form of spiking non-linearity should be related to the properties of short-term plasticity in any particular cell type.

EFFECTIVE-MASS THEORY FOR SUPERLATTICES GROWN ON (11N)-ORIENTED SUBSTRATES

An effective-mass formulation for superlattices grown on (11N)-oriented substrates is given. It is found that, for GaAs/AlxGa1-xAs superlattices, the hole subband structure and related properties are sensitive to the orientation because of the large anisotropy of the valence band. The energy-level positions for the heavy hole and the optical transition matrix elements for the light hole apparently change with orientation. The heavy- and light-hole energy levels at k parallel-to = 0 can be calculated separately by taking the classical effective mass in the growth direction. Under a uniaxial stress along the growth direction, the energy levels of the heavy and light holes shift down and up, respectively; at a critical stress, the first heavy- and light-hole energy levels cross over. The energy shifts caused by the uniaxial stress are largest for the (111) case and smallest for the (001) case. The optical transition matrix elements change substantially after the crossover of the first heavy- and light-hole energy has occurred.

Pressure effects on the thermodynamics of trans-decahydronaphthalene/polystyrene polymer solutions: Application of the Sanchez-Lacombe lattice fluid theory

The cloud-point temperatures (T-cl's) of trans-decahydronaphthalene(TD)/polystyrene (PS, (M) over bar (w) = 270 000) solutions were determined by light scattering measurements over a range of temperatures (1-16degreesC), pressures (100-900 bar), and compositions (4.2-21.6 vol.-% polymer). The system phase separates upon cooling and T-cl was found to increase with rising pressure for constant composition. In the absence of special effects, this finding indicates positive excess volume for the mixing. Special attention was paid to the demixing temperatures as a function of pressure for different polymer solutions and the plots in the T-phi plane (where phi signifies volume fractions). The cloud-point curves of polymer solutions under different pressures were observed for different compositions, which demonstrated that pressure has a greater effect on the TD/PS solutions when far from the critical point as opposed to near the critical point. The Sanchez-Lacombe lattice fluid theory (SLLFT) was used to calculate the spinodals, the binodals, the Flory-Huggins (FH) interaction parameter, the enthalpy of mixing, and the volume changes of mixing. The calculated results show that modified PS scaling parameters can describe the thermodynamics of the TD/PS system well. Moreover the SLLFT describes the experimental results well.

Thermodynamics of PMMA/SAN blends: Application of the Sanchez-Lacombe lattice fluid theory

Polymer blends of poly(methyl methacrylate) (PMMA) and poly(styrene-co-acrylonitrile) (SAN) with an acrylonitrile content of about 30 wt % were prepared by means of solution-casting and characterized by virtue of pressure-volume-temperature (PVT) dilatometry. The Sanchez-Lacombe (SL) lattice fluid theory was used to calculate the spinodals, the binodals, the Flory-Huggins (FH) interaction parameter, the enthalpy of the mixing, the volume change of the mixing, and the combinatorial and vacancy entropies of the mixing for the PMMA/SAN system. A new volume-combining rule was used to evaluate the close-packed volume per mer, upsilon*, of the PMMA/SAN blends. The calculated results showed that the new and the original volume-combining rules had a slight influence on the FH interaction parameter, the enthalpy of the mixing, and the combinatorial entropy of the mixing. Moreover, the spinodals and the binodals calculated with the SL theory by means of the new volume-combining rule could coincide with the measured data for the PMMA/SAN system with a lower critical solution temperature, whereas those obtained by means of the original one could not.

Application of the Sanchez-Lacombe lattice fluid theory to the system pvme/ps and model calculations

Cloud point curves reported in the literature for five representatives of the system poly(vinyl methyl ether)/polystyrene were evaluated theoretically by means of the Sanchez-Lacombe lattice fluid theory. The measured phase separation behavior can be described within experimental error using only one adjustable parameter (quantifying the interaction between the unlike mers). The Flory-Huggins interaction parameters calculated from this theoretical description depend in good approximation linearly on composition (volume fractions) and on the inverse temperature. An evaluation of these data yields a maximum heat effect which is almost one order of magnitude less (ca. -0.25 J/cm(3)) than obtained via Hess's cycle (dissolution of the components and of the blend) from calorimetric measurements. Model calculations on the basis of the present theory demonstrate that the critical points shift to a different extent upon a certain relative change in the molar mass of the blend components. The sensitivity of the calculated phase diagrams against changes in the scaling parameter decreases in the following order: interaction energies between unlike mers, differences in the scaling temperatures, pressures and densities.

Implementation of a Theory of Edge Detection

This report describes the implementation of a theory of edge detection, proposed by Marr and Hildreth (1979). According to this theory, the image is first processed independently through a set of different size filters, whose shape is the Laplacian of a Gaussian, ***. Zero-crossings in the output of these filters mark the positions of intensity changes at different resolutions. Information about these zero-crossings is then used for deriving a full symbolic description of changes in intensity in the image, called the raw primal sketch. The theory is closely tied with early processing in the human visual systems. In this report, we first examine the critical properties of the initial filters used in the edge detection process, both from a theoretical and practical standpoint. The implementation is then used as a test bed for exploring aspects of the human visual system; in particular, acuity and hyperacuity. Finally, we present some preliminary results concerning the relationship between zero-crossings detected at different resolutions, and some observations relevant to the process by which the human visual system integrates descriptions of intensity changes obtained at different resolutions.

Design thinking - a critical analysis.

The latest buzz phrase to enter the world of design research is “Design Thinking”. But is this anything new and does it really have any practical or theoretical relevance to the design world? Many sceptics believe the term has more to do with business strategy and little to do with the complex process of designing products, services and systems. Moreover, many view the term as misleading and a cheap attempt to piggyback the world of business management onto design. This paper seeks to ask is design thinking anything new? Several authors have explicitly or implicitly articulated the term “Design Thinking” before, such as Peter Rowe’s seminal book “Design Thinking” [1] first published in 1987 and Herbert Simon’s “The Sciences of the Artificial” [2] first published in 1969. In Tim Brown’s “Change by Design” [3], design thinking is thought of as a system of three overlapping spaces rather than a sequence of orderly steps namely inspiration – the problem or opportunity that motivates the search for solutions; ideation – the process of generating, developing and testing ideas; and implementation – the path that leads from the design studio, lab and factory to the market. This paper seeks to examine and critically analyse the tenets of this new design thinking manifesto set against three case studies of modern design practice. As such, the paper will compare design thinking theory with the reality of design in practice.

Understanding nurses’ decision-making when managing weaning from mechanical ventilation: a study of novice and experienced critical care nurses in Scotland and Greece.

Aim and objectives To examine how nurses collect and use cues from respiratory assessment to inform their decisions as they wean patients from ventilatory support. Background Prompt and accurate identification of the patient's ability to sustain reduction of ventilatory support has the potential to increase the likelihood of successful weaning. Nurses' information processing during the weaning from mechanical ventilation has not been well-described. Design A descriptive ethnographic study exploring critical care nurses' decision-making processes when weaning mechanically ventilated patients from ventilatory support in the real setting. Methods Novice and expert Scottish and Greek nurses from two tertiary intensive care units were observed in real practice of weaning mechanical ventilation and were invited to participate in reflective interviews near the end of their shift. Data were analysed thematically using concept maps based on information processing theory. Ethics approval and informed consent were obtained. Results Scottish and Greek critical care nurses acquired patient-centred objective physiological and subjective information from respiratory assessment and previous knowledge of the patient, which they clustered around seven concepts descriptive of the patient's ability to wean. Less experienced nurses required more encounters of cues to attain the concepts with certainty. Subjective criteria were intuitively derived from previous knowledge of patients' responses to changes of ventilatory support. All nurses used focusing decision-making strategies to select and group cues in order to categorise information with certainty and reduce the mental strain of the decision task. Conclusions Nurses used patient-centred information to make a judgment about the patients' ability to wean. Decision-making strategies that involve categorisation of patient-centred information can be taught in bespoke educational programmes for mechanical ventilation and weaning. Relevance to clinical practice Advanced clinical reasoning skills and accurate detection of cues in respiratory assessment by critical care nurses will ensure optimum patient management in weaning mechanical ventilation