Topology Optimization of Micromachined Structures with Surface Micromachining Manufacturing Constraints

Topology optimization methods have been shown to have extensive application in the design of microsystems. However, their utility in practical situations is restricted to predominantly planar configurations due to the limitations of most microfabrication techniques in realizing structures with arbitrary topologies in the direction perpendicular to the substrate. This study addresses the problem of synthesizing optimal topologies in the out-of-plane direction while obeying the constraints imposed by surface micromachining. A new formulation that achieves this by defining a design space that implicitly obeys the manufacturing constraints with a continuous design parameterization is presented in this paper. This is in contrast to including manufacturing cost in the objective function or constraints. The resulting solutions of the new formulation obtained with gradient-based optimization directly provide the photolithographic mask layouts. Two examples that illustrate the approach for the case of stiff structures are included.

Guided wave based detection of damage in honeycomb core sandwich structures

We report on the Lamb wave type guided wave propagation in honeycomb core sandwich structures. An experimental study supported by theoretical evaluation of the guided wave characteristics is presented that proves the potential of Lamb wave type guided wave for detection of damage in sandwich structures. A sandwich panel is fabricated with planar dimension of 600 mm x 600 mm, having a core thickness of 7 mm, cell size of 5 mm and 0.1 mm thick aluminum face sheets. Thin piezoelectric patch actuators and sensors are used to excite and sense a frequency band limited guided wave with a central frequency. A linear phased array of piezoelectric patch actuators is used to achieve higher signal strength and directivity. Group velocity dispersion curves and corresponding frequency response of sensed signal are obtained experimentally. Linearity between the excitation signal amplitude and the corresponding sensed signal amplitude is found for certain range of parameters. The nature of damping present in the sandwich panel is monitored by measuring the sensor signal amplitude at various different distances measured from the center of the linear phased array. Indentation and low velocity impact induced damages of increasing diameter covering several honeycomb cells are created. Crushing of honeycomb core with rupture of face sheet is observed while introducing the damage. The damages are then detected experimentally by pitch-catch interrogation with guided waves and wavelet transform of the sensed signal. Signal amplitudes are analyzed for various different sizes of damages to differentiate the damage size/severity. Monotonic changes in the sensor signal amplitude due to increase in the damage size has been established successfully. With this approach it is possible to locate and monitor the damages with the help of phased array and by tracking the wave packets scattered from the damages. (C) 2012 Elsevier Ltd. All rights reserved.

Comparison of tertiary structures of proteins in protein-protein complexes with unbound forms suggests prevalence of allostery in signalling proteins

Abstract: Background: Most signalling and regulatory proteins participate in transient protein-protein interactions during biological processes. They usually serve as key regulators of various cellular processes and are often stable in both protein-bound and unbound forms. Availability of high-resolution structures of their unbound and bound forms provides an opportunity to understand the molecular mechanisms involved. In this work, we have addressed the question "What is the nature, extent, location and functional significance of structural changes which are associated with formation of protein-protein complexes?" Results: A database of 76 non-redundant sets of high resolution 3-D structures of protein-protein complexes, representing diverse functions, and corresponding unbound forms, has been used in this analysis. Structural changes associated with protein-protein complexation have been investigated using structural measures and Protein Blocks description. Our study highlights that significant structural rearrangement occurs on binding at the interface as well as at regions away from the interface to form a highly specific, stable and functional complex. Notably, predominantly unaltered interfaces interact mainly with interfaces undergoing substantial structural alterations, revealing the presence of at least one structural regulatory component in every complex. Interestingly, about one-half of the number of complexes, comprising largely of signalling proteins, show substantial localized structural change at surfaces away from the interface. Normal mode analysis and available information on functions on some of these complexes suggests that many of these changes are allosteric. This change is largely manifest in the proteins whose interfaces are altered upon binding, implicating structural change as the possible trigger of allosteric effect. Although large-scale studies of allostery induced by small-molecule effectors are available in literature, this is, to our knowledge, the first study indicating the prevalence of allostery induced by protein effectors. Conclusions: The enrichment of allosteric sites in signalling proteins, whose mutations commonly lead to diseases such as cancer, provides support for the usage of allosteric modulators in combating these diseases.

Small-Angle Neutron-Scattering Studies of Mixed Micellar Structures Made of Dimeric Surfactants Having Imidazolium and Ammonium Headgroups

Planar imidazolium cation based gemini surfactants 16-Im-n-Im-16], 2Br(-) (where n = 2, 3, 4, 5, 6, 8, 10, and 12), exhibit different morphologies and internal packing arrangements by adopting different supramolecular assemblies in aqueous media depending on their number of spacer methylene units (CH2)(n). Detailed measurements of the small-angle neutron-scattering (SANS) cross sections from different imidazolium-based surfactant micelles in aqueous media (D2O) are reported. The SANS data, containing the information of aggregation behavior of such surfactants in the molecular level, have been analyzed on the basis of the Hayter and Penfold model for the macro ion solution to compute the interparticle structure factor S(Q) taking into account the screened Coulomb interactions between the dimeric surfactant micelles. The characteristic changes in the SANS spectra of the dimeric surfactant with n = 4 due to variation of temperature have also been investigated. These data are then compared with the SANS characterization data of the corresponding gemini micelles containing tetrahedral ammonium ion based polar headgroups. The critical micellar concentration of each surfactant micelle (cmc) has been determined using pyrene as an extrinsic fluorescence probe. The variation of cmc as a function of spacer chain length has been explained in terms of conformational variation and progressive looping of the spacer into the micellar interior upon increasing the n values. Small-angle neutron-scattering (SANS) cross sections from different mixed micelles composed of surfactants with ammonium headgroups, 16-A(0), 16-Am-n-Am-16], 2Br(-) (where n = 4), 16-I-0, and 16-Im-n-Im-16], 2Br(-) (where n = 4), in aqueous media (D2O) have also been analyzed. The aggregate composition matches with that predicted from the ideal mixing model.

Planar triazinium cations from vanadyl-mediated ring cyclizations: the thiazole species for efficient nuclear staining and photocytotoxicity

Planar triazinium cationic species from vanadyl-assisted cyclization of 1-(2-thiazolylazo)-2-naphthol (H-TAN, 1), 1-(2-pyridylazo)-2-naphthol (H-PAN, 2), 2-(2'-thiazolylazo)-p-cresol (H-TAC, 3) and 6-(2'-thiazolylazo)- resorcinol (H-TAR, 5) were prepared and characterized. A dioxovanadium(V) species VO2(TAR)] (4) was also isolated. Compounds 1, 2 and 4 were structurally characterized. Both 1 and 2 have planar structures. Complex 4 has (VO3N2)-O-V coordination geometry. The cyclised triazinium compound forms a radical species within -0.06 to -0.29 V vs. SCE in DMF-0.1 M tetrabutylammonium perchlorate with a second response due to formation of an anionic species. A confocal microscopic study showed higher nuclear uptake for 1 having a fused thiazole moiety than 2 with a fused pyridine ring. The compounds showed a partial intercalative mode of binding to calf thymus DNA. Compound 1 showed plasmid DNA photo-cleavage activity under argon and photocytotoxicity in HeLa and MCF-7 cells with IC50 values of 15.1 and 3.4 mu M respectively in visible light of 400-700 nm, while being essentially non-toxic in the dark with IC50 values of 90.4 and 21.9 mu M. ATDDFT study was done to rationalize the experimental data.

Matrix crack detection in spatially random composite structures using fractal dimension

Fractal dimension based damage detection method is studied for a composite structure with random material properties. A composite plate with localized matrix crack is considered. Matrix cracks are often seen as the initial damage mechanism in composites. Fractal dimension based method is applied to the static deformation curve of the structure to detect localized damage. Static deflection of a cantilevered composite plate under uniform loading is calculated using the finite element method. Composite material shows spatially varying random material properties because of complex manufacturing processes. Spatial variation of material property is represented as a two dimensional homogeneous Gaussian random field. Karhunen-Loeve (KL) expansion is used to generate a random field. The robustness of fractal dimension based damage detection methods is studied considering the composite plate with spatial variation in material properties.

Study of structures, energies and vibrational frequencies of (O-2)(n)(+) (n=2-5) clusters by GGA and meta-GGA density functional methods

Using Generalized Gradient Approximation (GGA) and meta-GGA density functional methods, structures, binding energies and harmonic vibrational frequencies for the clusters O-4(+), O-6(+), O-8(+) and O-10(+) have been calculated. The stable structures of O-4(+), O-6(+), O-8(+) and O-10(+) have point groups D-2h, D-3h, D-4h, and D-5h optimized on the quartet, sextet, octet and dectet potential energy surfaces, respectively. Rectangular (D-2h) O-4(+) has been found to be more stable compared to trans-planar (C-2h) on the quartet potential energy surface. Cyclic structure (D-3h) of CA cluster ion has been calculated to be more stable than other structures. Binding energy (B.E.) of the cyclic O-6(+) is in good agreement with experimental measurement. The zero-point corrected B.E. of O-8(+) with D4h symmetry on the octet potential energy surface and zero-point corrected B.E. of O-10(+) with D-5h symmetry on the dectet potential energy surface are also in good agreement with experimental values. The B.E. value for O-4(+) is close to the experimental value when single point energy is calculated by Brueckner coupled-cluster method, BD(T). (C) 2014 Elsevier B.V. All rights reserved.

Smaller is Plastic: Polymorphic Structures and Mechanism of Deformation in Nanoscale hcp Metals

Using first-principles calculations, we establish the existence of highly-stable polymorphs of hcp metals (Ti, Mg, Be, La and Y) with nanoscale structural periodicity. They arise from heterogeneous deformation of the hcp structure occurring in response to large shear stresses localized at the basal planes separated by a few nanometers. Through Landau theoretical analysis, we show that their stability derives from nonlinear coupling between strains at different length scales. Such multiscale hyperelasticity and long-period structures constitute a new mechanism of size-dependent plasticity and its enhancement in nanoscale hcp metals.

Air microjet system for non-contact force application and the actuation of micro-structures

We demonstrate a non-contact technique to apply calibrated and localized forces in the micro-Newton to milli-Newton range using an air microjet. An electromagnetically actuated diaphragm controlled by a signal generator is used to generate the air microjet. With a nozzle diameter of 150 mu m, the microjet diameter was maintained to a maximum of 1 mm at a distance of 5 mm from the nozzle. The force generated by the microjet was measured using a commercial force sensor to determine the velocity profile of the jet. Axial flow velocities of up to 25 m s(-1) were obtained at distances as long as 6 mm. The microjet exerted a force up to 1 mu N on a poly dimethyl siloxane (PDMS) micropillar (50 mu m in diameter, 157 mu m in height) and 415 mu N on a PDMS membrane (3 mm in diameter, 28 mu m thick). We also demonstrate that from a distance of 6 mm our microjet can exert a peak pressure of 187 Pa with a total force of about 84 mu N on a flat surface with 8 V operating voltage. Out of the cleanroom fabrication and robust design make this system cost effective and durable.

Efficiency measurement of GaN-based quantum well and light-emitting diode structures grown on silicon substrates

The optical efficiency of GaN-based multiple quantum well (MQW) and light emitting diode (LED) structures grown on Si(111) substrates by metal-organic vapor phase epitaxy was measured and compared with equivalent structures on sapphire. The crystalline quality of the LED structures was comprehensively characterized using x-ray diffraction, atomic force microscopy, and plan-view transmission electron microscopy. A room temperature photoluminescence (PL) internal quantum efficiency (IQE) as high as 58% has been achieved in an InGaN/GaN MQW on Si, emitting at 460 nm. This is the highest reported PL-IQE of a c-plane GaN-based MQW on Si, and the radiative efficiency of this sample compares well with similar structures grown on sapphire. Processed LED devices on Si also show good electroluminescence (EL) performance, including a forward bias voltage of ∼3.5 V at 20 mA and a light output power of 1 mW at 45 mA from a 500 ×500 μm2 planar device without the use of any additional techniques to enhance the output coupling. The extraction efficiency of the LED devices was calculated, and the EL-IQE was then estimated to have a maximum value of 33% at a current density of 4 A cm-2, dropping to 30% at a current density of 40 A cm-2 for a planar LED device on Si emitting at 455 nm. The EL-IQE was clearly observed to increase as the structural quality of the material increased for devices on both sapphire and Si substrates. © 2011 American Institute of Physics.

The steady-state response of multidegree-of-freedom systems with a spatially localized nonlinearity

This thesis is concerned with the dynamic response of a General multidegree-of-freedom linear system with a one dimensional nonlinear constraint attached between two points. The nonlinear constraint is assumed to consist of rate-independent conservative and hysteretic nonlinearities and may contain a viscous dissipation element. The dynamic equations for general spatial and temporal load distributions are derived for both continuous and discrete systems. The method of equivalent linearization is used to develop equations which govern the approximate steady-state response to generally distributed loads with harmonic time dependence.

The qualitative response behavior of a class of undamped chainlike structures with a nonlinear terminal constraint is investigated. It is shown that the hardening or softening behavior of every resonance curve is similar and is determined by the properties of the constraint. Also examined are the number and location of resonance curves, the boundedness of the forced response, the loci of response extrema, and other characteristics of the response. Particular consideration is given to the dependence of the response characteristics on the properties of the linear system, the nonlinear constraint, and the load distribution.

Numerical examples of the approximate steady-state response of three structural systems are presented. These examples illustrate the application of the formulation and qualitative theory. It is shown that disconnected response curves and response curves which cross are obtained for base excitation of a uniform shear beam with a cubic spring foundation. Disconnected response curves are also obtained for the steady-state response to a concentrated load of a chainlike structure with a hardening hysteretic constraint. The accuracy of the approximate response curves is investigated.

First principles based multiparadigm modeling of electronic structures and dynamics

Electronic structures and dynamics are the key to linking the material composition and structure to functionality and performance.

An essential issue in developing semiconductor devices for photovoltaics is to design materials with optimal band gaps and relative positioning of band levels. Approximate DFT methods have been justified to predict band gaps from KS/GKS eigenvalues, but the accuracy is decisively dependent on the choice of XC functionals. We show here for CuInSe₂ and CuGaSe₂, the parent compounds of the promising CIGS solar cells, conventional LDA and GGA obtain gaps of 0.0-0.01 and 0.02-0.24 eV (versus experimental values of 1.04 and 1.67 eV), while the historically first global hybrid functional, B3PW91, is surprisingly the best, with band gaps of 1.07 and 1.58 eV. Furthermore, we show that for 27 related binary and ternary semiconductors, B3PW91 predicts gaps with a MAD of only 0.09 eV, which is substantially better than all modern hybrid functionals, including B3LYP (MAD of 0.19 eV) and screened hybrid functional HSE06 (MAD of 0.18 eV).

The laboratory performance of CIGS solar cells (> 20% efficiency) makes them promising candidate photovoltaic devices. However, there remains little understanding of how defects at the CIGS/CdS interface affect the band offsets and interfacial energies, and hence the performance of manufactured devices. To determine these relationships, we use the B3PW91 hybrid functional of DFT with the AEP method that we validate to provide very accurate descriptions of both band gaps and band offsets. This confirms the weak dependence of band offsets on surface orientation observed experimentally. We predict that the CBO of perfect CuInSe₂/CdS interface is large, 0.79 eV, which would dramatically degrade performance. Moreover we show that band gap widening induced by Ga adjusts only the VBO, and we find that Cd impurities do not significantly affect the CBO. Thus we show that Cu vacancies at the interface play the key role in enabling the tunability of CBO. We predict that Na further improves the CBO through electrostatically elevating the valence levels to decrease the CBO, explaining the observed essential role of Na for high performance. Moreover we find that K leads to a dramatic decrease in the CBO to 0.05 eV, much better than Na. We suggest that the efficiency of CIGS devices might be improved substantially by tuning the ratio of Na to K, with the improved phase stability of Na balancing phase instability from K. All these defects reduce interfacial stability slightly, but not significantly.

A number of exotic structures have been formed through high pressure chemistry, but applications have been hindered by difficulties in recovering the high pressure phase to ambient conditions (i.e., one atmosphere and room temperature). Here we use dispersion-corrected DFT (PBE-ulg flavor) to predict that above 60 GPa the most stable form of N₂O (the laughing gas in its molecular form) is a 1D polymer with an all-nitrogen backbone analogous to cis-polyacetylene in which alternate N are bonded (ionic covalent) to O. The analogous trans-polymer is only 0.03-0.10 eV/molecular unit less stable. Upon relaxation to ambient conditions both polymers relax below 14 GPa to the same stable non-planar trans-polymer, accompanied by possible electronic structure transitions. The predicted phonon spectrum and dissociation kinetics validate the stability of this trans-poly-NNO at ambient conditions, which has potential applications as a new type of conducting polymer with all-nitrogen chains and as a high-energy oxidizer for rocket propulsion. This work illustrates in silico materials discovery particularly in the realm of extreme conditions.

Modeling non-adiabatic electron dynamics has been a long-standing challenge for computational chemistry and materials science, and the eFF method presents a cost-efficient alternative. However, due to the deficiency of FSG representation, eFF is limited to low-Z elements with electrons of predominant s-character. To overcome this, we introduce a formal set of ECP extensions that enable accurate description of p-block elements. The extensions consist of a model representing the core electrons with the nucleus as a single pseudo particle represented by FSG, interacting with valence electrons through ECPs. We demonstrate and validate the ECP extensions for complex bonding structures, geometries, and energetics of systems with p-block character (C, O, Al, Si) and apply them to study materials under extreme mechanical loading conditions.

Despite its success, the eFF framework has some limitations, originated from both the design of Pauli potentials and the FSG representation. To overcome these, we develop a new framework of two-level hierarchy that is a more rigorous and accurate successor to the eFF method. The fundamental level, GHA-QM, is based on a new set of Pauli potentials that renders exact QM level of accuracy for any FSG represented electron systems. To achieve this, we start with using exactly derived energy expressions for the same spin electron pair, and fitting a simple functional form, inspired by DFT, against open singlet electron pair curves (H₂ systems). Symmetric and asymmetric scaling factors are then introduced at this level to recover the QM total energies of multiple electron pair systems from the sum of local interactions. To complement the imperfect FSG representation, the AMPERE extension is implemented, and aims at embedding the interactions associated with both the cusp condition and explicit nodal structures. The whole GHA-QM+AMPERE framework is tested on H element, and the preliminary results are promising.

Diffusion and chemical reaction in the porous structures of solid oxide fuel cells

The Rolls-Royce Integrated-Planar Solid Oxide Fuel Cell (IP-SOFC) consists of ceramic modules which have electrochemical cells printed on the outer surfaces. The cathodes are the outermost layer of each cell and are supplied with oxygen from air flowing over the outside of the module. The anodes are in direct contact with the ceramic structure and are supplied with fuel from internal gas channels. Natural gas is reformed into hydrogen for use by the fuel cells in a separate reformer module of similar design except that the fuel cells are replaced by a reforming catalyst layer. The performance of the modules is intrinsically linked to the behaviour of the gas flows within their porous structures. Because the porous layers are very thin, a one-dimensional flow model provides a good representation of the flow property variations between fuel channel and fuel cell or reforming catalyst. The multi-component convective-diffusive flows are simulated using a new theory of flow in porous material, the Cylindrical Pore Interpolation Model. The effects of the catalysed methane reforming and water-gas shift chemical reactions are also considered using appropriate kinetic models. It is found that the shift reaction, which is catalysed by the anode material, has certain beneficial effects on the fuel cell module performance. In the reformer module it was found that the flow resistance of the porous support structure makes it difficult to sustain a high methane conversion rate. Although the analysis is based on IP-SOFC geometry, the modelling approach and general conclusions are applicable to other types of SOFC.

The formation and electronic structures of 3d transition-metal atoms doped in silicon nanowires

Using first-principles methods, we systematically study the mechanism of defect formation and electronic structures for 3d transition-metal impurities (V, Cr, Mn, Fe, and Co) doped in silicon nanowires. We find that the formation energies of 3d transition-metal impurities with electrons or holes at the defect levels always increase as the diameters of silicon nanowires decrease, which suggests that self-purification, i.e., the difficulty of doping in silicon nanowires, should be an intrinsic effect. The calculated results show that the defect formation energies of Mn and Fe impurities are lower than those of V, Cr, and Co impurities in silicon nanowires. It indicates that Mn and Fe can easily occupy substitutional site in the interior of silicon nanowires. Moreover, they have larger localized moments, which means that they are good candidates for Si-based dilute magnetic semiconductor nanowires. The doping of Mn and Fe atom in silicon nanowires introduces a pair of energy levels with t(2) symmetry. One of which is dominated by 3d electrons of Mn or Fe, and the other by neighboring dangling bonds of Si vacancies. In addition, a set of nonbonding states localized on the transition-metal atom with e symmetry is also introduced. (C) 2008 American Institute of Physics. [DOI: 10.1063/1.3000445]