Growth, modification and integration of carbon nanotubes into molecular electronics

Molecules are the smallest possible elements for electronic devices, with active elements for such devices typically a few Angstroms in footprint area. Owing to the possibility of producing ultrahigh density devices, tremendous effort has been invested in producing electronic junctions by using various types of molecules. The major issues for molecular electronics include (1) developing an effective scheme to connect molecules with the present micro- and nano-technology, (2) increasing the lifetime and stabilities of the devices, and (3) increasing their performance in comparison to the state-of-the-art devices. In this work, we attempt to use carbon nanotubes (CNTs) as the interconnecting nanoelectrodes between molecules and microelectrodes. The ultimate goal is to use two individual CNTs to sandwich molecules in a cross-bar configuration while having these CNTs connected with microelectrodes such that the junction displays the electronic character of the molecule chosen. We have successfully developed an effective scheme to connect molecules with CNTs, which is scalable to arrays of molecular electronic devices. To realize this far reaching goal, the following technical topics have been investigated. 1. Synthesis of multi-walled carbon nanotubes (MWCNTs) by thermal chemical vapor deposition (T-CVD) and plasma-enhanced chemical vapor deposition (PECVD) techniques (Chapter 3). We have evaluated the potential use of tubular and bamboo-like MWCNTs grown by T-CVD and PE-CVD in terms of their structural properties. 2. Horizontal dispersion of MWCNTs with and without surfactants, and the integration of MWCNTs to microelectrodes using deposition by dielectrophoresis (DEP) (Chapter 4). We have systematically studied the use of surfactant molecules to disperse and horizontally align MWCNTs on substrates. In addition, DEP is shown to produce impurityfree placement of MWCNTs, forming connections between microelectrodes. We demonstrate the deposition density is tunable by both AC field strength and AC field frequency. 3. Etching of MWCNTs for the impurity-free nanoelectrodes (Chapter 5). We show that the residual Ni catalyst on MWCNTs can be removed by acid etching; the tip removal and collapsing of tubes into pyramids enhances the stability of field emission from the tube arrays. The acid-etching process can be used to functionalize the MWCNTs, which was used to make our initial CNT-nanoelectrode glucose sensors. Finally, lessons learned trying to perform spectroscopic analysis of the functionalized MWCNTs were vital for designing our final devices. 4. Molecular junction design and electrochemical synthesis of biphenyl molecules on carbon microelectrodes for all-carbon molecular devices (Chapter 6). Utilizing the experience gained on the work done so far, our final device design is described. We demonstrate the capability of preparing patterned glassy carbon films to serve as the bottom electrode in the new geometry. However, the molecular switching behavior of biphenyl was not observed by scanning tunneling microscopy (STM), mercury drop or fabricated glassy carbon/biphenyl/MWCNT junctions. Either the density of these molecules is not optimum for effective integration of devices using MWCNTs as the nanoelectrodes, or an electroactive contaminant was reduced instead of the ionic biphenyl species. 5. Self-assembly of octadecanethiol (ODT) molecules on gold microelectrodes for functional molecular devices (Chapter 7). We have realized an effective scheme to produce Au/ODT/MWCNT junctions by spanning MWCNTs across ODT-functionalized microelectrodes. A percentage of the resulting junctions retain the expected character of an ODT monolayer. While the process is not yet optimized, our successful junctions show that molecular electronic devices can be fabricated using simple processes such as photolithography, self-assembled monolayers and dielectrophoresis.

Comparison of magnetic properties and petrography between dykes and lava flows from La Cienega, New Mexico and Thunder Bay area, Canada

Data on the evolution of geomagnetic paleointensity are crucial for understanding the geodynamo and Earth’s thermal history. Although basaltic flows are preferred for paleointensity experiments, quickly cooled mafic dykes have also been used. However, the paleointensity values obtained from the dykes are systematically lower than those from lava flows. This bias may originate from the difference in cooling histories and resultant magnetic mineralogies of extrusive and intrusive rocks. To explore this hypothesis, the magnetic mineralogy of two feeder dyke-lave flow systems, from Thunder Bay (Canada) and La Cienega (New-Mexico), has been studied using magnetic and microscopy methods. Within each system, the flow and dyke show different stages of deuteric oxidation of titanomagnetite, but the oxidation stages also differ between the two systems. It is concluded that the tested hypothesis is viable, but the relationships between the magnetic and mineralogical properties of flows and dykes are complex and need a further investigation.

Automation of shear-wave splitting parameter determination of local earthquakes at Yellowstone : application as indicator of crustal stress and temporal variation

Shear-wave splitting can be a useful technique for determining crustal stress fields in volcanic settings and temporal variations associated with activity. Splitting parameters were determined for a subset of local earthquakes recorded from 2000-2010 at Yellowstone. Analysis was automated using an unsupervised cluster analysis technique to determine optimum splitting parameters from 270 analysis windows for each event. Six stations clearly exhibit preferential fast polarization values sub-orthogonal to the direction of minimum horizontal compression. Yellowstone deformation results in a local crustal stress field differing from the regional field dominated by NE-SW extension, and fast directions reflect this difference rotating around the caldera maintaining perpendicularity to the rim. One station exhibits temporal variations concordant with identified periods of caldera subsidence and uplift. From splitting measurements, we calculated a crustal anisotropy of ~17-23% and crack density ~0.12-0.17 possibly resulting from stress-aligned fluid filled microcracks in the upper crust and an active hydrothermal system.

Numerical and granulometric approaches to geophysical granular flows

Volcanic ash clouds can be fed by an upward-directed eruption column (Plinian column) or by elutriation from extensive pyroclastic-flows (coignimbrite cloud). For large-scale eruptions, there is considerable uncertainty about which mechanism is dominant. Here we analyze in a novel way a comprehensive grainsize database for pyroclastic deposits. We demonstrate that the Mount Pinatubo climactic eruption deposits were substantially derived from coignimbrite clouds, and not only by a Plinian cloud as generally thought. Coignimbrite ash-fall deposits are much richer in breathable <10 m ash (5–25 wt%) than pure Plinian ash at most distances from the source volcano. We also show that coignimbrite ash clouds, as at Pinatubo, are expected to be more water rich than Plinian clouds, leading to removal of more HCl prior to stratospheric injection, thereby reducing their atmospheric impact.

Charge and spin transport in nanoscale junction from first principles

Recently nanoscale junctions consisting of 0-D nanostructures (single molecule) or 1-D nanostructures (semiconducting nanowire) sandwiched between two metal electrodes are successfully fabricated and characterized. What lacks in the recent developments is the understanding of the mechanism behind the observed phenomena at the level of atoms and electrons. For example, the origin of observed switching effect in a semiconducting nanowire due to the influence of an external gate bias is not yet understood at the electronic structure level. On the same context, different experimental groups have reported different signs in tunneling magneto-resistance for the same organic spin valve structure, which has baffled researchers working in this field. In this thesis, we present the answers to some of these subtle questions by investigating the charge and spin transport in different nanoscale junctions. A parameter-free, single particle Green’s function approach in conjunction with a posteriori density functional theory (DFT) involving a hybrid orbital dependent functional is used to calculate the tunneling current in the coherent transport limit. The effect of spin polarization is explicitly incorporated to investigate spin transport in a nanoscale junction. Through the electron transport studies in PbS nanowire junction, a new orbital controlled mechanism behind the switching of the current is proposed. It can explain the switching behavior, not only in PbS nanowire, but in other lead-chalcogenide nanowires as well. Beside this, the electronic structure properties of this nanowire are studied using periodic DFT. The quantum confinement effect was investigated by calculating the bandgap of PbS nanowires with different diameters. Subsequently, we explain an observed semiconducting to metallic phase transition of this nanowire by calculating the bandgap of the nanowire under uniform radial strain. The compressive radial strain on the nanowire was found to be responsible for the metallic to semiconducting phase transition. Apart from studying one dimensional nanostructure, we also present transport properties in zero dimensional single molecular junctions. We proposed a new codoping approach in a single molecular carborane junction, where a cation and an anion are simultaneously doped to find the role of a single atom in the device. The main purpose was to build a molecular junction where a single atom can dictate the flow of electrons in a circuit. Recent observations of both positive and negative sign in tunneling magnetoresistance (TMR) the using same organic spin-valve structure hasmystified researchers. From our spin dependent transport studies in a prototypical organic molecular tunneling device, we found that a 3% change in metal-molecule interfacial distance can alter the sign of TMR. Changing the interfacial distance by 3%, the number of participating eigenstates as well as their orbital characteristic changes for anti-parallel configuration of the magnetization at the two electrodes, leading to the sign reversal of the TMR. Apart from this, the magnetic proximity effect under applied bias is investigated quantitatively, which can be used to understand the observed unexpectedmagnetismin carbon basedmaterials when they are in close proximity with magnetic substrates.

Fixed block configuration GDDs with block size 6 and (3, r)-regular graphs

Chapter 1 is used to introduce the basic tools and mechanics used within this thesis. Most of the definitions used in the thesis will be defined, and we provide a basic survey of topics in graph theory and design theory pertinent to the topics studied in this thesis. In Chapter 2, we are concerned with the study of fixed block configuration group divisible designs, GDD(n; m; k; λ1; λ2). We study those GDDs in which each block has configuration (s; t), that is, GDDs in which each block has exactly s points from one of the two groups and t points from the other. Chapter 2 begins with an overview of previous results and constructions for small group size and block sizes 3, 4 and 5. Chapter 2 is largely devoted to presenting constructions and results about GDDs with two groups and block size 6. We show the necessary conditions are sufficient for the existence of GDD(n, 2, 6; λ1, λ2) with fixed block configuration (3; 3). For configuration (1; 5), we give minimal or nearminimal index constructions for all group sizes n ≥ 5 except n = 10, 15, 160, or 190. For configuration (2, 4), we provide constructions for several families ofGDD(n, 2, 6; λ1, λ2)s. Chapter 3 addresses characterizing (3, r)-regular graphs. We begin with providing previous results on the well studied class of (2, r)-regular graphs and some results on the structure of large (t; r)-regular graphs. In Chapter 3, we completely characterize all (3, 1)-regular and (3, 2)-regular graphs, as well has sharpen existing bounds on the order of large (3, r)- regular graphs of a certain form for r ≥ 3. Finally, the appendix gives computational data resulting from Sage and C programs used to generate (3, 3)-regular graphs on less than 10 vertices.

Game-theoretic view on intermediated exchange

Intermediaries permeate modern economic exchange. Most classical models on intermediated exchange are driven by information asymmetry and inventory management. These two factors are of reduced significance in modern economies. This makes it necessary to develop models that correspond more closely to modern financial marketplaces. The goal of this dissertation is to propose and examine such models in a game theoretical context. The proposed models are driven by asymmetries in the goals of different market participants. Hedging pressure as one of the most critical aspects in the behavior of commercial entities plays a crucial role. The first market model shows that no equilibrium solution can exist in a market consisting of a commercial buyer, a commercial seller and a non-commercial intermediary. This indicates a clear economic need for non-commercial trading intermediaries: a direct trade from seller to buyer does not result in an equilibrium solution. The second market model has two distinct intermediaries between buyer and seller: a spread trader/market maker and a risk-neutral intermediary. In this model a unique, natural equilibrium solution is identified in which the supply-demand surplus is traded by the risk-neutral intermediary, whilst the market maker trades the remainder from seller to buyer. Since the market maker’s payoff for trading at the identified equilibrium price is zero, this second model does not provide any motivation for the market maker to enter the market. The third market model introduces an explicit transaction fee that enables the market maker to secure a positive payoff. Under certain assumptions on this transaction fee the equilibrium solution of the previous model applies and now also provides a financial motivation for the market maker to enter the market. If the transaction fee violates an upper bound that depends on supply, demand and riskaversity of buyer and seller, the market will be in disequilibrium.

Comparing satellite and ground-based observations of paroxysmal degassing events at Etna volcano, Italy

Mount Etna, Italy, is one of the most active volcanoes in the world, and is also regarded as one of the strongest volcanic sources of sulfur dioxide (SO2) emissions to the atmosphere. Since October 2004, an automated ultraviolet (UV) spectrometer network (FLAME) has provided ground-based SO2 measurements with high temporal resolution, providing an opportunity to validate satellite SO2 measurements at Etna. The Ozone Monitoring Instrument (OMI) on the NASA Aura satellite, which makes global daily measurements of trace gases in the atmosphere, was used to compare SO2 amount released by the volcano during paroxysmal lava-fountaining events from 2004 to present. We present the first comparison between SO2 emission rates and SO2 burdens obtained by the OMI transect technique and OMI Normalized Cloud-Mass (NCM) technique and the ground-based FLAME Mini-DOAS measurements. In spite of a good data set from the FLAME network, finding coincident OMI and FLAME measurements proved challenging and only one paroxysmal event provided a good validation for OMI. Another goal of this work was to assess the efficacy of the FLAME network in capturing paroxysmal SO2 emissions from Etna, given that the FLAME network is only operational during daylight hours and some paroxysms occur at night. OMI measurements are advantageous since SO2 emissions from nighttime paroxysms can often be quantified on the following day, providing improved constraints on Etna’s SO2 budget.

Analyzing the effects of vertical acceleration and deceleration in multi-phase, debris-flow sediment-column experiments / by Christine Marie Williams.

Water-saturated debris flows are among some of the most destructive mass movements. Their complex nature presents a challenge for quantitative description and modeling. In order to improve understanding of the dynamics of these flows, it is important to seek a simplified dynamic system underlying their behavior. Models currently in use to describe the motion of debris flows employ depth-averaged equations of motion, typically assuming negligible effects from vertical acceleration. However, in many cases debris flows experience significant vertical acceleration as they move across irregular surfaces, and it has been proposed that friction associated with vertical forces and liquefaction merit inclusion in any comprehensive mechanical model. The intent of this work is to determine the effect of vertical acceleration through a series of laboratory experiments designed to simulate debris flows, testing a recent model for debris flows experimentally. In the experiments, a mass of water-saturated sediment is released suddenly from a holding container, and parameters including rate of collapse, pore-fluid pressure, and bed load are monitored. Experiments are simplified to axial geometry so that variables act solely in the vertical dimension. Steady state equations to infer motion of the moving sediment mass are not sufficient to model accurately the independent solid and fluid constituents in these experiments. The model developed in this work more accurately predicts the bed-normal stress of a saturated sediment mass in motion and illustrates the importance of acceleration and deceleration.

Controlling electronic and magnetic properties of ultra narrow multilayered nanowires

Interest in the study of magnetic/non-magnetic multilayered structures took a giant leap since Grünberg and his group established that the interlayer exchange coupling (IEC) is a function of the non-magnetic spacer width. This interest was further fuelled by the discovery of the phenomenal Giant Magnetoresistance (GMR) effect. In fact, in 2007 Albert Fert and Peter Grünberg were awarded the Nobel Prize in Physics for their contribution to the discovery of GMR. GMR is the key property that is being used in the read-head of the present day computer hard drive as it requires a high sensitivity in the detection of magnetic field. The recent increase in demand for device miniaturization encouraged researchers to look for GMR in nanoscale multilayered structures. In this context, one dimensional(1-D) multilayerd nanowire structure has shown tremendous promise as a viable candidate for ultra sensitive read head sensors. In fact, the phenomenal giant magnetoresistance(GMR) effect, which is the novel feature of the currently used multilayered thin film, has already been observed in multilayered nanowire systems at ambient temperature. Geometrical confinement of the supper lattice along the 2-dimensions (2-D) to construct the 1-D multilayered nanowire prohibits the minimization of magnetic interaction- offering a rich variety of magnetic properties in nanowire that can be exploited for novel functionality. In addition, introduction of non-magnetic spacer between the magnetic layers presents additional advantage in controlling magnetic properties via tuning the interlayer magnetic interaction. Despite of a large volume of theoretical works devoted towards the understanding of GMR and IEC in super lattice structures, limited theoretical calculations are reported in 1-D multilayered systems. Thus to gauge their potential application in new generation magneto-electronic devices, in this thesis, I have discussed the usage of first principles density functional theory (DFT) in predicting the equilibrium structure, stability as well as electronic and magnetic properties of one dimensional multilayered nanowires. Particularly, I have focused on the electronic and magnetic properties of Fe/Pt multilayered nanowire structures and the role of non-magnetic Pt spacer in modulating the magnetic properties of the wire. It is found that the average magnetic moment per atom in the nanowire increases monotonically with an ~1/(N(Fe)) dependance, where N(Fe) is the number of iron layers in the nanowire. A simple model based upon the interfacial structure is given to explain the 1/(N(Fe)) trend in magnetic moment obtained from the first principle calculations. A new mechanism, based upon spin flip with in the layer and multistep electron transfer between the layers, is proposed to elucidate the enhancement of magnetic moment of Iron atom at the Platinum interface. The calculated IEC in the Fe/Pt multilayered nanowire is found to switch sign as the width of the non-magnetic spacer varies. The competition among short and long range direct exchange and the super exchange has been found to play a key role for the non-monotonous sign in IEC depending upon the width of the Platinum spacer layer. The calculated magnetoresistance from Julliere's model also exhibit similar switching behavior as that of IEC. The universality of the behavior of exchange coupling has also been looked into by introducing different non-magnetic spacers like Palladium, Copper, Silver, and Gold in between magnetic Iron layers. The nature of hybridization between Fe and other non-magnetic spacer is found to dictate the inter layer magnetic interaction. For example, in Fe/Pd nanowire the d-p hybridization in two spacer layer case favors anti-ferromagnetic (AFM) configuration over ferromagnetic (FM) configuration. However, the hybridization between half-filled Fe(d) and filled Cu(p) state in Fe/Cu nanowire favors FM coupling in the 2-spacer system.

Geophysical and geological analysis of the Collingwood Member of the Trenton Formation

The Collingwood Member is a mid to late Ordovician self-sourced reservoir deposited across the northern Michigan Basin and parts of Ontario, Canada. Although it had been previously studied in Canada, there has been relatively little data available from the Michigan subsurface. Recent commercial interest in the Collingwood has resulted in the drilling and production of several wells in the state of Michigan. An analysis of core samples, measured laboratory data, and petrophysical logs has yielded both a quantitative and qualitative understanding of the formation in the Michigan Basin. The Collingwood is a low permeability and low porosity carbonate package that is very high in organic content. It is composed primarily of a uniformly fine grained carbonate matrix with lesser amounts of kerogen, silica, and clays. The kerogen content of the Collingwood is finely dispersed in the clay and carbonate mineral phases. Geochemical and production data show that both oil and gas phases are present based on regional thermal maturity. The deposit is richest in the north-central part of the basin with thickest deposition and highest organic content. The Collingwood is a fairly thin deposit and vertical fractures may very easily extend into the surrounding formations. Completion and treatment techniques should be designed around these parameters to enhance production.

Simulating dioxane transport in a heterogeneous glacial aquifer system (Washtenaw County, Michigan) using publicly available models and data

The primary challenge in groundwater and contaminant transport modeling is obtaining the data needed for constructing, calibrating and testing the models. Large amounts of data are necessary for describing the hydrostratigraphy in areas with complex geology. Increasingly states are making spatial data available that can be used for input to groundwater flow models. The appropriateness of this data for large-scale flow systems has not been tested. This study focuses on modeling a plume of 1,4-dioxane in a heterogeneous aquifer system in Scio Township, Washtenaw County, Michigan. The analysis consisted of: (1) characterization of hydrogeology of the area and construction of a conceptual model based on publicly available spatial data, (2) development and calibration of a regional flow model for the site, (3) conversion of the regional model to a more highly resolved local model, (4) simulation of the dioxane plume, and (5) evaluation of the model's ability to simulate field data and estimation of the possible dioxane sources and subsequent migration until maximum concentrations are at or below the Michigan Department of Environmental Quality's residential cleanup standard for groundwater (85 ppb). MODFLOW-2000 and MT3D programs were utilized to simulate the groundwater flow and the development and movement of the 1, 4-dioxane plume, respectively. MODFLOW simulates transient groundwater flow in a quasi-3-dimensional sense, subject to a variety of boundary conditions that can simulate recharge, pumping, and surface-/groundwater interactions. MT3D simulates solute advection with groundwater flow (using the flow solution from MODFLOW), dispersion, source/sink mixing, and chemical reaction of contaminants. This modeling approach was successful at simulating the groundwater flows by calibrating recharge and hydraulic conductivities. The plume transport was adequately simulated using literature dispersivity and sorption coefficients, although the plume geometries were not well constrained.

Sulfoxides as an intramolecular sulfenylating agent for indoles and diverse applications of the sulfide-sulfoxide redox cycle in organic chemistry

This dissertation involves study of various aspects of sulfoxide chemistry. Specifically designed t-butyl and propanenitrile sulfoxides tethered to indole-2-carboxamide were used as a source of intramolecular sulfenylating agents to synthesize novel indolo[3,2-b]-1-5-benzothiazepinones which are structurally analogous to the other biologically active benzothiazepinones. This study reveals that the intramolecular cyclization of sulfoxide follows an electrophilic sulfenylation (Sulfoxide Electrophilic Sulfenylation, SES) reaction pathway. Evidence of the absence of sulfenic acid as a transient reactive intermediate in such intramolecular cyclization is also provided. In another study, sulfoxide was used as a “protecting group” of thioether to synthesize 8-membered, indole substituted, thiazocine-2-acetic acid derivative via Ring Closing Metathesis (RCM). Protection (oxidation) of inert (to RCM) sulfide to sulfoxide followed by RCM produced cyclized product in good yields. Deprotection (reduction) of sulfoxide was achieved using Lawessons Reagent (L.R.). Application of the sulfide-sulfoxide redox cycle to solve the existing difficulties in using RCM methodology to thioethers is illustrated. A new design of a “molecular brake”, based on the sulfide-sulfoxide redox cycle is described. N-Ar rotation in simple isoindolines is controlled by the oxidation state of the proximate sulfur atom. Sulfide [S(II)] shows “free” [brake OFF] N-Ar rotation whereas sulfoxide displayed hindered [brake ON] N-Ar rotation. The semi-empirical molecular orbital (PM3) calculations revealed concerted pyramidalization of amidic nitrogen with N-Ar rotation.

Planar magneto-photonic and gradient-photonic structures : crystals and metamaterials.

In the field of photonics, two new types of material structures, photonic crystals and metamaterials, are presently of great interest. Both are studied in the present work, which focus on planar magnetic materials in the former and planar gradient metamaterials in the latter. These planar periodic structures are easy to handle and integrate into optical systems. The applications are promising field for future optical telecommunication systems and give rise to new optical, microwave and radio technologies. The photonic crystal part emphasizes the utilization of magnetic material based photonic crystals due to its remarkable magneto-optical characteristics. Bandgaps tuning by magnetic field in bismuth-gadolinium-substituted lutetium iron garnet (Bi0.8 Gd0.2 Lu2.0 Fe5 O12) based one- dimensional photonic crystals are investigated and demonstrated in this work. Magnetic optical switches are fabricated and tested. Waveguide formulation for band structure in magneto photonic crystals is developed. We also for the first time demonstrate and test two- dimensional magneto photonic crystals optical. We observe multi-stopbands in two- dimensional photonic waveguide system and study the origin of multi-stopbands. The second part focus on studying photonic metamaterials and planar gradient photonic metamaterial design. We systematically study the effects of varying the geometry of the fishnet unit cell on the refractive index in optical frequency. It is the first time to design and demonstrate the planar gradient structure in the high optical frequency. Optical beam bending using planar gradient photonic metamaterials is observed. The technologies needed for the fabrication of the planar gradient photonic metamaterials are investigated. Beam steering devices, shifter, gradient optical lenses and etc. can be derived from this design.

Edge coloring BIBDS and constructing MOELRs

Chapter 1 is used to introduce the basic tools and mechanics used within this thesis. Some historical uses and background are touched upon as well. The majority of the definitions are contained within this chapter as well. In Chapter 2 we consider the question whether one can decompose λ copies of monochromatic Kv into copies of Kk such that each copy of the Kk contains at most one edge from each Kv. This is called a proper edge coloring (Hurd, Sarvate, [29]). The majority of the content in this section is a wide variety of examples to explain the constructions used in Chapters 3 and 4. In Chapters 3 and 4 we investigate how to properly color BIBD(v, k, λ) for k = 4, and 5. Not only will there be direct constructions of relatively small BIBDs, we also prove some generalized constructions used within. In Chapter 5 we talk about an alternate solution to Chapters 3 and 4. A purely graph theoretical solution using matchings, augmenting paths, and theorems about the edgechromatic number is used to develop a theorem that than covers all possible cases. We also discuss how this method performed compared to the methods in Chapters 3 and 4. In Chapter 6, we switch topics to Latin rectangles that have the same number of symbols and an equivalent sized matrix to Latin squares. Suppose ab = n2. We define an equitable Latin rectangle as an a × b matrix on a set of n symbols where each symbol appears either [b/n] or [b/n] times in each row of the matrix and either [a/n] or [a/n] times in each column of the matrix. Two equitable Latin rectangles are orthogonal in the usual way. Denote a set of ka × b mutually orthogonal equitable Latin rectangles as a k–MOELR(a, b; n). We show that there exists a k–MOELR(a, b; n) for all a, b, n where k is at least 3 with some exceptions.