Graphene-MoS2 Hybrid Technology for Large-Scale Two-Dimensional Electronics

Two-dimensional (2D) materials have generated great interest in the last few years as a new toolbox for electronics. This family of materials includes, among others, metallic graphene, semiconducting transition metal dichalcogenides (such as MoS2) and insulating Boron Nitride. These materials and their heterostructures offer excellent mechanical flexibility, optical transparency and favorable transport properties for realizing electronic, sensing and optical systems on arbitrary surfaces. In this work, we develop several etch stop layer technologies that allow the fabrication of complex 2D devices and present for the first time the large scale integration of graphene with molybdenum disulfide (MoS2) , both grown using the fully scalable CVD technique. Transistor devices and logic circuits with MoS2 channel and graphene as contacts and interconnects are constructed and show high performances. In addition, the graphene/MoS2 heterojunction contact has been systematically compared with MoS2-metal junctions experimentally and studied using density functional theory. The tunability of the graphene work function significantly improves the ohmic contact to MoS2. These high-performance large-scale devices and circuits based on 2D heterostructure pave the way for practical flexible transparent electronics in the future. The authors acknowledge financial support from the Office of Naval Research (ONR) Young Investigator Program, the ONR GATE MURI program, and the Army Research Laboratory. This research has made use of the MI.

Device Physics of Two-Dimensional Crystalline Materials

Thesis (Ph.D.)--University of Washington, 2016-07

Tile-Based Two-Dimensional Phase Unwrapping for Digital Holography Using a Modular Framework

A variety of physical and biomedical imaging techniques, such as digital holography, interferometric synthetic aperture radar (InSAR), or magnetic resonance imaging (MRI) enable measurement of the phase of a physical quantity additionally to its amplitude. However, the phase can commonly only be measured modulo 2π, as a so called wrapped phase map. Phase unwrapping is the process of obtaining the underlying physical phase map from the wrapped phase. Tile-based phase unwrapping algorithms operate by first tessellating the phase map, then unwrapping individual tiles, and finally merging them to a continuous phase map. They can be implemented computationally efficiently and are robust to noise. However, they are prone to failure in the presence of phase residues or erroneous unwraps of single tiles. We tried to overcome these shortcomings by creating novel tile unwrapping and merging algorithms as well as creating a framework that allows to combine them in modular fashion. To increase the robustness of the tile unwrapping step, we implemented a model-based algorithm that makes efficient use of linear algebra to unwrap individual tiles. Furthermore, we adapted an established pixel-based unwrapping algorithm to create a quality guided tile merger. These original algorithms as well as previously existing ones were implemented in a modular phase unwrapping C++ framework. By examining different combinations of unwrapping and merging algorithms we compared our method to existing approaches. We could show that the appropriate choice of unwrapping and merging algorithms can significantly improve the unwrapped result in the presence of phase residues and noise. Beyond that, our modular framework allows for efficient design and test of new tile-based phase unwrapping algorithms. The software developed in this study is freely available.

High-temperature superconductivity in the two-dimensional t-J model: Gutzwiller wavefunction solution

A systematic diagrammatic expansion for Gutzwiller wavefunctions (DE-GWFs) proposed very recently is used for the description of the superconducting (SC) ground state in the two-dimensional square-lattice t-J model with the hopping electron amplitudes t (and t') between nearest (and next-nearest) neighbors. For the example of the SC state analysis we provide a detailed comparison of the method's results with those of other approaches. Namely, (i) the truncated DE-GWF method reproduces the variational Monte Carlo (VMC) results and (ii) in the lowest (zeroth) order of the expansion the method can reproduce the analytical results of the standard Gutzwiller approximation (GA), as well as of the recently proposed 'grand-canonical Gutzwiller approximation' (called either GCGA or SGA). We obtain important features of the SC state. First, the SC gap at the Fermi surface resembles a d(x2-y2) wave only for optimally and overdoped systems, being diminished in the antinodal regions for the underdoped case in a qualitative agreement with experiment. Corrections to the gap structure are shown to arise from the longer range of the real-space pairing. Second, the nodal Fermi velocity is almost constant as a function of doping and agrees semi-quantitatively with experimental results. Third, we compare the

Bumps and rings in a two-dimensional neural field: splitting and rotational instabilities

In this paper we consider instabilities of localised solutions in planar neural field firing rate models of Wilson-Cowan or Amari type. Importantly we show that angular perturbations can destabilise spatially localised solutions. For a scalar model with Heaviside firing rate function we calculate symmetric one-bump and ring solutions explicitly and use an Evans function approach to predict the point of instability and the shapes of the dominant growing modes. Our predictions are shown to be in excellent agreement with direct numerical simulations. Moreover, beyond the instability our simulations demonstrate the emergence of multi-bump and labyrinthine patterns. With the addition of spike-frequency adaptation, numerical simulations of the resulting vector model show that it is possible for structures without rotational symmetry, and in particular multi-bumps, to undergo an instability to a rotating wave. We use a general argument, valid for smooth firing rate functions, to establish the conditions necessary to generate such a rotational instability. Numerical continuation of the rotating wave is used to quantify the emergent angular velocity as a bifurcation parameter is varied. Wave stability is found via the numerical evaluation of an associated eigenvalue problem.

Green's function retrieval through cross-correlations in a two -dimensional complex reverberating medium

International audience

Radiation from a D-dimensional collision of shock waves: two-dimensional reduction and Carter–Penrose diagram

We analyze the causal structure of the two-dimensional (2D) reduced background used in the perturbative treatment of a head-on collision of two D-dimensional Aichelburg–Sexl gravitational shock waves. After defining all causal boundaries, namely the future light-cone of the collision and the past light-cone of a future observer, we obtain characteristic coordinates using two independent methods. The first is a geometrical construction of the null rays which define the various light cones, using a parametric representation. The second is a transformation of the 2D reduced wave operator for the problem into a hyperbolic form. The characteristic coordinates are then compactified allowing us to represent all causal light rays in a conformal Carter–Penrose diagram. Our construction holds to all orders in perturbation theory. In particular, we can easily identify the singularities of the source functions and of the Green’s functions appearing in the perturbative expansion, at each order, which is crucial for a successful numerical evaluation of any higher order corrections using this method.

Discrete breathers in a two-dimensional hexagonal Fermi-Pasta-Ulam lattice

We consider a two-dimensional Fermi-Pasta-Ulam (FPU) lattice with hexagonal symmetry. Using asymptotic methods based on small amplitude ansatz, at third order we obtain a eduction to a cubic nonlinear Schr{\"o}dinger equation (NLS) for the breather envelope. However, this does not support stable soliton solutions, so we pursue a higher-order analysis yielding a generalised NLS, which includes known stabilising terms. We present numerical results which suggest that long-lived stationary and moving breathers are supported by the lattice. We find breather solutions which move in an arbitrary direction, an ellipticity criterion for the wavenumbers of the carrier wave, symptotic estimates for the breather energy, and a minimum threshold energy below which breathers cannot be found. This energy threshold is maximised for stationary breathers, and becomes vanishingly small near the boundary of the elliptic domain where breathers attain a maximum speed. Several of the results obtained are similar to those obtained for the square FPU lattice (Butt \& Wattis, {\em J Phys A}, {\bf 39}, 4955, (2006)), though we find that the square and hexagonal lattices exhibit different properties in regard to the generation of harmonics, and the isotropy of the generalised NLS equation.

Two-dimensional Stokes flow driven by elliptical paddles

A fast and accurate numerical technique is developed for solving the biharmonic equation in a multiply connected domain, in two dimensions. We apply the technique to the computation of slow viscous flow (Stokes flow) driven by multiple stirring rods. Previously, the technique has been restricted to stirring rods of circular cross section; we show here how the prior method fails for noncircular rods and how it may be adapted to accommodate general rod cross sections, provided only that for each there exists a conformal mapping to a circle. Corresponding simulations of the flow are described, and their stirring properties and energy requirements are discussed briefly. In particular the method allows an accurate calculation of the flow when flat paddles are used to stir a fluid chaotically.

Spin injection in two-dimensional layered materials and local magnetoresistance side-effects

180 p.

Exact solutions for the agricultural and the two-dimensional packing problems

Objectives and study method: The objective of this study is to develop exact algorithms that can be used as management tools for the agricultural production planning and to obtain exact solutions for two of the most well known twodimensional packing problems: the strip packing problem and the bin packing problem. For the agricultural production planning problem we propose a new hierarchical scheme of three stages to improve the current agricultural practices. The objective of the first stage is to delineate rectangular and homogeneous management zones into the farmer’s plots considering the physical and chemical soil properties. This is an important task because the soil properties directly affect the agricultural production planning. The methodology for this stage is based on a new method called “Positions and Covering” that first generates all the possible positions in which the plot can be delineated. Then, we use a mathematical model of linear programming to obtain the optimal physical and chemical management zone delineation of the plot. In the second stage the objective is to determine the optimal crop pattern that maximizes the farmer’s profit taken into account the previous management zones delineation. In this case, the crop pattern is affected by both management zones delineation, physical and chemical. A mixed integer linear programming is used to solve this stage. The objective of the last stage is to determine in real-time the amount of water to irrigate in each crop. This stage takes as input the solution of the crop planning stage, the atmospheric conditions (temperature, radiation, etc.), the humidity level in plots, and the physical management zones of plots, just to name a few. This procedure is made in real-time during each irrigation period. A linear programming is used to solve this problem. A breakthrough happen when we realize that we could propose some adaptations of the P&C methodology to obtain optimal solutions for the two-dimensional packing problem and the strip packing. We empirically show that our methodologies are efficient on instances based on real data for both problems: agricultural and two-dimensional packing problems. Contributions and conclusions: The exact algorithms showed in this study can be used in the making-decision support for agricultural planning and twodimensional packing problems. For the agricultural planning problem, we show that the implementation of the new hierarchical approach can improve the farmer profit between 5.27% until 8.21% through the optimization of the natural resources. An important characteristic of this problem is that the soil properties (physical and chemical) and the real-time factors (climate, humidity level, evapotranspiration, etc.) are incorporated. With respect to the two-dimensional packing problems, one of the main contributions of this study is the fact that we have demonstrate that many of the best solutions founded in literature by others approaches (heuristics approaches) are the optimal solutions. This is very important because some of these solutions were up to now not guarantee to be the optimal solutions.

Efficient photo-induced dielectrophoretic particle trapping on Fe:LiNbO3 for arbitrary two dimensional patterning

1D and 2D patterning of uncharged micro- and nanoparticles via dielectrophoretic forces on photovoltaic z-cut Fe:LiNbO3 have been investigated for the first time. The technique has been successfully applied with dielectric micro-particles of CaCO3 (diameter d = 1-3 ?m) and metal nanoparticles of Al (d = 70 nm). At difference with previous experiments in x- and y-cut, the obtained patterns locally reproduce the light distribution with high fidelity. A simple model is provided to analyse the trapping process. The results show the remarkably good capabilities of this geometry for high quality 2D light-induced dielectrophoretic patterning overcoming the important limitations presented by previous configurations.

Topological insulating phases from two-dimensional nodal loop semimetals

Starting from a minimal model for a two-dimensional nodal loop semimetal, we study the effect of chiral mass gap terms. The resulting Dirac loop anomalous Hall insulator’s Chern number is the phase-winding number of the mass gap terms on the loop.We provide simple lattice models, analyze the topological phases, and generalize a previous index characterizing topological transitions. The responses of the Dirac loop anomalous Hall and quantum spin Hall insulators to a magnetic field’s vector potential are also studied both in weak- and strong-field regimes, as well as the edge states in a ribbon geometry.

Characterization of the volatile profile of Brazilian Merlot wines through comprehensive two dimensional gas chromatography time-of-flight mass spectrometric detection.

Wine aroma is an important characteristic and may be related to certain specific parameters, such as raw material and production process. The complexity of Merlot wine aroma was considered suitable for comprehensive two-dimensional gas chromatography (GCGC), as this technique offers superior performance when compared to one-dimensional gas chromatography (1D-GC). The profile of volatile compounds of Merlot wine was, for the first time, qualitatively analyzed by HS-SPME-GCxGC with a time-of-flight mass spectrometric detector (TOFMS), resulting in 179 compounds tentatively identified by comparison of experimental GCxGC retention indices and mass spectra with literature 1D-GC data and 155 compounds tentatively identified only by mass spectra comparison. A set of GCGC experimental retention indices was also, for the first time, presented for a specific inverse set of columns. Esters were present in higher number (94), followed by alcohols (80), ketones (29), acids (29), aldehydes (23), terpenes (23), lactones (16), furans (14), sulfur compounds (9), phenols (7), pyrroles (5), C13-norisoprenoids (3), and pyrans (2). GCxGC/TOFMS parameters were improved and optimal conditions were: a polar (polyethylene glycol)/medium polar (50% phenyl 50% dimethyl arylene siloxane) column set, oven temperature offset of 10ºC, 7 s as modulation period and 1.4 s of hot pulse duration. Co-elutions came up to 138 compounds in 1D and some of them were resolved in 2D. Among the coeluted compounds, thirty-three volatiles co-eluted in both 1D and 2D and their tentative identification was possible only due to spectral deconvolution. Some compounds that might have important contribution to aroma notes were included in these superimposed peaks. Structurally organized distribution of compounds in the 2D space was observed for esters, aldehydes and ketones, alcohols, thiols, lactones, acids and also inside subgroups, as occurred with esters and alcohols. The Fischer Ratio was useful for establishing the analytes responsible for the main differences between Merlot and non-Merlot wines. Differentiation among Merlot wines and wines of other grape varieties were mainly perceived through the following components: ethyl dodecanoate, 1-hexanol, ethyl nonanoate, ethyl hexanoate, ethyl decanoate, dehydro-2-methyl-3(2H)thiophenone, 3-methyl butanoic acid, ethyl tetradecanoate, methyl octanoate, 1,4 butanediol, and 6-methyloctan-1-ol.

Tailoring the resonance of bilayer graphene sheets by interlayer sp3 bonds

Graphene-based resonators are envisioned to build the ultimate limit of two-dimensional nanoelectromechanical system due to their ultrasensitive detection of mass, force, pressure and charge. However, such application has been greatly impeded by their extremely low quality factor. In the present work, we explore, using the large-scale molecular dynamics simulation, the possibility of tailoring the resonance properties of a bilayer graphene sheet (GS) with interlayer sp3 bonds. For the bilayer GS resonator with interlayer sp3 bonds, we discovered that the sp3 bonds can either degrade or enhance the resonance properties of the resonator depending on their density and location. It is found that the distribution of sp3 bonds only along the edges of either pristine or hydrogenated bilayer GS, leads to a greatly enhanced quality factor. A quality factor of ~1.18×105 is observed for a 3.07×15.31 nm2 bilayer GS resonator with sp3 bonds, which is more than 30 times larger comparing with that of a pristine bilayer GS. The present study demonstrates that the resonance properties of a bilayer GS resonator can be tuned by introducing sp3 bonds. This finding provides a useful guideline for the synthesis of the bilayer GS for its application as a resonator component.