Large linear magnetoresistance in a GaAs/AlGaAs heterostructure

We report non-saturating linear magnetoresistance (MR) in a two-dimensional electron system (2DES) at a GaAs/AlGaAs heterointerface in the strongly insulating regime. We achieve this by driving the gate voltage below the pinch-off point of the device and operating it in the non-equilibrium regime with high source-drain bias. Remarkably, the magnitude of MR is as large as 500% per Tesla with respect to resistance at zero magnetic field, thus dwarfing most non-magnetic materials which exhibit this linearity. Its primary advantage over most other materials is that both linearity and the enormous magnitude are retained over a broad temperature range (0.3 K to 10 K), thus making it an attractive candidate for cryogenic sensor applications.

Probing Photoexcited Carriers in a Few-Layer MoS2 Laminate by Time-Resolved Optical Pump-Terahertz Probe Spectroscopy

We report the dynamics of photoinduced carriers in a free-standing MoS2 laminate consisting of a few layers (1-6 layers) using time-resolved optical pump-terahertz probe spectroscopy. Upon photoexcitation with the 800 nm pump pulse, the terahertz conductivity increases due to absorption by the photoinduced charge carriers. The relaxation of the non-equilibrium carriers shows fast as well as slow decay channels, analyzed using a rate equation model incorporating defect-assisted Auger scattering of photoexcited electrons, holes, and excitons. The fast relaxation time occurs due to the capture of electrons and holes by defects via Auger processes, resulting in nonradiative recombination. The slower relaxation arises since the excitons are bound to the defects, preventing the defect-assisted Auger recombination of the electrons and the holes. Our results provide a comprehensive understanding of the non-equilibrium carrier kinetics in a system of unscreened Coulomb interactions, where defect-assisted Auger processes dominate and should be applicable to other 2D systems.

Theoretical insights on the electro-thermal transport properties of monolayer MoS2 with line defects

Two dimensional (2D) materials demonstrate several novel electrical, mechanical, and thermal properties which are quite distinctive to those of their bulk form. Among many others, one important potential application of the 2D material is its use in the field of energy harvesting. Owing to that, here we present a detailed study on electrical as well as thermal transport of monolayer MoS2, in quasi ballistic regime. Besides the perfect monolayer in its pristine form, we also consider various line defects which have been experimentally observed in mechanically exfoliated MoS2 samples. For calculating various parameters related to the electrical transmission, we employ the non-equilibrium Green's function-density functional theory combination. However, to obtain the phonon transmission, we take help of the parametrized Stillinger-Weber potential which can accurately delineate the inter-atomic interactions for the monolayer MoS2. Due to the presence of line defects, we observed significant reductions in both the charge carrier and the phonon transmissions through a monolayer MoS2 flake. Moreover, we also report a comparative analysis showing the temperature dependency of the thermoelectric figure of merit values, as obtained for the perfect as well as the other defective 2D samples. (C) 2016 AIP Publishing LLC.

Dynamics Of Fine Particles In Liquid-Solid Fluidized Beds

On the basis of the Local Equilibrium Model (LEM), fine particles with large Richardson-Zaki exponent n show, under certain conditions during bed expansion and collapse, different dynamic behavior from particles with small n. For an expansion process there may be a concentration discontinuity propagating upward from the distributor, and, on the contrary, for a collapse process there may be a progressively broadening and upward-propagating continuous transition zone instead of discontinuity. The predictions of the bed height variation and the discontinuity trace have been validated experimentally. (c) 2007 Chinese Society of Particuology and Institute of Process Engineering, Chinese Academy of Sciences. Published by Elsevier B.V. All rights reserved.

Thickness And Component Distributions Of Yttrium-Titanium Alloy Films In Electron-Beam Physical Vapor Deposition

Thickness and component distributions of large-area thin films are an issue of international concern in the field of material processing. The present work employs experiments and direct simulation Monte Carlo (DSMC) method to investigate three-dimensional low-density, non-equilibrium jets of yttrium and titanium vapor atoms in an electron-beams physical vapor deposition (EBPVD) system furnished with two or three electron-beams, and obtains their deposition thickness and component distributions onto 4-inch and 6-inch mono-crystal silicon wafers. The DSMC results are found in excellent agreement with our measurements, such as evaporation rates of yttrium and titanium measured in-situ by quartz crystal resonators, deposited film thickness distribution measured by Rutherford backscattering spectrometer (RBS) and surface profilometer and deposited film molar ratio distribution measured by RBS and inductively coupled plasma atomic emission spectrometer (ICP-AES). This can be taken as an indication that a combination of DSMC method with elaborate measurements may be satisfactory for predicting and designing accurately the transport process of EBPVD at the atomic level.

Research on laser welding of cast Ni-based superalloy K418 turbo disk and alloy steel 42CrMo shaft

Exploratory experiments of laser welding cast Ni-based superalloy K418 turbo disk and alloy steel 42CrMo shaft were conducted. Microstructure of the welded seam was characterized by optical microscopy (OM), scanning electron microscopy (SEM), X-ray diffraction (XRD), energy-dispersive spectrometer (EDS). Mechanical properties of the welded seam were evaluated by microhardness and tensile strength testing. The corresponding mechanisms were discussed in detail. Results showed that the laser-welded seam had non-equilibrium solidified microstructures consisting of FeCr0.29Ni0.16C0.06 austenite solid solution dendrites as the dominant and some fine and dispersed Ni3Al gamma' phase and Laves particles as well as little amount of MC short stick or particle-like carbides distributed in the interdendritic regions. The average microhardness of the welded seam was relatively uniform and lower than that of the base metal due to partial dissolution and suppression of the strengthening phase gamma' to some extent. About 88.5% tensile strength of the base metal was achieved in the welded joint because of a non-full penetration welding and the fracture mechanism was a mixture of ductility and brittleness. The existence of some Laves particles in the welded seam also facilitated the initiation and propagation of the microcracks and microvoids and hence, the detrimental effects of the tensile strength of the welded joint. The present results stimulate further investigation on this field. (c) 2006 Elsevier B.V. All rights reserved.

Distribution of Temperature in GlidArc Discharge Plasma

The GlidArc discharge is one of the main generation methods of non-equilibrium plasma near atmospheric pressures. In general, Gliding Arc discharge is driven by gas flow [1] in axial direction or by magnetic field in circumferential direction. [2] In this paper, a GlidArc discharge driven by rotating-gas-flow in circumferential direction is presented. The principle of the plasma generator is analyzed. The distribution of the temperature in axial direction is measured by a digital thermometer for three different gases. The experimental set-up of the GlidArc plasma is shown in Fig.1. It consists of a center electrode, an outside electrode, a power supply and a gas supply. The shortest distance between the electrodes is 2-3 mm. When a power supply with 10000 volts is attached to the electrodes, the arc will be ignited at the shortest distance. The small plasma column is rotated by the rotating gas flow in circumferential direction and then the rotating arc is driven towards the exit of the setup by the gas flow.

Evaluating approaches and technologies for monitoring organic contaminants in the aquatic environment, Ann Arbor, MI, July 21-23, 2006: workshop proceedings

The Alliance for Coastal Technologies (ACT) convened a workshop on Evaluating Approaches and Technologies for Monitoring Organic Contaminants in the Aquatic Environment in Ann Arbor, MI on July 21-23, 2006. The primary objectives of this workshop were to: 1) identify the priority management information needs relative to organic contaminant loading; 2) explore the most appropriate approaches to estimating mass loading; and 3) evaluate the current status of the sensor technology. To meet these objectives, a mixture of leading research scientists, resource managers, and industry representatives were brought together for a focused two-day workshop. The workshop featured four plenary talks followed by breakout sessions in which arranged groups of participants where charged to respond to a series of focused discussion questions. At present, there are major concerns about the inadequacies in approaches and technologies for quantifying mass emissions and detection of organic contaminants for protecting municipal water supplies and receiving waters. Managers use estimates of land-based contaminant loadings to rivers, lakes, and oceans to assess relative risk among various contaminant sources, determine compliance with regulatory standards, and define progress in source reduction. However, accurately quantifying contaminant loading remains a major challenge. Loading occurs over a range of hydrologic conditions, requiring measurement technologies that can accommodate a broad range of ambient conditions. In addition, in situ chemical sensors that provide a means for acquiring continuous concentration measurements are still under development, particularly for organic contaminants that typically occur at low concentrations. Better approaches and strategies for estimating contaminant loading, including evaluations of both sampling design and sensor technologies, need to be identified. The following general recommendations were made in an effort to advance future organic contaminant monitoring: 1. Improve the understanding of material balance in aquatic systems and the relationship between potential surrogate measures (e.g., DOC, chlorophyll, particle size distribution) and target constituents. 2. Develop continuous real-time sensors to be used by managers as screening measures and triggers for more intensive monitoring. 3. Pursue surrogate measures and indicators of organic pollutant contamination, such as CDOM, turbidity, or non-equilibrium partitioning. 4. Develop continuous field-deployable sensors for PCBs, PAHs, pyrethroids, and emerging contaminants of concern and develop strategies that couple sampling approaches with tools that incorporate sensor synergy (i.e., measure appropriate surrogates along with the dissolved organics to allow full mass emission estimation).[PDF contains 20 pages]

The hydro-kinetic theory of liquid helium II

In Part I the kinetic theory of excitations in flowing liquid He II is developed to a higher order than that carried out previously, by Landau and Khalatnikov, in order to demonstrate the existence of non-equilibrium terms of a new nature in the hydrodynamic equations. It is then shown that these terms can lead to spontaneous destabilization in counter currents when the relative velocity of the normal and super fluids exceeds a critical value that depends on the temperature, but not on geometry. There are no adjustable parameters in the theory. The critical velocities are estimated to be in the 14-20 m/sec range for T ≤ 2.0° K, but tend to zero as T → T_λ. The possibility that these critical velocities may be related to the experimentally observed "intrinsic" critical velocities is discussed.

Part II consists of a semi-classical investigation of rotonquantized vortex line interactions. An essentially classical model is used for the collision and the behavior of the roton in the vortex field is investigated in detail. From this model it is possible to derive the HVBK mutual friction terms that appear in the phenomenalogical equations of motion for rotating liquid He II. Estimates of the Hall and Vinen B and B' coefficients are in good agreement with experiments. The claim is made that the theory does not contain any arbitrary adjustable parameters.

Geometric integration applied to moving mesh methods and degenerate Lagrangians

Moving mesh methods (also called r-adaptive methods) are space-adaptive strategies used for the numerical simulation of time-dependent partial differential equations. These methods keep the total number of mesh points fixed during the simulation, but redistribute them over time to follow the areas where a higher mesh point density is required. There are a very limited number of moving mesh methods designed for solving field-theoretic partial differential equations, and the numerical analysis of the resulting schemes is challenging. In this thesis we present two ways to construct r-adaptive variational and multisymplectic integrators for (1+1)-dimensional Lagrangian field theories. The first method uses a variational discretization of the physical equations and the mesh equations are then coupled in a way typical of the existing r-adaptive schemes. The second method treats the mesh points as pseudo-particles and incorporates their dynamics directly into the variational principle. A user-specified adaptation strategy is then enforced through Lagrange multipliers as a constraint on the dynamics of both the physical field and the mesh points. We discuss the advantages and limitations of our methods. The proposed methods are readily applicable to (weakly) non-degenerate field theories---numerical results for the Sine-Gordon equation are presented.

In an attempt to extend our approach to degenerate field theories, in the last part of this thesis we construct higher-order variational integrators for a class of degenerate systems described by Lagrangians that are linear in velocities. We analyze the geometry underlying such systems and develop the appropriate theory for variational integration. Our main observation is that the evolution takes place on the primary constraint and the 'Hamiltonian' equations of motion can be formulated as an index 1 differential-algebraic system. We then proceed to construct variational Runge-Kutta methods and analyze their properties. The general properties of Runge-Kutta methods depend on the 'velocity' part of the Lagrangian. If the 'velocity' part is also linear in the position coordinate, then we show that non-partitioned variational Runge-Kutta methods are equivalent to integration of the corresponding first-order Euler-Lagrange equations, which have the form of a Poisson system with a constant structure matrix, and the classical properties of the Runge-Kutta method are retained. If the 'velocity' part is nonlinear in the position coordinate, we observe a reduction of the order of convergence, which is typical of numerical integration of DAEs. We also apply our methods to several models and present the results of our numerical experiments.

Studies of phonon anharmonicity in solids

Today our understanding of the vibrational thermodynamics of materials at low temperatures is emerging nicely, based on the harmonic model in which phonons are independent. At high temperatures, however, this understanding must accommodate how phonons interact with other phonons or with other excitations. We shall see that the phonon-phonon interactions give rise to interesting coupling problems, and essentially modify the equilibrium and non-equilibrium properties of materials, e.g., thermodynamic stability, heat capacity, optical properties and thermal transport of materials. Despite its great importance, to date the anharmonic lattice dynamics is poorly understood and most studies on lattice dynamics still rely on the harmonic or quasiharmonic models. There have been very few studies on the pure phonon anharmonicity and phonon-phonon interactions. The work presented in this thesis is devoted to the development of experimental and computational methods on this subject.

Modern inelastic scattering techniques with neutrons or photons are ideal for sorting out the anharmonic contribution. Analysis of the experimental data can generate vibrational spectra of the materials, i.e., their phonon densities of states or phonon dispersion relations. We obtained high quality data from laser Raman spectrometer, Fourier transform infrared spectrometer and inelastic neutron spectrometer. With accurate phonon spectra data, we obtained the energy shifts and lifetime broadenings of the interacting phonons, and the vibrational entropies of different materials. The understanding of them then relies on the development of the fundamental theories and the computational methods.

We developed an efficient post-processor for analyzing the anharmonic vibrations from the molecular dynamics (MD) calculations. Currently, most first principles methods are not capable of dealing with strong anharmonicity, because the interactions of phonons are ignored at finite temperatures. Our method adopts the Fourier transformed velocity autocorrelation method to handle the big data of time-dependent atomic velocities from MD calculations, and efficiently reconstructs the phonon DOS and phonon dispersion relations. Our calculations can reproduce the phonon frequency shifts and lifetime broadenings very well at various temperatures.

To understand non-harmonic interactions in a microscopic way, we have developed a numerical fitting method to analyze the decay channels of phonon-phonon interactions. Based on the quantum perturbation theory of many-body interactions, this method is used to calculate the three-phonon and four-phonon kinematics subject to the conservation of energy and momentum, taking into account the weight of phonon couplings. We can assess the strengths of phonon-phonon interactions of different channels and anharmonic orders with the calculated two-phonon DOS. This method, with high computational efficiency, is a promising direction to advance our understandings of non-harmonic lattice dynamics and thermal transport properties.

These experimental techniques and theoretical methods have been successfully performed in the study of anharmonic behaviors of metal oxides, including rutile and cuprite stuctures, and will be discussed in detail in Chapters 4 to 6. For example, for rutile titanium dioxide (TiO₂), we found that the anomalous anharmonic behavior of the B_1g mode can be explained by the volume effects on quasiharmonic force constants, and by the explicit cubic and quartic anharmonicity. For rutile tin dioxide (SnO₂), the broadening of the B_2g mode with temperature showed an unusual concave downwards curvature. This curvature was caused by a change with temperature in the number of down-conversion decay channels, originating with the wide band gap in the phonon dispersions. For silver oxide (Ag₂O), strong anharmonic effects were found for both phonons and for the negative thermal expansion.

Investigation of bulk indium antimonide as heterodyne detector for the submillimeter wavelength region

Bulk n-lnSb is investigated at a heterodyne detector for the submillimeter wavelength region. Two modes or operation are investigated: (1) the Rollin or hot electron bolometer mode (zero magnetic field), and (2) the Putley mode (quantizing magnetic field). The highlight of the thesis work is the pioneering demonstration or the Putley mode mixer at several frequencies. For example, a double-sideband system noise temperature of about 510K was obtained using a 812 GHz methanol laser for the local oscillator. This performance is at least a factor or 10 more sensitive than any other performance reported to date at the same frequency. In addition, the Putley mode mixer achieved system noise temperatures of 250K at 492 GHz and 350K at 625 GHz. The 492 GHz performance is about 50% better and the 625 GHz is about 100% better than previous best performances established by the Rollin-mode mixer. To achieve these results, it was necessary to design a totally new ultra-low noise, room-temperature preamp to handle the higher source impedance imposed by the Putley mode operation. This preamp has considerably less input capacitance than comparably noisy, ambient designs.

In addition to advancing receiver technology, this thesis also presents several novel results regarding the physics of n-lnSb at low temperatures. A Fourier transform spectrometer was constructed and used to measure the submillimeter wave absorption coefficient of relatively pure material at liquid helium temperatures and in zero magnetic field. Below 4.2K, the absorption coefficient was found to decrease with frequency much faster than predicted by Drudian theory. Much better agreement with experiment was obtained using a quantum theory based on inverse-Bremmstrahlung in a solid. Also the noise of the Rollin-mode detector at 4.2K was accurately measured and compared with theory. The power spectrum is found to be well fit by a recent theory of non- equilibrium noise due to Mather. Surprisingly, when biased for optimum detector performance, high purity lnSb cooled to liquid helium temperatures generates less noise than that predicted by simple non-equilibrium Johnson noise theory alone. This explains in part the excellent performance of the Rollin-mode detector in the millimeter wavelength region.

Again using the Fourier transform spectrometer, spectra are obtained of the responsivity and direct detection NEP as a function of magnetic field in the range 20-110 cm^-1. The results show a discernable peak in the detector response at the conduction electron cyclotron resonance frequency tor magnetic fields as low as 3 KG at bath temperatures of 2.0K. The spectra also display the well-known peak due to the cyclotron resonance of electrons bound to impurity states. The magnitude of responsivity at both peaks is roughly constant with magnet1c field and is comparable to the low frequency Rollin-mode response. The NEP at the peaks is found to be much better than previous values at the same frequency and comparable to the best long wavelength results previously reported. For example, a value NEP=4.5x10^-13/Hz^1/2 is measured at 4.2K, 6 KG and 40 cm^-1. Study of the responsivity under conditions of impact ionization showed a dramatic disappearance of the impurity electron resonance while the conduction electron resonance remained constant. This observation offers the first concrete evidence that the mobility of an electron in the N=0 and N=1 Landau levels is different. Finally, these direct detection experiments indicate that the excellent heterodyne performance achieved at 812 GHz should be attainable up to frequencies of at least 1200 GHz.

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.

Programming chemical kinetics: engineering dynamic reaction networks with DNA strand displacement

Over the last century, the silicon revolution has enabled us to build faster, smaller and more sophisticated computers. Today, these computers control phones, cars, satellites, assembly lines, and other electromechanical devices. Just as electrical wiring controls electromechanical devices, living organisms employ "chemical wiring" to make decisions about their environment and control physical processes. Currently, the big difference between these two substrates is that while we have the abstractions, design principles, verification and fabrication techniques in place for programming with silicon, we have no comparable understanding or expertise for programming chemistry.

In this thesis we take a small step towards the goal of learning how to systematically engineer prescribed non-equilibrium dynamical behaviors in chemical systems. We use the formalism of chemical reaction networks (CRNs), combined with mass-action kinetics, as our programming language for specifying dynamical behaviors. Leveraging the tools of nucleic acid nanotechnology (introduced in Chapter 1), we employ synthetic DNA molecules as our molecular architecture and toehold-mediated DNA strand displacement as our reaction primitive.

Abstraction, modular design and systematic fabrication can work only with well-understood and quantitatively characterized tools. Therefore, we embark on a detailed study of the "device physics" of DNA strand displacement (Chapter 2). We present a unified view of strand displacement biophysics and kinetics by studying the process at multiple levels of detail, using an intuitive model of a random walk on a 1-dimensional energy landscape, a secondary structure kinetics model with single base-pair steps, and a coarse-grained molecular model that incorporates three-dimensional geometric and steric effects. Further, we experimentally investigate the thermodynamics of three-way branch migration. Our findings are consistent with previously measured or inferred rates for hybridization, fraying, and branch migration, and provide a biophysical explanation of strand displacement kinetics. Our work paves the way for accurate modeling of strand displacement cascades, which would facilitate the simulation and construction of more complex molecular systems.

In Chapters 3 and 4, we identify and overcome the crucial experimental challenges involved in using our general DNA-based technology for engineering dynamical behaviors in the test tube. In this process, we identify important design rules that inform our choice of molecular motifs and our algorithms for designing and verifying DNA sequences for our molecular implementation. We also develop flexible molecular strategies for "tuning" our reaction rates and stoichiometries in order to compensate for unavoidable non-idealities in the molecular implementation, such as imperfectly synthesized molecules and spurious "leak" pathways that compete with desired pathways.

We successfully implement three distinct autocatalytic reactions, which we then combine into a de novo chemical oscillator. Unlike biological networks, which use sophisticated evolved molecules (like proteins) to realize such behavior, our test tube realization is the first to demonstrate that Watson-Crick base pairing interactions alone suffice for oscillatory dynamics. Since our design pipeline is general and applicable to any CRN, our experimental demonstration of a de novo chemical oscillator could enable the systematic construction of CRNs with other dynamic behaviors.

The generalized second law of thermodynamics in generalized gravity theories

We investigate the generalized second law of thermodynamics (GSL) in generalized theories of gravity. We examine the total entropy evolution with time including the horizon entropy, the non-equilibrium entropy production, and the entropy of all matter, field and energy components. We derive a universal condition to protect the generalized second law and study its validity in different gravity theories. In Einstein gravity (even in the phantom-dominated universe with a Schwarzschild black hole), Lovelock gravity and braneworld gravity, we show that the condition to keep the GSL can always be satisfied. In f ( R) gravity and scalar-tensor gravity, the condition to protect the GSL can also hold because the temperature should be positive, gravity is always attractive and the effective Newton constant should be an approximate constant satisfying the experimental bounds.