Outer crust of a cold non-accreting magnetar

The outer-crust structure and composition of a cold, non-accreting magnetar are studied. We model the outer crust to be made of fully equilibrated matter where ionized nuclei form a Coulomb crystal embedded in an electron gas. The main effects of the strong magnetic field are those of quantizing the electron motion in Landau levels and of modifying the nuclear single-particle levels producing, on average, an increased binding of nucleons in nuclei present in the Coulomb lattice. The effect of a homogeneous and constant magnetic field on nuclear masses has been predicted by using a covariant density functional in which induced currents and axial deformation due to the presence of a magnetic field that breaks time-reversal symmetry have been included self-consistently in the nucleon and meson equations of motion. Although not yet observed, for Ba 1016 G both effects contribute to produce different compositions - odd-mass nuclei are frequently predicted - and to increase the neutron-drip pressure as compared to a typical neutron star. Specifically, in such a regime, the magnetic-field effects on nuclei favor the appearance of heavier nuclei at low pressures. As B increases, such heavier nuclei are also preferred up to larger pressures. For the most extreme magnetic field considered, B=1018 G, and for the models studied, almost the whole outer crust is made of 4092Zr52.

Solar magnetic carpet I: simulation of synthetic magnetograms

This paper describes a new 2D model for the photospheric evolution of the magnetic carpet. It is the first in a series of papers working towards constructing a realistic 3D non-potential model for the interaction of small-scale solar magnetic fields. In the model, the basic evolution of the magnetic elements is governed by a supergranular flow profile. In addition, magnetic elements may evolve through the processes of emergence, cancellation, coalescence and fragmentation. Model parameters for the emergence of bipoles are based upon the results of observational studies. Using this model, several simulations are considered, where the range of flux with which bipoles may emerge is varied. In all cases the model quickly reaches a steady state where the rates of emergence and cancellation balance. Analysis of the resulting magnetic field shows that we reproduce observed quantities such as the flux distribution, mean field, cancellation rates, photospheric recycle time and a magnetic network. As expected, the simulation matches observations more closely when a larger, and consequently more realistic, range of emerging flux values is allowed (4×1016 - 1019 Mx). The model best reproduces the current observed properties of the magnetic carpet when we take the minimum absolute flux for emerging bipoles to be 4×1016 Mx. In future, this 2D model will be used as an evolving photospheric boundary condition for 3D non-potential modeling.

Solar magnetic carpet II: coronal interactions of small-scale magnetic fields

This paper is the second in a series of studies working towards constructing a realistic, evolving, non-potential coronal model for the solar magnetic carpet. In the present study, the interaction of two magnetic elements is considered. Our objectives are to study magnetic energy build-up, storage and dissipation as a result of emergence, cancellation, and flyby of these magnetic elements. In the future these interactions will be the basic building blocks of more complicated simulations involving hundreds of elements. Each interaction is simulated in the presence of an overlying uniform magnetic field, which lies at various orientations with respect to the evolving magnetic elements. For these three small-scale interactions, the free energy stored in the field at the end of the simulation ranges from 0.2 – 2.1×1026 ergs, whilst the total energy dissipated ranges from 1.3 – 6.3×1026 ergs. For all cases, a stronger overlying field results in higher energy storage and dissipation. For the cancellation and emergence simulations, motion perpendicular to the overlying field results in the highest values. For the flyby simulations, motion parallel to the overlying field gives the highest values. In all cases, the free energy built up is sufficient to explain small-scale phenomena such as X-ray bright points or nanoflares. In addition, if scaled for the correct number of magnetic elements for the volume considered, the energy continually dissipated provides a significant fraction of the quiet Sun coronal heating budget.

Solar magnetic carpet III: coronal modelling of synthetic magnetograms

This article is the third in a series working towards the construction of a realistic, evolving, non-linear force-free coronal-field model for the solar magnetic carpet. Here, we present preliminary results of 3D time-dependent simulations of the small-scale coronal field of the magnetic carpet. Four simulations are considered, each with the same evolving photospheric boundary condition: a 48-hour time series of synthetic magnetograms produced from the model of Meyer et al. ( Solar Phys. 272, 29, 2011). Three simulations include a uniform, overlying coronal magnetic field of differing strength, the fourth simulation includes no overlying field. The build-up, storage, and dissipation of magnetic energy within the simulations is studied. In particular, we study their dependence upon the evolution of the photospheric magnetic field and the strength of the overlying coronal field. We also consider where energy is stored and dissipated within the coronal field. The free magnetic energy built up is found to be more than sufficient to power small-scale, transient phenomena such as nanoflares and X-ray bright points, with the bulk of the free energy found to be stored low down, between 0.5 - 0.8 Mm. The energy dissipated is currently found to be too small to account for the heating of the entire quiet-Sun corona. However, the form and location of energy-dissipation regions qualitatively agree with what is observed on small scales on the Sun. Future MHD modelling using the same synthetic magnetograms may lead to a higher energy release.

The storage and dissipation of magnetic energy in the quiet sun corona determined from SDO/HMI magnetograms

In recent years, higher cadence, higher resolution observations have revealed the quiet-Sun photosphere to be complex and rapidly evolving. Since magnetic fields anchored in the photosphere extend up into the solar corona, it is expected that the small-scale coronal magnetic field exhibits similar complexity. For the first time, the quiet-Sun coronal magnetic field is continuously evolved through a series of non-potential, quasi-static equilibria, deduced from magnetograms observed by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory, where the photospheric boundary condition which drives the coronal evolution exactly reproduces the observed magnetograms. The build-up, storage, and dissipation of magnetic energy within the simulations is studied. We find that the free magnetic energy built up and stored within the field is sufficient to explain small-scale, impulsive events such as nanoflares. On comparing with coronal images of the same region, the energy storage and dissipation visually reproduces many of the observed features. The results indicate that the complex small-scale magnetic evolution of a large number of magnetic features is a key element in explaining the nature of the solar corona.

A finite-volume shallow layer method for the MHD instabilities in an aluminium production cell

An electrolytic cell for Aluminium production contains molten metal and molten electrolyte, which are subject to high dc-currents and magnetic fields. Lorentz forces arising from the cross product of current and magnetic field may amplify natural gravity waves at the interface between the two fluids, leading to short circuits in extreme cases. The external magnetic field and current distribution in the production cell is computed through a detailed finite element analysis at Torino Polytechnic. The results are then used to compute the magnetohydrodynamic and thermal effects in the aluminium/electrolyte bath. Each cell has lateral dimensions of 6m x 2m, whilst the bath depth is only 30cm. the electrically resistive electrolyte path, which is critical in the operation of the cell, has layer depth of only a few centimetres below each carbon anode. Because the shallow dimensions of the liquid layer a finite-volume shallow-layer technique has been used at Greenwich to compute the resulting flow-field and interface perturbations. The information obtained from this method, i.e. depth averaged velocities and aluminium/electrolyte interface position is then embedded in the three-dimensional finite volume code PHYSICA and will be used to compute the heat transfer and phase change in the cell.

Ionisation and discharge in cloud-forming atmospheres of brown dwarfs and extrasolar planets

Brown dwarfs and giant gas extrasolar planets have cold atmospheres with rich chemical compositions from which mineral cloud particles form. Their properties, like particle sizes and material composition, vary with height, and the mineral cloud particles are charged due to triboelectric processes in such dynamic atmospheres. The dynamics of the atmospheric gas is driven by the irradiating host star and/or by the rotation of the objects that changes during its lifetime. Thermal gas ionisation in these ultra-cool but dense atmospheres allows electrostatic interactions and magnetic coupling of a substantial atmosphere volume. Combined with a strong magnetic field , a chromosphere and aurorae might form as suggested by radio and x-ray observations of brown dwarfs. Non-equilibrium processes like cosmic ray ionisation and discharge processes in clouds will increase the local pool of free electrons in the gas. Cosmic rays and lighting discharges also alter the composition of the local atmospheric gas such that tracer molecules might be identified. Cosmic rays affect the atmosphere through air showers in a certain volume which was modelled with a 3D Monte Carlo radiative transfer code to be able to visualise their spacial extent. Given a certain degree of thermal ionisation of the atmospheric gas, we suggest that electron attachment to charge mineral cloud particles is too inefficient to cause an electrostatic disruption of the cloud particles. Cloud particles will therefore not be destroyed by Coulomb explosion for the local temperature in the collisional dominated brown dwarf and giant gas planet atmospheres. However, the cloud particles are destroyed electrostatically in regions with strong gas ionisation. The potential size of such cloud holes would, however, be too small and might occur too far inside the cloud to mimic the effect of, e.g. magnetic field induced star spots.

L'orientation et la propriété de mémoire de forme des polymères cristallins liquides à chaînes latérales covalents et supramoléculaires

In many studies of the side-chain liquid crystalline polymers (SCLCPs) bearing azobenzene mesogens as pendant groups, obtaining the orientation of azobenzene mesogens at a macroscopic scale as well as its control is important, because it impacts many properties related to the cooperative motion characteristic of liquid crystals and the trans-cis photoisomerization of the azobenzene molecules. Various means can be used to align the mesogens in the polymers, including rubbed surface, mechanical stretching or shearing, and electric or magnetic field. In the case of azobenzene-containing SCLCPs, another method consists in using linearly polarized light (LPL) to induce orientation of azobenzene mesogens perpendicular to the polarization direction of the excitation light, and such photoinduced orientation has been the subject of numerous studies. In the first study realized in this thesis (Chapter 1), we carried out the first systematic investigation on the interplay of the mechanically and optically induced orientation of azobenzene mesogens as well as the effect of thermal annealing in a SCLCP and a diblock copolymer comprising two SCLCPs bearing azobenzene and biphenyl mesogens, respectively. Using a supporting-film approach previously developed by our group, a given polymer film can be first stretched in either the nematic or smectic phase to yield orientation of azobenzene mesogens either parallel or perpendicular to the strain direction, then exposed to unpolarized UV light to erase the mechanically induced orientation upon the trans–cis isomerization, followed by linearly polarized visible light for photoinduced reorientation as a result of the cis–trans backisomerization, and finally heated to different LC phases for thermal annealing. Using infrared dichroism to monitor the change in orientation degree, the results of this study have unveiled complex and different orientational behavior and coupling effects for the homopolymer of poly{6-[4-(4-methoxyphenylazo)phenoxy]hexyl methacrylate} (PAzMA) and the diblock copolymer of PAzMA-block- poly{6-[4-(4-cyanophenyl) phenoxy]hexyl methacrylate} (PAzMA-PBiPh). Most notably for the homopolymer, the stretching-induced orientation exerts no memory effect on the photoinduced reorientation, the direction of which is determined by the polarization of the visible light regardless of the mechanically induced orientation direction in the stretched film. Moreover, subsequent thermal annealing in the nematic phase leads to parallel orientation independently of the initial mechanically or photoinduced orientation direction. By contrast, the diblock copolymer displays a strong orientation memory effect. Regardless of the condition used, either for photoinduced reorientation or thermal annealing in the liquid crystalline phase, only the initial stretching-induced perpendicular orientation of azobenzene mesogens can be recovered. The reported findings provide new insight into the different orientation mechanisms, and help understand the important issue of orientation induction and control in azobenzene-containing SCLCPs. The second study presented in this thesis (Chapter 2) deals with supramolecular side-chain liquid crystalline polymers (S-SCLCPs), in which side-group mesogens are linked to the chain backbone through non-covalent interactions such as hydrogen bonding. Little is known about the mechanically induced orientation of mesogens in S-SCLCPs. In contrast to covalent SCLCPs, free-standing, solution-cast thin films of a S-SCLCP, built up with 4-(4’-heptylphenyl) azophenol (7PAP) H-bonded to poly(4-vinyl pyridine) (P4VP), display excellent stretchability. Taking advantage of this finding, we investigated the stretching-induced orientation and the viscoelastic behavior of this S-SCLCP, and the results revealed major differences between supramolecular and covalent SCLCPs. For covalent SCLCPs, the strong coupling between chain backbone and side-group mesogens means that the two constituents can mutually influence each other; the lack of chain entanglements is a manifestation of this coupling effect, which accounts for the difficulty in obtaining freestanding and mechanically stretchable films. Upon elongation of a covalent SCLCP film cast on a supporting film, the mechanical force acts on the coupled polymer backbone and mesogenic side groups, and the latter orients cooperatively and efficiently (high orientation degree), which, in turn, imposes an anisotropic conformation of the chain backbone (low orientation degree). In the case of the S-SCLCP of P4VP-7PAP, the coupling between the side-group mesogens and the chain backbone is much weakened owing to the dynamic dissociation/association of the H-bonds linking the two constituents. The consequence of this decoupling is readily observable from the viscoelastic behavior. The average molecular weight between entanglements is basically unchanged in both the smectic and isotropic phase, and is similar to non-liquid crystalline samples. As a result, the S-SCLCP can easily form freestanding and stretchable films. Furthermore, the stretching induced orientation behavior of P4VP-7PAP is totally different. Stretching in the smectic phase results in a very low degree of orientation of the side-group mesogens even at a large strain (500%), while the orientation of the main chain backbone develops steadily with increasing the strain, much the same way as amorphous polymers. The results imply that upon stretching, the mechanical force is mostly coupled to the polymer backbone and leads to its orientation, while the main chain orientation exerts little effect on orienting the H-bonded mesogenic side groups. This surprising finding is explained by the likelihood that during stretching in the smectic phase (at relatively higher temperatures) the dynamic dissociation of the H-bonds allow the side-group mesogens to be decoupled from the chain backbone and relax quickly. In the third project (Chapter 3), we investigated the shape memory properties of a S-SCLCP prepared by tethering two azobenzene mesogens, namely, 7PAP and 4-(4'-ethoxyphenyl) azophenol (2OPAP), to P4VP through H-bonding. The results revealed that, despite the dynamic nature of the linking H-bonds, the supramolecular SCLCP behaves similarly to covalent SCLCP by exhibiting a two-stage thermally triggered shape recovery process governed by both the glass transition and the LC-isotropic phase transition. The ability for the supramolecular SCLCP to store part of the strain energy above T[subscript g] in the LC phase enables the triple-shape memory property. Moreover, thanks to the azobenzene mesogens used, which can undergo trans-cis photoisomerization, exposure the supramolecular SCLCP to UV light can also trigger the shape recovery process, thus enabling the remote activation and the spatiotemporal control of the shape memory. By measuring the generated contractile force and its removal upon turning on and off the UV light, respectively, on an elongated film under constant strain, it seems that the optically triggered shape recovery stems from a combination of a photothermal effect and an effect of photoplasticization or of an order-disorder phase transition resulting from the trans-cis photoisomerization of azobenzene mesogens.

Disentangling the magnetic force noise contribution in LISA Pathfinder

Magnetically-induced forces on the inertial masses on-board LISA Path finder are expected to be one of the dominant contributions to the mission noise budget, accounting for up to 40%. The origin of this disturbance is the coupling of the residual magnetization and susceptibility of the test masses with the environmental magnetic field. In order to fully understand this important part of the noise model, a set of coils and magnetometers are integrated as a part of the diagnostics subsystem. During operations a sequence of magnetic excitations will be applied to precisely determine the coupling of the magnetic environment to the test mass displacement using the on-board magnetometers. Since no direct measurement of the magnetic field in the test mass position will be available, an extrapolation of the magnetic measurements to the test mass position will be carried out as a part of the data analysis activities. In this paper we show the first results on the magnetic experiments during an end-to-end LISA Path finder simulation, and we describe the methods under development to map the magnetic field on-board.

Engineering thin films of magnetic alloys and semiconductor oxides at the nanoscale

The thesis aims to exploit properties of thin films for applications such as spintronics, UV detection and gas sensing. Nanoscale thin films devices have myriad advantages and compatibility with Si-based integrated circuits processes. Two distinct classes of material systems are investigated, namely ferromagnetic thin films and semiconductor oxides. To aid the designing of devices, the surface properties of the thin films were investigated by using electron and photon characterization techniques including Auger electron spectroscopy (AES), X-ray photoelectron spectroscopy (XPS), grazing incidence X-ray diffraction (GIXRD), and energy-dispersive X-ray spectroscopy (EDS). These are complemented by nanometer resolved local proximal probes such as atomic force microscopy (AFM), magnetic force microscopy (MFM), electric force microscopy (EFM), and scanning tunneling microscopy to elucidate the interplay between stoichiometry, morphology, chemical states, crystallization, magnetism, optical transparency, and electronic properties. Specifically, I studied the effect of annealing on the surface stoichiometry of the CoFeB/Cu system by in-situ AES and discovered that magnetic nanoparticles with controllable areal density can be produced. This is a good alternative for producing nanoparticles using a maskless process. Additionally, I studied the behavior of magnetic domain walls of the low coercivity alloy CoFeB patterned nanowires. MFM measurement with the in-plane magnetic field showed that, compared to their permalloy counterparts, CoFeB nanowires require a much smaller magnetization switching field , making them promising for low-power-consumption domain wall motion based devices. With oxides, I studied CuO nanoparticles on SnO2 based UV photodetectors (PDs), and discovered that they promote the responsivity by facilitating charge transfer with the formed nanoheterojunctions. I also demonstrated UV PDs with spectrally tunable photoresponse with the bandgap engineered ZnMgO. The bandgap of the alloyed ZnMgO thin films was tailored by varying the Mg contents and AES was demonstrated as a surface scientific approach to assess the alloying of ZnMgO. With gas sensors, I discovered the rf-sputtered anatase-TiO2 thin films for a selective and sensitive NO2 detection at room temperature, under UV illumination. The implementation of UV enhances the responsivity, response and recovery rate of the TiO2 sensor towards NO2 significantly. Evident from the high resolution XPS and AFM studies, the surface contamination and morphology of the thin films degrade the gas sensing response. I also demonstrated that surface additive metal nanoparticles on thin films can improve the response and the selectivity of oxide based sensors. I employed nanometer-scale scanning probe microscopy to study a novel gas senor scheme consisting of gallium nitride (GaN) nanowires with functionalizing oxides layer. The results suggested that AFM together with EFM is capable of discriminating low-conductive materials at the nanoscale, providing a nondestructive method to quantitatively relate sensing response to the surface morphology.

Magnetoresistance in a High Mobility Two-Dimensional Electron System as a Function of Sample Geometry

In a high mobility two-dimensional electron gas (2DEG) realized in a GaAs/Al0.3Ga0.7As quantum well we observe changes in the Shubnikov-de Haas oscillations (SdHO) and in the Hall resistance for different sample geometries. We observe for each sample geometry a strong negative magnetoresistance around zero magnetic field which consists of a peak around zero magnetic field and of a huge magnetoresistance at larger fields. The peak around zero magnetic field is left unchanged for different geometries.

INFLUÊNCIA DA SAÚDE BUCAL SOBRE A QUALIDADE DE VIDA DE ADOLESCENTES ESCOLARES

Objetivo: Investigar o impacto da saúde bucal em relação à qualidade de vida de adolescentes escolares, associando-o às condições sociodemográficas. Métodos: Estudo de campo transversal e quantitativo desenvolvido em 2012, no município de Sumé-PB, com 184 adolescentes na faixa etária de 15 a 19 anos. Para avaliar o impacto, aplicou-se o questionário Oral Health Impact Profile (OHIP-14) aos participantes, enquanto para a obtenção dos dados referentes às condições sociodemográficas, os pais ou responsáveis responderam a um segundo questionário. Utilizou-se o teste Qui-quadrado para associar o impacto da saúde bucal sobre a qualidade de vida e as variáveis sociodemográficas pesquisadas, sendo considerados significativos com p<0,05. Resultados: Em geral, o impacto foi considerado fraco em 167 pesquisados (90,8%). “Dor física” foi a dimensão na qualidade de vida mais afetada pelas questões bucais entre aquelas que resultaram em impacto médio (22,8%; n=42). Apenas as variáveis “Situação do imóvel” e “Acomodação” associaram-se ao impacto geral (p<0,05). Os reduzidos percentuais de impacto geral forte (1,1%; n=2) relacionaram-se aos adolescentes cujas mães só estudaram até o ensino fundamental, ou às famílias que vivem com um salário mínimo ou menos (1,1%). Conclusão: Observou-se que as condições de saúde bucal apresentaram um impacto negativo fraco na qualidade de vida dos adolescentes investigados. As análises das condições sociodemográficas dos indivíduos relacionadas ao impacto geral da qualidade de vida relacionada à saúde oral associaram-se as variáveis “Situação do imóvel” e “Acomodação”.

Examining the growth and evolution of magnetic fields in clusters of galaxies: a numerical perspective

Magnetic fields are ubiquitous in galaxy cluster atmospheres and have a variety of astrophysical and cosmological consequences. Magnetic fields can contribute to the pressure support of clusters, affect thermal conduction, and modify the evolution of bubbles driven by active galactic nuclei. However, we currently do not fully understand the origin and evolution of these fields throughout cosmic time. Furthermore, we do not have a general understanding of the relationship between magnetic field strength and topology and other cluster properties, such as mass and X-ray luminosity. We can now begin to answer some of these questions using large-scale cosmological magnetohydrodynamic (MHD) simulations of the formation of galaxy clusters including the seeding and growth of magnetic fields. Using large-scale cosmological simulations with the FLASH code combined with a simplified model of the acceleration of cosmic rays responsible for the generation of radio halos, we find that the galaxy cluster frequency distribution and expected number counts of radio halos from upcoming low-frequency sur- veys are strongly dependent on the strength of magnetic fields. Thus, a more complete understanding of the origin and evolution of magnetic fields is necessary to understand and constrain models of diffuse synchrotron emission from clusters. One favored model for generating magnetic fields is through the amplification of weak seed fields in active galactic nuclei (AGN) accretion disks and their subsequent injection into cluster atmospheres via AGN-driven jets and bubbles. However, current large-scale cosmological simulations cannot directly include the physical processes associated with the accretion and feedback processes of AGN or the seeding and merging of the associated SMBHs. Thus, we must include these effects as subgrid models. In order to carefully study the growth of magnetic fields in clusters via AGN-driven outflows, we present a systematic study of SMBH and AGN subgrid models. Using dark-matter only cosmological simulations, we find that many important quantities, such as the relationship between SMBH mass and galactic bulge velocity dispersion and the merger rate of black holes, are highly sensitive to the subgrid model assumptions of SMBHs. In addition, using MHD calculations of an isolated cluster, we find that magnetic field strengths, extent, topology, and relationship to other gas quantities such as temperature and density are also highly dependent on the chosen model of accretion and feedback. We use these systematic studies of SMBHs and AGN inform and constrain our choice of subgrid models, and we use those results to outline a fully cosmological MHD simulation to study the injection and growth of magnetic fields in clusters of galaxies. This simulation will be the first to study the birth and evolution of magnetic fields using a fully closed accretion-feedback cycle, with as few assumptions as possible and a clearer understanding of the effects of the various parameter choices.

Magnetic and Acoustic Investigations of Turbulent Spherical Couette Flow

Title of dissertation: MAGNETIC AND ACOUSTIC INVESTIGATIONS OF TURBULENT SPHERICAL COUETTE FLOW Matthew M. Adams, Doctor of Philosophy, 2016 Dissertation directed by: Professor Daniel Lathrop Department of Physics This dissertation describes experiments in spherical Couette devices, using both gas and liquid sodium. The experimental geometry is motivated by the Earth's outer core, the seat of the geodynamo, and consists of an outer spherical shell and an inner sphere, both of which can be rotated independently to drive a shear flow in the fluid lying between them. In the case of experiments with liquid sodium, we apply DC axial magnetic fields, with a dominant dipole or quadrupole component, to the system. We measure the magnetic field induced by the flow of liquid sodium using an external array of Hall effect magnetic field probes, as well as two probes inserted into the fluid volume. This gives information about possible velocity patterns present, and we extend previous work categorizing flow states, noting further information that can be extracted from the induced field measurements. The limitations due to a lack of direct velocity measurements prompted us to work on developing the technique of using acoustic modes to measure zonal flows. Using gas as the working fluid in our 60~cm diameter spherical Couette experiment, we identified acoustic modes of the container, and obtained excellent agreement with theoretical predictions. For the case of uniform rotation of the system, we compared the acoustic mode frequency splittings with theoretical predictions for solid body flow, and obtained excellent agreement. This gave us confidence in extending this work to the case of differential rotation, with a turbulent flow state. Using the measured splittings for this case, our colleagues performed an inversion to infer the pattern of zonal velocities within the flow, the first such inversion in a rotating laboratory experiment. This technique holds promise for use in liquid sodium experiments, for which zonal flow measurements have historically been challenging.

Radio Frequency Superconducting Quantum Interference Device Meta-atoms and Metamaterials: Experiment, Theory, and Analysis

Metamamterials are 1D, 2D or 3D arrays of articial atoms. The articial atoms, called "meta-atoms", can be any component with tailorable electromagnetic properties, such as resonators, LC circuits, nano particles, and so on. By designing the properties of individual meta-atoms and the interaction created by putting them in a lattice, one can create a metamaterial with intriguing properties not found in nature. My Ph. D. work examines the meta-atoms based on radio frequency superconducting quantum interference devices (rf-SQUIDs); their tunability with dc magnetic field, rf magnetic field, and temperature are studied. The rf-SQUIDs are superconducting split ring resonators in which the usual capacitance is supplemented with a Josephson junction, which introduces strong nonlinearity in the rf properties. At relatively low rf magnetic field, a magnetic field tunability of the resonant frequency of up to 80 THz/Gauss by dc magnetic field is observed, and a total frequency tunability of 100% is achieved. The macroscopic quantum superconducting metamaterial also shows manipulative self-induced broadband transparency due to a qualitatively novel nonlinear mechanism that is different from conventional electromagnetically induced transparency (EIT) or its classical analogs. A near complete disappearance of resonant absorption under a range of applied rf flux is observed experimentally and explained theoretically. The transparency comes from the intrinsic bi-stability and can be tuned on/ off easily by altering rf and dc magnetic fields, temperature and history. Hysteretic in situ 100% tunability of transparency paves the way for auto-cloaking metamaterials, intensity dependent filters, and fast-tunable power limiters. An rf-SQUID metamaterial is shown to have qualitatively the same behavior as a single rf-SQUID with regards to dc flux, rf flux and temperature tuning. The two-tone response of self-resonant rf-SQUID meta-atoms and metamaterials is then studied here via intermodulation (IM) measurement over a broad range of tone frequencies and tone powers. A sharp onset followed by a surprising strongly suppressed IM region near the resonance is observed. This behavior can be understood employing methods in nonlinear dynamics; the sharp onset, and the gap of IM, are due to sudden state jumps during a beat of the two-tone sum input signal. The theory predicts that the IM can be manipulated with tone power, center frequency, frequency difference between the two tones, and temperature. This quantitative understanding potentially allows for the design of rf-SQUID metamaterials with either very low or very high IM response.