United States

Magnetic navigation systems (MNSs) capable of generating various types of magnetic fields have been widely studied for the effective manipulation of biomedical magnetic robots [1-4]. Previous researchers proposed a geometrically and electrically efficient MNS comprised of a simple triad of electromagnetic coils (TEC) that can generate two-dimensional (2D) magnetic robot motions [3]. However, the TEC suffered from the restricted working area of the magnetic robot due to the nonlinearly distributed magnetic field throughout the system. In this research, we propose a method to precisely generate 2D magnetic torque and force of a magnetic robot regarding the robot&#39;s current position inside the TEC. Since the TEC is composed of three identical circular coils located at the vertices of an equilateral triangle, the magnetic field at a certain position inside the TEC can be expressed in terms of the three coils&#39; currents and a position vector by using Jacobian and transformation matrices [4]. From this expression, we established several constraint equations to optimally manipulate the magnetic robot’s 2D motions at any position within the TEC. Fig. 1 shows several simulated magnetic fields of the TEC to propel a magnetic robot (red dots) in the direction of the x-axis. It shows that the TEC can generate a relatively uniform magnetic force within the nearby area of a magnetic robot, making it possible to move the magnetic robot in the desired direction (Fig. 2). We then constructed the TEC and a closed-loop visual tracking system to demonstrate that the TEC can manipulate a magnetic robot inside the TEC in a closed-loop manner. The results showed that the proposed method and system can improve the magnetic robot’s admissible working area in the TEC effectively and precisely.References: [1] M. Kummer, J. Abbott, and B. E. Kratochvil et al., IEEE Transactions on Robotics, Vol. 26, pp. 1006-1017, 2010. [2] Q. Cao, X. Han, B. Zhang, and L. Li et al., IEEE Transactions on Applied Superconductivity, Vol. 22, 2012. [3] H. Lee and S. Jeon, AIP Advances. Vol. 10, 2020. [4] G. Go, H. Choi, and S. Jeong et al., Sensors and Actuators, Vol. 205, P. 215-223, 2014.

MMM 2020

Development of a Position-Based Precise Magnetic Torque and Force Generation Method Using a Triad of Electromagnetic Coils for Biomedical Magnetic Robot Applications

The 65th Annual Conference on Magnetism and Magnetic Materials (MMM 2020) was held in a virtual form. The COVID-19 pandemic, with its associated uncertainty about international travel and the ability to gather in large groups, meant that the Steering Committee decided to cancel the face-to-face meeting that had been planned in Palm Beach, Florida November 2-6, 2020.

We were reluctant to give up the benefits of meeting face-to-face, and know that shifting to a virtual conference is a significant change for this community. Nevertheless, we are first and foremost concerned about the health and safety of all our Conference participants and their loved ones on returning home after travelling.

The MMM 2020 Steering Committee was working exceptionally hard, along with our sponsors, the IEEE Magnetics Society and AIP Publishing, to make this event serve the needs of our participants, who live in time zones all around the globe, in the best possible way. 

You need a ticket to explore this event.

If you have already registered for MMM 2020, and you have already created an account, then you
just need to login (above).

If you have already registered for MMM 2020, but you have not yet created an account, then you
need to click on &quot;create an account&quot; (above) and then login.

Forgot your password? Go to
https://www.underline.io/log-in

If you have not yet registered for MMM 2020 but would like to do so, go to
https://magnetism.org/register

Still need help? Contact our team!


Magnetic navigation systems (MNSs) capable of generating various types of magnetic fields have been widely studied for the effective manipulation of biomedical magnetic robots [1-4]. Previous researchers proposed a geometrically and electrically efficient MNS comprised of a simple triad of electromagnetic coils (TEC) that can generate two-dimensional (2D) magnetic robot motions [3]. However, the TEC suffered from the restricted working area of the magnetic robot due to the nonlinearly distributed magnetic field throughout the system. In this research, we propose a method to precisely generate 2D magnetic torque and force of a magnetic robot regarding the robot's current position inside the TEC. Since the TEC is composed of three identical circular coils located at the vertices of an equilateral triangle, the magnetic field at a certain position inside the TEC can be expressed in terms of the three coils' currents and a position vector by using Jacobian and transformation matrices [4]. From this expression, we established several constraint equations to optimally manipulate the magnetic robot’s 2D motions at any position within the TEC. Fig. 1 shows several simulated magnetic fields of the TEC to propel a magnetic robot (red dots) in the direction of the x-axis. It shows that the TEC can generate a relatively uniform magnetic force within the nearby area of a magnetic robot, making it possible to move the magnetic robot in the desired direction (Fig. 2). We then constructed the TEC and a closed-loop visual tracking system to demonstrate that the TEC can manipulate a magnetic robot inside the TEC in a closed-loop manner. The results showed that the proposed method and system can improve the magnetic robot’s admissible working area in the TEC effectively and precisely.References: [1] M. Kummer, J. Abbott, and B. E. Kratochvil et al., IEEE Transactions on Robotics, Vol. 26, pp. 1006-1017, 2010. [2] Q. Cao, X. Han, B. Zhang, and L. Li et al., IEEE Transactions on Applied Superconductivity, Vol. 22, 2012. [3] H. Lee and S. Jeon, AIP Advances. Vol. 10, 2020. [4] G. Go, H. Choi, and S. Jeong et al., Sensors and Actuators, Vol. 205, P. 215-223, 2014.

poster

Rubber bands can be found in many children’s toys. For example, one winds up a rubber band and when released it drives an airplane or a racing car. If instead we wind up a tiny magnetic structure, where the magnetization twists and stores energy, then that energy will be released upon unwinding. In this theoretical work we consider a very soft, low-exchange film (80-400 atomic layers of Cu-doped Permalloy [1]), which is pinned at the bottom and coupled to a high-exchange ferromagnetic film on top. We apply a magnetic field which is rotated in the plane of the film. This generates a slow twisting of the soft material with the high-exchange ferromagnet acting as a handle as is seen in Fig. 1. We use both a local self-consistent mean-field theory [2] as well as a dynamic approach with the Landau-Lifshitz equation, with results matching well for both. With these models we have been able to generate multiple twists (over 10) in these films. At some point an instability occurs that prevents additional twists. Often the instability causes one or more twists to be released from the structure. The tightness of the twists as well as how many twists are formed before snapping depends on the length of the structure, the strength of the applied field, the size of the of the ferromagnetic cap, the exchange coupling strength, and the saturation magnetization. For example, increasing the exchange coupling decreases the number of twists that are stable in a given structure. The number of allowed twists also increases linearly with the length of the sample. Analytic energy minimization calculations, using a constant turn angle approximation, also allow us to estimate how many twists will fit into a given structure, matching reasonably well with numerical outcomes. The more twists in a given structure, the more energy it can hold. Once the twist is released or snapped, the system will unwind with frequencies in the gigahertz range. We study the frequency distribution as a function of the number of twists. Based on energy storage and this frequency generation we discuss potential applications.References: [1] N. Sorensen, R. E. Camley and Z. Celinski, J.Magn Magn Mat, vol. 477, p. 344, (2019). [2] R. E. Camley, Phys Rev B, vol. 35, p. 3608, (1987).

Mini Magnetic Propellers: Stability Considerations in Twisted Magnetic Structures

We analyze in this work the low temperature magnetic behavior of a series of Co/Ag multilayers formed by samples having reduced Co content per period (either one monolayer or half a monolayer) separated by several Ag monolayers (either 8 or 16). The multilayers were grown by molecular beam epitaxy on a clean MgO (001) surface, alternating Ag and Co deposition up to complete either 32 or 64 periods (in order to maintain constant in the series the nominal number of Co atoms per unit surface of the sample) and protected by a 3nm Ag cap layer. The samples were studied by X-ray reflectivity and diffraction, TEM, SQUID magnetometry, and ac susceptometry. From the results of our measurements we identified the following facts: 1) the occurrence in all the samples of stacking sequences corresponding to the nominal ones 2) the spatial discontinuity of the Co layers that corresponded to a quasi-monodisperse in-plane distribution of Co particles embedded in a Ag(001) matrix; 3) the observation of hysteresis at temperatures (always below 20 K) that depended on both the number of Co and Ag monolayers per period; 4) an ac field frequency dependence of the real part of the susceptibility, χ’, maxima temperatures corresponding to a Vogel-Fulcher behavior with a frequency shift parameter of the order in all the samples of 4 x 10-2; 5) the collapse (Figure 1 ), for each sample, of the experimental data corresponding to the ac field frequency, f, and temperature, T, variations of the imaginary part of the ac susceptibility, χ’’, in a single curve, G(x) with x = f[(T/Tg) − 1]−Ζν (where Tg and Ζν are the critical temperature and exponent, respectively), and verifying T χ’’ (f, T) = [(T/Tg) – 1]β G(x) [1]. Taken together, these results allow us to conclude that our three multilayers experience a phase transition (of the paramagnetic to superspin glass type) at the corresponding critical temperature and that the exponent β in the universal dynamic scaling law takes values dependent on the number of Co monolayers per period.References: [1] P. Beauvillain, C. Dupas, J.P. Renard, and P. Veillet, Physical Review B 29 (7), 4086 (1984).

Critical Magnetic Behavior in [Ag8/Co0.5]x64, [Ag16/Co1]x32 and [Ag8/Co1]x32 Epitaxial Multilayers

The study of magnetic systems, spin glasses, and their states, is an actively developed field of statistical physics. Spin glass models are studied in various scientific contexts, including experimental physics, condensed matter physics, theoretical physics, mathematical statistical physics and, more recently, probability theory [1 - 3]. The problem of finding the ground states (GS), i.e. states with the lowest energy is a very difficult task due to the macroscopic degeneracy and frustrations that are exist in various realistic spin glass models. In fact, it’s well known that finding the GS, of a spin glass in a three-dimensional lattice is an NP-complete problem, which means that this challenge is at least as difficult as the hardest problems of practical interest. The system of interest is the Edwards-Anderson (EA) spin glass model [4] in two-dimensional (2D) lattice, with bimodal distribution of bonds. In present work, for simulation we used a combination of our Hybrid Multispin Method (HMM) [5] and Long Short-Term Memory (LSTM) - an artificial recurrent neural network (RNN) architecture [6] with Boltzmann distribution. HMM’s advantages are the ability to increase the number of simultaneously gathered statistical-thermodynamic parameters without noticeable speed and efficiency loss, and the ability to search for configurations of the ground and low-energy states. The usage of machine learning in statistical physics began relatively recently but is developing rapidly. Neural networks have become popular tool due to the high speed of their training and accuracy of predictions. LSTM with Boltzmann distribution is useful for certain types of prediction that require the network to retain information for longer periods of time that traditional RNNs encounter. In our research, MC algorithm generated a training data set which is configurations of SG system. The trained predicted system configurations, which are using for the next computing cycle for HMM. Such iterative approach allows us to efficiently simulate the spin-glass systems even in the low-temperature regime, avoiding the critical slowing down that plague of usual MC simulation. This work was supported by the state task of the MSHE of Russia No. 0657-2020-0005.References: [1] A. Callison et al., New J. Phys., Vol. 21, 123022 (2019) [2] B. McNaughton et al., Phys. Rev. E, Vol. 101, 053312 (2020) [3] L. A. Fernandez, E. Marinari, V. Martin-Mayor, I. Paga, and J. J. Ruiz-Lorenzo, Phys. Rev. Vol. 100, 184412 (2019) [4] S. Edwards, P. Anderson, Phys. F: Metal Phys., Vol. 5 (5), p. 965 (1975) [5] A.G. Makarov et al., JETP Letters., Vol. 110 (10), pp. 702-706 (2019) [6] S. Hochreiter, J. Schmidhuber, Neural computation., Vol. 9(8), pp. 1735-1780 (1997)

Investigation of the spin glass density of states using a combination of Monte Carlo method and neural network

We are studying the effect that illumination by coherent resonant x-rays may have on magnetic domain memory (MDM) in a [Co / Pd] / IrMn multilayers [1-3]. MDM is the ability of the magnetic domains to retain their exact same domain topology upon field cycling. Earlier studies have suggested that under higher doses of x-ray illumination, the material may lose its existing MDM. To investigate this possible effect we have used both x-ray resonant magnetic scattering (XRMS) along with magneto-transport measurements [4,5] to track the exchange bias while the sample is illuminated with x-rays, as illustrated in Figure 1. Fig.1a shows magneto-transport data collected on the [Co/Pd] / IrMn multilayer from which we hope to extract information about exchange bias. A loss of exchange bias would indicate that the x-rays illumination dose may alter the strength of the internal exchange couplings between the ferromagnetic and antiferromagnetic layers. Ultimately it would alter the amount of MDM, which is measured via XRMS as illustrated in Fig. 1b, and represented in a correlation map such as Fig.1c.References: [1] K. Chesnel, A. Safsten, M. Rytting, and E.E. Fullerton, Nature Communications 7, (2016). [2] K. Chesnel, B. Wilcken, M. Rytting, S.D. Kevan, and E.E. Fullerton, New Journal of Physics 15, 023016 (2013). [3] K. Chesnel, J. Nelson, B. Wilcken, and S.D. Kevan, Journal of Synchrotron Radiation 19, 293 (2012). [4] C. Hurd, Hall Effect in Metals and Alloys (Springer, 2012). [5] L.J. van der Pauw, A Method of Measuring Specific Resistivity and Hall Effect of Discs of Arbitrary Shape (1958).

The effect of x-ray illumination on magnetic domain memory in [Co/Pd] / IrMn multilayers

One of the most important goals in the current research in molecular magnetism is to reach a deep understanding of the relaxation dynamics of Molecular NanoMagnets (MNMs), being this aspect of paramount importance in the design of new systems acting as classical [1, 2] or quantum bits [3]. Despite being the main source of relaxation in MNMs, very limited experimental investigations on their phonons have been done so far, yielding no information about their dispersions. Indeed, the key to construct a reliable model of phonon-induced relaxation dynamics in MNMs is to have access to phonon dispersions and eigenvectors, which are necessary for a quantitative evaluation of spin-phonon coupling coefficients. Here we exploit the four-Dimensional Inelastic Neutron Scattering (4D-INS) technique [4] to directly measure for the first time phonon dispersions in a prototypical molecular qubit, [VO(acac)2]. Vanadyl-based molecules are in fact currently emerging as archetypes of a new generation of molecular qubits with long coherence times even at high temperatures [3]. 4D-INS is a very powerful technique to study phonons, because it enables a direct access to phonon eigenvalues and eigenvectors, by measuring the four-dimensional scattering function in large portions of the reciprocal space. We detected both acoustic and optical branches in single crystals of [VO(acac)2] along different directions in the reciprocal space and we compared the results of this challenging experiment with state-of-the-art density functional theory (DFT) calculations (Fig.1). Simulations reproduce important features of the 4D-INS data, such as the presence of low-lying optical modes and anti-crossings between acoustic and optical phonons, providing insights on their effect on spin-phonon couplings [5].References: [1] C.A.P. Goodwin, F. Ortu, D. Reta, et al., Nature Vol. 548, p. 439, (2017). [2] F.-S. Guo, B. M. Day, Y.-C. Chen, et al., Science Vol. 362, p. 1400, (2018). [3] M. Atzori, L. Tesi, E. Morra, et al., J. Am. Chem. Soc. Vol. 138, p. 2154, (2016). [4] E. Garlatti, A. Chiesa, T. Guidi, et al., Eur. J. Inorg. Chem., p. 1106, (2019). [5] E. Garlatti, L. Tesi, A. Lunghi, et al., Nat. Commun. Vol. 11, p. 1751 (2020).

Unveiling Phonons in a Molecular Qubit with Four-Dimensional Inelastic Neutron Scattering and Density Functional Theory

The achievement of quantum error correction is a fundamental step towards the implementation of a scalable platform for quantum computation, in which logical units are made robust against errors. Here we show that magnetic molecules consisting of weakly interacting and individually addressable spins ½ centers are a very promising platform to implement repetition codes for quantum error-correction, where several spins ½ are used to encode an error protected logical qubit. In particular, we demonstrate that these challenging requirements can be fulfilled by the [Er-Ce-Er] molecular trinuclear coordination compound. We show that this complex can be exploited to efficiently implement the three-qubit quantum error correction code protecting against pure dephasing, the most important error in magnetic molecules [1]. The system is characterized by combining magnetometry, low-temperature specific heat and electron paramagnetic resonance measurements on both the trimer and the properly synthesized [LuCeLu] and [ErLaEr] analogues, which contain only one of the two types of magnetic ions. The resulting large and distinct g values, together with weak but non-negligible dipolar interactions, enable us to selectively address different transitions by a proper sequence of microwave pulses, implementing the three-qubit repetition code. The parameters obtained from the experimental characterization are finally used to perform thorough numerical simulations of the system dynamics, including the effect of pure dephasing in a Lindblad formalism. These show a remarkable gain of the quantum error correction procedure, thus proving the efficiency of the device in suppressing pure dephasing.References: [1] E. Macaluso, M. Rubin, D. Aguilà, A. Chiesa, L. Barrios, J. I. Martínez, P. J. Alonso, O. Roubeau, F. Luis, G. Aromí and S. Carretta, submitted.

A Heterometallic [LnLn’Ln] Lanthanide Complex as a Qubit with Embedded Quantum Error Correction

Low-dimensional (d=1,2) metallic systems with electron-hole symmetry are dominated at low temperatures by spin-flip fluctuations and enhanced antiferromagnetic interaction. They show, however, no long-range order at non-zero temperatures. The electronic structure is nevertheless substantially altered by the antiferromagnetic fluctuations and significantly deviates from the Fermi liquid. These changes are difficult to capture by the standard theoretical many-body techniques. We show that only many-body approaches including renormalizations of the interaction strength can successfully describe the low-temperature behavior of low-dimensional metals. We use the standard many-body perturbation expansion within which we resort to the reduced parquet equations introducing a two-particle self-consistency suppressing spurious phase transitions. They were introduced by the authors to describe quantum criticality in strongly correlated electron systems, Refs. [1-3]. We choose the one-dimensional Hubbard model and show that the solution at zero temperature and at half-filling contains an antiferromagnetic quantum critical point even if there is no long-range antiferromagnetic order at non-zero temperatures. Antiferromagnetic fluctuations open an energy gap in the spectral function that is proportional to the bare interaction strength. This behavior can be derived only if the quantum critical point and the renormalization of the interaction strength are treated consistently in both the thermodynamic response and spectral functions. That is, a balance between quantum and thermal fluctuations is maintained and consistency between one-particle and two-particle functions is accomplished.References: [1] V. Janiš, P. Zalom, V. Pokorný, and A. Klíč, Phys. Rev. B 100, 195114 (2019). [2] V. Janiš and A. Klíč, JPS Conference Proceedings 30, 011124 (2020). [3] V. Janiš, A. Klíč, and J. Yan, preprint arXiv:2006.08492 (2020).

Antiferromagnetic Fluctuations in One-Dimensional Hubbard Model

For nearly 100 years the localized-exchange-interaction, as represented by the Heisenberg Hamiltonian H = Σ J(r) Si. Sj, has been viewed as the primary underlying cause of magnetism in magnetic materials. An inescapable consequence of the Heisenberg Hamiltonian is that the magnetic dispersion relation at small k must go as ω∝k2 (for sufficiently short-range J(r)). Empirical evidence over the past 60 years has persistently given indications, however, that the measured dispersion-relation exponents in various ferromagnetic systems were invariably smaller than 2. The empirical evidence thus appears to be in direct conflict with the applicability of the Heisenberg Hamiltonian to explain spin ordering in real ferromagnets. The very first attempt to determine ω(k) at small k, using Spin Wave Resonance (SWR) showed near linear dependence,1 but the belief in the validity of theory was so strong that a way was always found to force agreement (e.g., by renumbering1,2 the resonances so as to obtain the expected ω∝k2 dependence). Such ad hoc renumbering, as well as postulation of magnetization and anisotropy variations within the films, etc., have been used to this day to insist that ω must be ∝k2. Neutron scattering (NS) has not been precise enough until recently to measure ω∝kd at small k. But a recent measurement3 on CrSiTe3 showed that d is indeed <2, despite the authors' attempt to fit the data to ω∝k2. Indirect evidence from precision NMR measurements had previously shown that empirically4 only ω∝ k; k1.5; and k1.75 for 3D ferromagnets, and ∝ k and k1.5 for 2D ferromagnets, appear to characterize simple ferromagnets. These behaviors may be understood as resulting from the exchange interaction being a collective-exchange effect, instead of an interatomic-exchange effect. If d<2, then 2D magnets can order spontaneously (the Mermin-Wagner theorem no longer applies).. This is consonant with the existence of high Tc magnetic monolayers and quasi-2D magnets, no matter the strength of the anisotropy fields. A full replacement model Hamiltonian is yet to be discovered. However, a simple spin-wave scheme can already be formulated from just the small-k approximations and is in excellent agreement with experimental evidence.References: 1. M.H. Seavey and P.E.Tannenwald, Phys Rev Lett 1, 168 (1958); J Appl Phys, 30, S227 (1959) - SWR in permalloy 2. P.K. Schwob and G.E. Everett, J Phys Colloque C1, 1066 (1971). Also, P.K. Schwob et al, PR B 10, 165 (1974) - SWR in EuS - 3. T.J. Williams et al, PR B 92, 144404 (2015) - NS in CrSiTe3 - 4. N. Bykovetz, J Appl Phys 55, 2062 (1984). Please note the obvious typographical error in the abstract, where n=0 should be replaced by n=∞

Is it Time to Sideline the Renowned Heisenberg Hamiltonian?

Ferromagnetic materials have played a significant role in magnetic memory devices. However, due to inherent limitations associated with stray field generation and vulnerability to external stimulation, antiferromagnetic materials are increasingly investigated as alternatives. Specifically, metallic tetragonal antiferromagnets with easy-plane magnetism are appealing because of the degenerate domains with 90° difference of Néel vector for binary information storage and the possibility of electrical switching. To understand the switching process, the magnetocrystalline anisotropy is essential, acting as an energy barrier between two degenerate states. In this study, we used density functional theory (DFT) to unravel the origin of the magnetocrystalline anisotropy in Fe2As, one of the metallic tetragonal antiferromagnets with easy-plane magnetism. The tetragonal structure allows two possible anisotropies out-of-plane and in-plane. Here, we calculate the spin-orbit interaction (SOI) in a DFT framework and the magnetic dipole-dipole interaction (MDD) based on the DFT ground state structure using a classical model. The computed out-of-plane |K1| anisotropy coefficient is 830 kJ/m3, composed to one-third of MDD and two-thirds of SOI, both of which exhibit the same two-fold symmetry. Conversely, the MDD contribution for in-plane anisotropy is almost imperceptible, leading to an in-plane |K22| anisotropy coefficient of 280 J/m3 with four-fold symmetry, due to SOI. The lowest frequency of a spin-wave is calculated using our anisotropy result, and it is predicted as 0.64 THz. Based on the frequency and anisotropy, the fast dynamics and high speed of switching are anticipated. However, low thermal stability might be observed due to the small in-plane anisotropy barrier. ** The Illinois Materials Research Science and Engineering Center is supported by the National Science Foundation MRSEC program under NSF Award Number DMR-1720633.

Ab initio Study of Magnetocrystalline Anisotropy in Antiferromagnetic Fe2As

The physics of the chiral magnet MnSi touches onto several fundamental questions in condensed-matter physics, from the stabilisation of exotic states like chiral skyrmions [1] to the interplay between localised and itinerant magnetism. Under mechanical pressure, or chemical pressure in the form of Fe doping in Mn1-xFexSi, a non-Fermi liquid behaviour sets-in which has been attributed to a chiral spin liquid and quantum fluctuations [2]. The origin and nature of these phenomena are still an open question, which we address by combining small angle neutron scattering and neutron spin echo spectroscopy [3-5]. Our results show that the destabilisation of the long range helical periodicity, close to the quantum critical points, is driven by the modification of the electronic state but with different microscopic mechanisms for mechanical and chemical pressure. Mechanical pressure enhances the itinerant electron character of MnSi leading to a softening of the magnetic moment. On the other hand, chemical doping increases the weight of achiral RKKY interactions. The resulting frustration explains both the destabilisation of long range helical periodicity and the enhancement of the precursor region with increasing doping.References: 1. S. Mühlbauer et al., Science 323, 915 (2009). 2. C. Pfleiderer, et al., Nature 414, 427 (2001) 3. L. J. Bannenberg et al., Phys. Rev B 100, 054447 (2019). 4. L. J. Bannenberg, et al., Phys. Rev. B 98, 184431 (2018). 5. C. Pappas, et al., in preparation

Premium content

Development of a Position-Based Precise Magnetic Torque and Force Generation Method Using a Triad of Electromagnetic Coils for Biomedical Magnetic Robot Applications

Next from MMM 2020

Mini Magnetic Propellers: Stability Considerations in Twisted Magnetic Structures

Similar lecture

THz Magneto-Optical Investigation of Quadrupolar Spin-Lattice Effects in Magnetically Frustrated Tb2Ti2O7