Interface-assisted spintronics: Tailoring at the molecular scale

Organic molecules adsorbed on magnetic surfaces offer the possibility to merge the concepts of molecular electronics with spintronics to build future nanoscale data storage, sensing, and computing multifunctional devices. In order to engineer the functionalities of such hybrid spintronic devices, an understanding of the electronic and magnetic properties of the interface between carbon-based aromatic materials and magnetic surfaces is essential. In this article, we discuss recent progress in the study of spin-dependent chemistry and physics associated with the above molecule-ferromagnet interface by combining state-of-the-art experiments and theoretical calculations. The magnetic properties such as molecular magnetic moment, electronic interface spin-polarization, magnetic anisotropy, and magnetic exchange coupling can be specifically tuned by an appropriate choice of the organic material and the magnetic substrate. These reports suggest a gradual shift in research toward an emerging subfield of interface-assisted molecular spintronics.

Fabrication of large-area PbSe films at the organic-aqueous interface and their near-infrared photoresponse

An organic-aqueous interfacial reaction at room temperature has been employed to synthesize large-area self-assembled films consisting of PbSe single crystallites. The use of the films for the low-cost fabrication of IR-photodetectors has been explored. (111)-oriented single crystallites of PbSe self-assemble to form robust large-area films. The near-infrared photoresponse of the film measured at room temperature showed large responsivity and gain owing to trap-associated mechanisms. Low-cost, mild reaction conditions and tunability of the nature of deposits make the present strategy useful for synthesizing large-area films of functional materials for possible opto-electronic applications.

A spectral multiscale method for wave propagation analysis: Atomistic-continuum coupled simulation

In this paper, we present a new multiscale method which is capable of coupling atomistic and continuum domains for high frequency wave propagation analysis. The problem of non-physical wave reflection, which occurs due to the change in system description across the interface between two scales, can be satisfactorily overcome by the proposed method. We propose an efficient spectral domain decomposition of the total fine scale displacement along with a potent macroscale equation in the Laplace domain to eliminate the spurious interfacial reflection. We use Laplace transform based spectral finite element method to model the macroscale, which provides the optimum approximations for required dynamic responses of the outer atoms of the simulated microscale region very accurately. This new method shows excellent agreement between the proposed multiscale model and the full molecular dynamics (MD) results. Numerical experiments of wave propagation in a 1D harmonic lattice, a 1D lattice with Lennard-Jones potential, a 2D square Bravais lattice, and a 2D triangular lattice with microcrack demonstrate the accuracy and the robustness of the method. In addition, under certain conditions, this method can simulate complex dynamics of crystalline solids involving different spatial and/or temporal scales with sufficient accuracy and efficiency. (C) 2014 Elsevier B.V. All rights reserved.

Nonlinear preconditioning for efficient and accurate interface capturing in simulation of multicomponent compressible flows

Single fluid schemes that rely on an interface function for phase identification in multicomponent compressible flows are widely used to study hydrodynamic flow phenomena in several diverse applications. Simulations based on standard numerical implementation of these schemes suffer from an artificial increase in the width of the interface function owing to the numerical dissipation introduced by an upwind discretization of the governing equations. In addition, monotonicity requirements which ensure that the sharp interface function remains bounded at all times necessitate use of low-order accurate discretization strategies. This results in a significant reduction in accuracy along with a loss of intricate flow features. In this paper we develop a nonlinear transformation based interface capturing method which achieves superior accuracy without compromising the simplicity, computational efficiency and robustness of the original flow solver. A nonlinear map from the signed distance function to the sigmoid type interface function is used to effectively couple a standard single fluid shock and interface capturing scheme with a high-order accurate constrained level set reinitialization method in a way that allows for oscillation-free transport of the sharp material interface. Imposition of a maximum principle, which ensures that the multidimensional preconditioned interface capturing method does not produce new maxima or minima even in the extreme events of interface merger or breakup, allows for an explicit determination of the interface thickness in terms of the grid spacing. A narrow band method is formulated in order to localize computations pertinent to the preconditioned interface capturing method. Numerical tests in one dimension reveal a significant improvement in accuracy and convergence; in stark contrast to the conventional scheme, the proposed method retains its accuracy and convergence characteristics in a shifted reference frame. Results from the test cases in two dimensions show that the nonlinear transformation based interface capturing method outperforms both the conventional method and an interface capturing method without nonlinear transformation in resolving intricate flow features such as sheet jetting in the shock-induced cavity collapse. The ability of the proposed method in accounting for the gravitational and surface tension forces besides compressibility is demonstrated through a model fully three-dimensional problem concerning droplet splash and formation of a crownlike feature. (C) 2014 Elsevier Inc. All rights reserved.

Stabilizing the Native Trimer of HIV-1 Env by Destabilizing the Heterodimeric Interface of the gp41 Postfusion Six-Helix Bundle

The HIV-1 envelope glycoprotein (Env) is a trimer of gp120-gp41 heterodimers and is essential for viral entry. The gp41 subunit in native, prefusion trimeric Env exists in a metastable conformation and attains a stable six-helix bundle (6-HB) conformation comprised of a trimer of N-heptad repeat (NHR) and C-heptad repeat (CHR) heterodimers, that drives the fusion of viral and cellular membranes. We attempted to stabilize native Env trimers by incorporation of mutations at the NHR-CHR interface that disrupt the postfusion 6-HB of gp41. The mutations V570D and I573D stabilize native Env of the HIV-1 JRFL strain and occlude nonneutralizing epitopes to a greater extent than the previously identified I559P mutation that is at the interface of the NHR trimers in the 6-HB. The mutations prevent soluble-CD4 (sCD4)-induced gp120 shedding and 6-HB formation. In the context of cell surface-expressed JRFL Env, introduction of a previously reported additional disulfide between residues A501 and T605 perturbs the native conformation, though this effect is partially alleviated by furin coexpression. The data suggest that positions 570 and 573 are surface proximal in native Env and that the NHR homotrimeric coiled coil in native Env terminates before or close to residue 573. Aspartic acid substitutions at these positions stabilize native trimers through destabilization of the postfusion 6-HB conformation. These mutations can be used to stabilize Env in a DNA vaccine format. IMPORTANCE The major protein on the surface of HIV-1 is the envelope (Env) glycoprotein. Env is a trimer of gp120-gp41 heterodimers. gp120 is involved in receptor/coreceptor binding and gp41 in the fusion of viral and cellular membranes. Like many other viral fusion proteins, the gp41 subunit in native trimeric Env exists in a metastable conformation. gp41 readily forms a stable six-helix bundle (6-HB) conformation comprised of a trimer of N-heptad repeat (NHR) and C-heptad repeat (CHR) heterodimers that drives fusion of viral and cellular membranes. While it is expected that native Env is a good immunogen, its metastability results in exposure of immunodominant nonneutralizing epitopes. In the present study, we stabilize native Env trimers by incorporation of a number of different mutations at the NHR-CHR interface that disrupt the postfusion 6-HB of gp41. The stabilized constructs described here can be incorporated into DNA vaccine candidates.

Finite element analysis of free vibration of the delaminated composite plate with variable kinematic multilayered plate elements

Composite laminates are prone to delamination. Implementation of delamination in the Carrera Unified Formulation frame work using nine noded quadrilateral MITC9 element is discussed in this article. MITC9 element is devoid of shear locking and membrane locking. Delaminated as well as healthy structure is analyzed for free mode vibration. The results from the present work are compared with the available experimental or/and research article or/and the three dimensional finite element simulations. The effect of different kinds and different percentages of area of delamination on the first three natural frequencies of the structure is discussed. The presence of open-mode delamination mode shape for large delaminations within the first three natural frequencies is discussed. Also, the switching of places between the second bending mode, with that of the first torsional mode frequency is discussed. Results obtained from different ordered theories are compared in the presence of delamination. Advantage of layerwise theories as compared to equivalent single layer theories for very large delaminations is stated. The effect of different kinds of delamination and their effect on the second bending and first torsional mode shape is discussed. (C) 2014 Elsevier Ltd. All rights reserved.

Numerical modeling of the non-isothermal liquid droplet impact on a hot solid substrate

The heat transfer from a solid phase to an impinging non-isothermal liquid droplet is studied numerically. A new approach based on an arbitrary Lagrangian-Eulerian (ALE) finite element method for solving the incompressible Navier Stokes equations in the liquid and the energy equation within the solid and the liquid is presented. The novelty of the method consists in using the ALE-formulation also in the solid phase to guarantee matching grids along the liquid solid interface. Moreover, a new technique is developed to compute the heat flux without differentiating the numerical solution. The free surface and the liquid solid interface of the droplet are represented by a moving mesh which can handle jumps in the material parameter and a temperature dependent surface tension. Further, the application of the Laplace-Beltrami operator technique for the curvature approximation allows a natural inclusion of the contact angle. Numerical simulation for varying Reynold, Weber, Peclet and Biot numbers are performed to demonstrate the capabilities of the new approach. (C) 2014 Elsevier Ltd. All rights reserved.

Interface-assisted molecular spintronics

Molecular spintronics, a field that utilizes the spin state of organic molecules to develop magneto-electronic devices, has shown an enormous scientific activity for more than a decade. But, in the last couple of years, new insights in understanding the fundamental phenomena of molecular interaction on magnetic surfaces, forming a hybrid interface, are presenting a new pathway for developing the subfield of interface-assisted molecular spintronics. The recent exploration of such hybrid interfaces involving carbon based aromatic molecules shows a significant excitement and promise over the previously studied single molecular magnets. In the above new scenario, hybridization of the molecular orbitals with the spin-polarized bands of the surface creates new interface states with unique electronic and magnetic character. This study opens up a molecular-genome initiative in designing new handles to functionalize the spin dependent electronic properties of the hybrid interface to construct spin-functional tailor-made devices. Through this article, we review this subject by presenting a fundamental understanding of the interface spin-chemistry and spin-physics by taking support of advanced computational and spectroscopy tools to investigate molecular spin responses with demonstration of new interface phenomena. Spin-polarized scanning tunneling spectroscopy is favorably considered to be an important tool to investigate these hybrid interfaces with intra-molecular spatial resolution. Finally, by addressing some of the recent findings, we propose novel device schemes towards building interface tailored molecular spintronic devices for applications in sensor, memory, and quantum computing.

Mixed finite elements for electromagnetic analysis

The occurrence of spurious solutions is a well-known limitation of the standard nodal finite element method when applied to electromagnetic problems. The two commonly used remedies that are used to address this problem are (i) The addition of a penalty term with the penalty factor based on the local dielectric constant, and which reduces to a Helmholtz form on homogeneous domains (regularized formulation); (ii) A formulation based on a vector and a scalar potential. Both these strategies have some shortcomings. The penalty method does not completely get rid of the spurious modes, and both methods are incapable of predicting singular eigenvalues in non-convex domains. Some non-zero spurious eigenvalues are also predicted by these methods on non-convex domains. In this work, we develop mixed finite element formulations which predict the eigenfrequencies (including their multiplicities) accurately, even for nonconvex domains. The main feature of the proposed mixed finite element formulation is that no ad-hoc terms are added to the formulation as in the penalty formulation, and the improvement is achieved purely by an appropriate choice of finite element spaces for the different variables. We show that the formulation works even for inhomogeneous domains where `double noding' is used to enforce the appropriate continuity requirements at an interface. For two-dimensional problems, the shape of the domain can be arbitrary, while for the three-dimensional ones, with our current formulation, only regular domains (which can be nonconvex) can be modeled. Since eigenfrequencies are modeled accurately, these elements also yield accurate results for driven problems. (C) 2014 Elsevier Ltd. All rights reserved.

Finite element simulations of notch tip fields in magnesium single crystals

Recent experiments using three point bend specimens of Mg single crystals have revealed that tensile twins of {10 (1) over bar2}-type form profusely near a notch tip and enhance the fracture toughness through large plastic dissipation. In this work, 3D finite element simulations of these experiments are carried out using a crystal plasticity framework which includes slip and twinning to gain insights on the mechanics of fracture. The predicted load-displacement curves, slip and tensile twinning activities from finite element analysis corroborate well with the experimental observations. The numerical results are used to explore the 3D nature of the crack tip stress, plastic slip and twin volume fraction distributions near the notch root. The occurrence of tensile twinning is rationalized from the variation of normal stress ahead of the notch tip. Further, deflection of the crack path at twin-twin intersections observed in the experiments is examined from an energy standpoint by modeling discrete twins close to the notch root.

Large deformation dynamic finite element analysis of delaminated composite plates using contact-impact conditions

A new C-0 composite plate finite element based on Reddy's third order theory is used for large deformation dynamic analysis of delaminated composite plates. The inter-laminar contact is modeled with an augmented Lagrangian approach. Numerical results show that the widely used ``unconditionally stable'' beta-Newmark method presents instability problems in the transient simulation of delaminated composite plate structures with large deformation. To overcome this instability issue, an energy and momentum conserving composite implicit time integration scheme presented by Bathe and Baig is used. It is found that a proper selection of the penalty parameter is very crucial in the contact simulation. (C) 2014 Elsevier Ltd. All rights reserved.

A general methodology for calculating mixed mode stress intensity factors and fracture toughness of solder joints with interfacial cracks

Solder joints in electronic packages undergo thermo-mechanical cycling, resulting in nucleation of micro-cracks, especially at the solder/bond-pad interface, which may lead to fracture of the joints. The fracture toughness of a solder joint depends on material properties, process conditions and service history, as well as strain rate and mode-mixity. This paper reports on a methodology for determining the mixed-mode fracture toughness of solder joints with an interfacial starter-crack, using a modified compact mixed mode (CMM) specimen containing an adhesive joint. Expressions for stress intensity factor (K) and strain energy release rate (G) are developed, using a combination of experiments and finite element (FE) analysis. In this methodology, crack length dependent geometry factors to convert for the modified CMM sample are first obtained via the crack-tip opening displacement (CTOD)-based linear extrapolation method to calculate the under far-field mode I and II conditions (f(1a) and f(2a)), (ii) generation of a master-plot to determine a(c), and (iii) computation of K and G to analyze the fracture behavior of joints. The developed methodology was verified using J-integral calculations, and was also used to calculate experimental fracture toughness values of a few lead-free solder-Cu joints. (C) 2014 Elsevier Ltd. All rights reserved.

Modelling of transient heat conduction with diffuse interface methods

We present a survey on different numerical interpolation schemes used for two-phase transient heat conduction problems in the context of interface capturing phase-field methods. Examples are general transport problems in the context of diffuse interface methods with a non-equal heat conductivity in normal and tangential directions to the interface. We extend the tonsorial approach recently published by Nicoli M et al (2011 Phys. Rev. E 84 1-6) to the general three-dimensional (3D) transient evolution equations. Validations for one-dimensional, two-dimensional and 3D transient test cases are provided, and the results are in good agreement with analytical and numerical reference solutions.

Formation of hierarchical structures of Fe2O3 by the liquid-liquid interface technique

We report the formation of dendritic hierarchical structures of alpha-Fe2O3 and nanostructures of Fe2O3 by the simple liquid-liquid interface method. The morphology of thin films determined by high-resolution scanning electron microscopy shows nanorods, nanosheets and dendritic Fe2O3. The identification of phases of iron oxide structures is carried out by using XRD and XPS studies. XRD and XPS measurements point out the highly crystalline dendritic alpha-Fe2O3 phase and the mixed phase of alpha- and gamma-Fe2O3 nanostructures. The magnetic measurement also suggests the presence of a mixed phase in the sample grown for 72 hours.

The effect of gamma-gamma ` interface on the tensile and shear strengths of nickel-based superalloys: A first-principles study

First-principles density functional theory has been used to evaluate the shear and cleavage strength in terms of Griffith work and generalized stacking fault energy (GSF) of (001) plane for gamma, gamma' and gamma-gamma' system as a function of distance from the gamma/gamma' interface. Calculation of Griffith work suggests higher cleavage energy for bulk gamma as compared to gamma' while the GSF calculation suggests higher shear strength for bulk gamma' as compared to gamma. It has been found that the shear strength of the cubic plane of the gamma/gamma' interface is marginally lower than those of bulk gamma and gamma' phases. (C) 2014 Elsevier B.V. All rights reserved.