A. Structure and morphology

For combustion synthesis of BaLi₂Ti₆O₁₄ titanate, several trials were made to optimize the amount of lithium precursor, fuel, pre-calcination temperature, and final calcination temperature. To compensate for possible lithium volatisation at high temperature, a slightly excess amount (0–5%) of lithium precursor can be used. Varying the Li precursor content, the final BaLi₂Ti₆O₁₄ product was formed along with impurities like Li₄Ti₇O₁₆, BaTiO₃, and TiO₂ (Supplementary Material, Fig. S1). Using a pre-calcination temperature of 500 °C and final annealing at 800 °C, the phase pure titanate product was obtained. Despite several attempts, we were unable to get pure titanate phase by combustion synthesis employing alternate fuels such as ascorbic acid, citric acid, stearic acid, sucrose, and urea. The exothermicity of the combustion reaction can be altered by varying the oxidant to fuel ratio. In the current work, stoichiometric amounts of Ba/Li/Ti precursors (as oxidants) and glycine (as fuel) were used involving the final annealing at 800 °C for 2 h.

Rietveld refinement confirmed the phase purity of combustion synthesized BaLi₂Ti₆O₁₄ (χ² = 3.82%, R factor = 6.69%) (Fig. 1). When compared to solid-state synthesis involving 950 °C calcination for 12 h, solution combustion synthesis yielded the desired titanate with the lowest reported calcination temperature (800 °C) and calcination duration (2 h). This titanate stabilizes into an orthorhombic crystal structure (symmetry: Cmca) with lattice parameters a = 16.565(2) Å, b = 11.265(5) Å, and c = 11.584(0) Å similar to earlier reports.^{Reference Dambournet, Belharouak and Amine6,Reference Lin, Li, Shao, Shui, Wang, Long, Ren and Shu13} The atomic coordinates and bond valence sums of all atoms gave a GII value of 0.147(2) in agreement with the value of 0.102(8) found for single crystal literature data. It validates the stability of the chemical environment of the refined structural model for further study. Table I lists the crystallographic parameters, atomic coordinates, and bond valence sum values obtained from Rietveld refinement.

FIG. 1.

Rietveld refinement of XRD (Cu K_α) pattern of BaLi₂Ti₆O₁₄ formed by solution combustion synthesis. The experimental data points (red hollow dots), calculated pattern (black line), their difference (blue line), and Bragg reflections of the orthorhombic Cmca phase (black bars) are shown. The inset shows the crystal structure along c axis of the refined cell having partly filled Li tunnels surrounded by rutile tunnels transversely reinforced by TiO₆ pairs (green). (Ba = blue balls, TiO₆ octahedra = distorted green squares, and Li = black balls).

TABLE I.

Crystallographic parameters determined by the Rietveld refinement of high-resolution XRD data (λ = 1.5418 Å) of BaLi₂Ti₆O₁₄ at 25 °C. (U _iso = equivalent isotropic displacement parameters) (BVS = bond valence sum).

^a The high calculated theoretical density (4.09 g/cm³) in BaLi₂Ti₆O₁₄ can favor high tap density energy density.

^b GII is the root mean square bond valence mismatch of the bond valence sums.

The three-dimensional framework of BaLi₂Ti₆O₁₄ is built from layers TiO₆ octahedra having tunnels partly filled with Li ions. One set of TiO₆ octahedra containing Ti(1) and Ti(2) share opposite edges of their equatorial plane to form a rutile-like chain along c direction. Another similar set of TiO₆ octahedra reinforce these rutile chains by alternately sharing adjacent edges along b and c axis. It forms corner-sharing blocks of 3 edge-sharing TiO₆ octahedra (along a direction) having zig-zag tunnels partially filled with Li ions. Barium atoms occupy 11 coordinated 8f sites between the layers. Despite high unit cell volume, mean Ba–O bond distance (2.87 Å) in our case was similar to the value (2.889 Å) obtained from the single crystal study,^{Reference Koseva, Chaminade, Gravereau, Pechev, Peshev and Etourneau11} both of which were significantly lower than the sum of the ionic radii (2.97 Å). This stems from the residual compressive strain, which was further validated by a high bond valence value of Ba (2.11) (Table I).

The short annealing time led to less aggressive grain growth so as to form small primary particles in the size range of 150–300 nm (Fig. 2). They agglomerated to form secondary micrometric (<5 μm) particles. Characteristic to combustion reaction involving gaseous by-products, porous morphology was developed that can be beneficial for electrolyte wetting and efficient Li⁺ diffusion (with shorter diffusion length) during battery cycling. The BET specific surface area was found to be ∼8 m²/g (Supplementary Material, Fig. S2). HRTEM images [Fig. 2(f)] show $\left( {\bar{1}\bar{1}\bar{1}} \right)$ plane with a dislocation-free interplanar spacing of 7.1725 Å. The corresponding selected area diffraction pattern [SAED, Fig. 2(h)] confirmed phase purity with distinct points indexed to $\left[ {\bar{1}\bar{1}\bar{1}} \right]$ zone axis. Figure 2(g) shows (311) and (112) planes with an interplanar spacing of 4.54 Å and 4.94 Å, respectively. Due to their small difference, several edge dislocations can be seen. Elemental mapping [Fig. 2(j)] showed a uniform distribution of constituent Ba, Ti, and O elements in BaLi₂Ti₆O₁₄ titanate.

FIG. 2.

(a–c) SEM images of as-prepared BaLi₂Ti₆O₁₄ showing porous morphology. (d and e) TEM images indicating nanometric particle size. (f and g) High resolution TEM images showing fringes indexed to atomic plane and dislocations showing extra half plane at the interface between two planes. (h) SAED pattern indexed to the orthorhombic Cmca crystal structure. (i and j) Bright field STEM images and elemental mapping showing homogenous distribution of barium, titanium, and oxygen.

B. BVSE analysis

During lithiation, BaLi₂Ti₆O₁₄ reduces one (or more) Ti from Ti⁴⁺ to Ti³⁺ while intercalating xLi⁺ to form BaLi_2+xTi₆O₁₄. It involves diffusion of Li⁺ ions that is controlled by the crystal structure and surrounding potential. As Li⁺ ions intercalate at x = 0 in BaLi_2+xTi₆O₁₄ proceeding through partially delithiated structures (0 < x < 6), the chemical environment around Li⁺ gets gradually modified and determines its iterative distribution and movement. The key is to build an accurate lowest energy chemical environment from experimentally reliable and a hence relevant local structure consisting of refined atomic positions. Environment modeling can be done using several ways such as computationally expensive ab initio methods. Similar results can be obtained by using simpler BVSE calculations.^{Reference Adams and Rao22–Reference Chen and Adams24} BVSE method, working on force field of the chemical environment, intuitively identifies interstitial or vacancy sites of the lowest BVSE of lithium (E _BVSE Li⁺), stitching them by intersite hopping barrier energies into one or multidimensional threads of lithium ion migration pathways. Based on blocking due to defects/impurity or change of barrier energies with composition x, different pathways can activate. For consistency, migration pathways found by computationally expensive ab initio methods have been found to match those by the bond valence approach in suites of electrode materials.

In the pristine delithiated state (x = 0), tunnels directed along c axis host 16 Li⁺ ions (2 per formula unit as Z = 8) completely occupying single 16 g empty tetrahedral sites, along with two interstitial sites i2 and i3. These tunnels are created by TiO₆ rutile chains interconnected alternately above and below by sets of edge-sharing TiO₆ blocks (Fig. 1, inset). Inside the tunnels, along with i2 and i3 sites, Li1 sites form a one-dimensional (1D) zigzag migration pathway Li1–i2–Li1–i3 that interconnect to neighboring 1D channels by Li1–i4–Li1 2D paths forming a two-dimensional (2D) network spread in the ac or (010) plane [Fig. 3(a) showing isosurface as pathways along with position of the environmentally most perturbing ion]. The 2D nets are linked in the bc plane (100) by Li1–i3–i3–Li1 three dimensional (3D) paths [Fig. 3(b)]. For a clearer understanding, complete Li⁺ ion migration pathways overlapping with the crystal structure is illustrated in Fig. 3(c).

FIG. 3.

BVSE model showing lithium ion migration pathways in BaLi₂Ti₆O₁₄. The pathways shown are two isosurfaces (red and yellow) (a) (010) plane having one dimensional and two dimensional paths involving i2, i3, and i4 interstitials in half of the unit cell (to avoid overlap), (b) (100) plane having three dimensional pathways through i4 and (c) (001) plane showing 1D path running along the lithium filled tunnel, rutile chain (one labelled TiO₆), and (d) energy landscape showing migration barriers across different interstitials connected 1D, 2D, and 3D paths. The barrier for i3 access by Li1 is smaller than i2, hence as the energy corresponding to the isosurface in red is enough for i3 access but just exceeds the i2 barrier, i2 is not connected to lithium while i3 is. All illustrations have been presented using VESTA.^{Reference Momma and Izumi25}

From the energy plot in Fig. 3(d), it can be seen that i2 and i3 interstitials have dissimilar site energies unlike the case of SrLi₂Ti₆O₁₄.^{Reference Dayamani, Shinde, Chaupatnaik, Rao, Adams and Barpanda9} Despite the compressive strain of the M (M = Sr/Ba) atom in BaLi₂Ti₆O₁₄ (bond valence of Ba = 2.11) that is absent in SrLi₂Ti₆O₁₄ (bond valence of Sr = 1.75), the larger ionic size of Ba significantly increases the b and c lattice parameters in comparison to a. Ba is closer towards the i2 site while being far and opposite from i3. This smaller Ba-i2 distance [Fig. 3(a)] increases the i2 energy significantly from i3. Hence the 2D links are difficult to operate in the case of Ba with only one interstitial as compared to two interstitials in Sr. Concomitant to the case of Sr, with BVSE ignoring the mobile Li–Li repulsions, though interstitials i2 and i3 exist at energies much lower than Li1, stay unfilled because of Li1 stabilization in i1 by reduced repulsion among mobile Li⁺ ions. (The nomenclature is ix if the site is unfilled and Li_x if the site is filled by Li).

The migration pathways calculated for BaLi₂Ti₆O₁₄ have low barriers of 0.414 eV for 1D links with higher barriers of 1.281 and 1.574 eV for the 2D and 3D links [Fig. 3(d)]. Hence, most likely it is a one-dimensional ionic conductor, rarely activating the 2D and 3D pathways in cases when 1D pathways are blocked. It is in contrast to isostructural SrLi₂Ti₆O₁₄ which is at discretion to freely use 1D, 2D, and 3D pathways with barrier energies lower than 0.69 eV making it a 3D ion conductor. Following the energy landscape, a maximum of four Li⁺ will insert sequentially in i3, i2, i4, and i5 sites, each of which can hold one lithium. This is if progressive lithium insertion does not drastically energize the chemical environment. Otherwise, based on the lowest energy at that composition, Li⁺ ions will dynamically rearrange altering the energy landscape. Overall, BaLi₂Ti₆O₁₄ was found to be a one-dimensional Li⁺ conductor. The ionic conductivity was further analyzed by experimental ac conductivity measurement as described below.

C. AC (ionic) conductivity

Unlike graphitic anode where Li⁺ ions intercalate along with electrons without reducing anything, ionic conductivity is important in anodes involving transition metal ions. Ionic conductivity is more important than electronic conductivity as diffusion governs/limits the rate of (de)lithiation during battery cycling. Ionic conductivity was measured using BaLi₂Ti₆O₁₄ pellets having 94.24% (3.83 g/cm³) of theoretical density (4.069 g/cm³). Figure 4(a) shows the high temperature complex Nyquist plots from 200 to 300 °C at intervals of 10° C taken to understand the process involved during ionic conduction and correlate it with that obtained from the earlier BVSE analysis. For temperatures up to 250 °C, depressed semicircles were observed. A tail at the low-frequency end suddenly appeared in Nyquist plots above 250 °C indicating fast ion mobility which gets blocked as ionic charges accumulate at the sample–electrode interface. All depressed semicircles diminished at higher temperatures indicating an increase in conductivity.

FIG. 4.

Electrical measurements of AC conductivity. (a) Complex impedance plots for BaLi₂Ti₆O₁₄ recorded in the frequency range of 40 Hz–110 MHz at high temperatures and the corresponding fit, using ZView, as per the equivalent circuit in the inset, (b) AC conductivity versus frequency and fit based on Jonscher’s power law. The individual resolved and combined bulk resistance values from fits in (a) and Jonscher’s dc conductivity values from fits in (b) have been used to calculate the activation energy in graphs (c) and (d).

The observed values were fitted to an equivalent circuit consisting of a pair of parallel Debye-like resistor-constant phase element (representing impedance contribution from grain and grain boundary) and a constant phase element (representing the impedance contribution from the electrode–grain interface) connected in series. Though semicircles of impedance spectra seem as a single depressed semicircle, two overlapping semicircles from grain contribution at higher frequency and grain boundary contribution at lower frequency regions exist, which can be resolved slightly more clearly at higher temperatures. This difficulty in resolution indicates overlapping time constants for grain and grain boundary relaxation. The bulk resistance R was calculated from a parallel combination of resistance values of grain (R _g) and grain boundary (R _gb) obtained after fitting. Bulk conductivity (σ) in S/cm was calculated using σ = 1/ρ = l/(R × A), where l is the thickness of the pellet and A is the disc area of the pellet in cm. The grain, grain boundary, and bulk conductivity at 300 °C were 2.16 × 10⁻⁴ S/cm, 2.44 × 10⁻⁵ S/cm, and 2.41 × 10⁻⁴ S/cm, respectively. The activation energy for ion movement was calculated from Arrhenius relation of the temperature dependence of dc conductivity:

(2)$$\sigma = {A_0}\exp \left( { - {{{E_{\rm{a}}}} \over {kT}}} \right)\quad ,$$

where A ₀ is the pre-exponential factor, E _a is the activation energy in eV, k is the Boltzman’s constant in eV/K, and T is the temperature in K.

The individual grain and grain boundary conductivity values (in S/cm) at different temperatures obtained after fitting the impedance spectra have been plotted versus 1000/T [Fig. 4(c)]. The grain conductivity shows a large jump above 250 °C with a corresponding change in slope while the grain boundary conductivity follows a single trend. Two temperature regions of grain and entire temperature region of grain boundary conductivity were fitted. Grain boundary activation energy of 0.82 eV was obtained. Grain activation energy corresponding to 0.7 and 1.268 eV for high and low temperature regions was found. It indicates the occurrence of an internal switch in ionic conduction mechanism within the grain from 2D to 1D, increasing ion mobility, possibly from clearing away of defects in the 1D pathways at a higher temperature. It is justified from an excellent agreement of the activation energy obtained from BVSE where calculated values for the 1D path (0.414 eV) match the high-temperature values (0.7 eV), while values for the 2D path (1.281 eV) match the low-temperature values (1.268 eV) of grain conductivity.

Figure 4(b) shows the variation of ac conductivity with frequency at 10 °C intervals from 200 to 300 °C. AC conductivity increases with an increase in temperature indicating higher mobility, while at temperatures above 240 °C, it varies with frequency as per Jonscher’s universal dielectric response depicting localized ionic motion.^{Reference Jonscher26}

(3)$$\sigma \left( \omega \right) = {\sigma _{{\rm{dc}}}} + A{\omega ^n}\quad ,$$

where σ_dc is the low frequency or dc conductivity, A is the preexponential temperature-dependent factor, and n is the degree of interaction between mobile ions. Solid black lines in Fig. 4(b) show the fitting of experimental data to Jonscher’s law at each temperature from which the magnitudes of n and σ_dc are obtained. Plotted versus 1000/T [Fig. 4(d)], the Jonscher’s dc conductivity σ_dc shows a sharp change in slope at 250 °C. The activation energies of Jonscher’s dc and bulk conductivity are 0.57 and 0.68 eV at high temperature and 3.22 and 1.25 eV at low temperature. Except for the low temperature value of Jonscher’s dc conductivity, the change of activation energy signals a change in ion conduction mechanism from 2D to 1D as discussed earlier.

D. Electrochemical performance

The diffusion behavior in BaLi₂Ti₆O₁₄ was further analyzed by cyclic voltammetry (CV). Figure 5(a) shows potentiostatic cyclic voltammograms at scan rates ranging from 0.1 to 1 mV/s. Pristine BaLi₂Ti₆O₁₄ reduces while intercalating Li⁺ at 1.37–1.29 V (R1) and 1.16–1.11 V (R2) and oxidizes while deintercalating Li⁺ at 1.17–1.28 V (O1) and 1.48–1.55 V (O2). O1 can be resolved to two peaks which merge at higher scan rates. Figure 5(b) shows the variation of current for four redox peaks (R1, R2, O1, and O2) with square root of the scan rate in (V/s)^0.5. Linearly fit with zero intercepts R1, R2, O1, and O2 gave slopes of −0.076, −0.028, 0.030, and 0.109, respectively. These slopes were used to calculate the diffusion coefficient as per Randal Sevchik equation.^{Reference Ding, Xu, Yao, Wegner, Fang, Chen and Lieberwirth27}

(4)$${i_{\rm{p}}} = 2.69 \times {10^5} \times {n^{1.5}}A{D^{0.5}}C{\nu^{0.5}}\quad ,$$

where the i _p, n, A, D, C, and ν represent the peak current in A, number of electrons utilized in the peak, electrode area of 2.324 cm², diffusion coefficient in cm²/s, the concentration of the diffusing ion for each peak in mol/cm³, and scan rate in V/s. C is the concentration of lithium ion in a unit cell of BaLi₂Ti₆O₁₄. For 1 Li⁺ ion in a unit cell, C is the ratio of moles of Li to the volume of the unit cell (2164.17 ų) calculated from refinement, i.e., 0.006137 mol/cm³. C can then be written as

(5)$$C = 0.006137\left( {\# \;{\rm{of}}\;{\rm{L}}{{\rm{i}}^ + }} \right)\;{\rm{mol/c}}{{\rm{m}}^3}\quad .$$

FIG. 5.

(a) CV of BaLi₂Ti₆O₁₄ versus lithium half-cell cycled at current rates from 0.1 to 1 mV/s. Four distinct redox peaks were observed named as R1, R2, O1, and O2 in the sequence of their appearance. (b) Linear fit of the four peak currents to square root of the scan rate. (c) Schematic of calculated diffusion coefficient from slope obtained from fits in (b).

All Li⁺ ions in the unit cell are considered. This includes pre-existing two lithium ions in BaLi₂Ti₆O₁₄ in addition to n lithiums exchanged corresponding to n electrons during the redox process. The reason why the two pre-existing structural lithiums are included is because the BVSE analysis confirmed their mobility as they were part of 1D migration pathways having a low energy barrier of 0.414 eV. The number of electrons, n, at each peak depends on the number of Li⁺ ions exchanged during that peak. n for R1 and O2 is assumed to be 1.28 based on the intake of around 52 mA h/g observed in 1st charge at 1.41 V in 4a/4b/8c sites as discussed later. Similarly, n for R2 and O1 is assumed to be 1.24 based on 50 mA h/g charge capacity stored in the slopes other than during and after the plateau in the 1st charge. Thus, the cumulative numbers of lithium ions present in the unit cell in R1, R2, O1, and O2 are 3.38 (2 + 1.28), 4.62 (2 + 1.28 + 1.24), 3.38 (4.62 − 1.24), and 2 (3.38 − 1.28). Although the O1 and R2 peaks are smaller than the sharp peaks, the area under the curve remains the same confirming equal charge exchange (close to 1.3 lithiums). Based on the above discussions, the diffusion coefficients of lithium ions at R1, R2, O1, and O2 are compared in Fig. 5(c). It is highest for the O2 peak as it has the least number of lithiums (2) amongst all redox peaks.

Diffusion coefficient values in the range of 10⁻¹¹ to 10⁻¹² cm²/s [Fig. 5(c)] were obtained for combustion prepared BaLi₂Ti₆O₁₄. The ionic diffusivity can be enhanced by four standard ways: (i) changing crystallinity, (ii) doping, (iii) conductive coating, and (iv) particle downsizing. Combustion synthesis delivers carbon-coated small particles, which is favorable to obtain high diffusion coefficient. Following the diffusion analyses, the electrochemical performance was measured in standard Li- and Na-half-cell architectures.

The galvanostatic voltage–capacity profiles of combustion-made BaLi₂Ti₆O₁₄ are shown in Fig. 6 cycled (i) at C/10 rate between a narrow potential window (1–2 V) and (ii) at C/20 rate between a wider potential window (0.5–2 V). In the former case, out of 120 mA h/g stored in the first cycle, 16 mA h/g was lost to the parasitic reaction, 62.6 mA h/g (1.55 Li⁺) was stored in the 4a/4b/8c sites at 1.41 V (inset) while 41 mA h/g (1 Li⁺) was stored in the vacant 8f sites. In the subsequent cycles, out of 87 mA h/g (2.12 Li⁺) extracted during the first charge, 53 mA h/g (1.3 Li⁺) could be extracted at 1.41 V [the inset of Fig. 6(a)] from the partly lithium filled 4a/4b/8c sites. Excellent capacity retention was recorded for over 120 cycles.

FIG. 6.

Cycle life performance of BaLi₂Ti₆O₁₄ versus lithium (a) at C/10 rate for a narrow voltage window of 1–2 V and (b) at C/20 for a wider voltage window of 0.5–2 V. A C/20 slow formation cycle precedes both tests. The insets show capacity retention, dQ/dV versus V, and first cycle losses.

When cycled in the wide 1.5 V potential window (0.5–2 V), the first discharge of 176 mA h/g was observed [the inset of Fig. 6(b)], out of which the first 17 mA h/g (0.4 Li⁺) was lost to the parasitic reaction, next 63 mA h/g (1.57 Li⁺) was stored in vacant 4a/4b/8c sites of lithium tunnels at 1.41 V while remaining 96 mA h/g (2.38 Li⁺) was stored collectively at (i) 2 different voltage regions in 8f sites and (ii) consumed to possible SEI formation. 111 mA h/g of capacity (2.75 Li⁺) could be extracted in the first charge, out of which only 51.78 mA h/g (1.28 Li⁺) could be removed at 1.43 V (low polarization of 0.02 V) from the partly filled 4a/4b/8c sites. 78% of total first charge capacity (87 mA h/g or 2.12 Li⁺) could be retained or reversibly extracted in the 50th charge cycle [Fig. 6(b) inset]. These results indicated (a) the same 1.3 Li⁺ can be extracted in the first plateau independent of potential window indicating reversibility in 4a/4b/8c sites remains unaffected and (b) reducing lower cut-off voltage significantly affects Li⁺ reversibility/capacity retention whose origins could be gradual lithium loss to unstable SEI formation and/or trapping in the inherent weakly reversible 8f sites.

Following, the rate performance of combustion prepared BaLi₂Ti₆O₁₄ was assessed by varying the current rate from C/10 to 3C. The corresponding capacity retention and voltage–capacity profiles are summarized in Fig. 7. Charge capacities of 96 mA h/g (2.38 Li⁺), 80 mA h/g (1.98 Li⁺), 71 mA h/g (1.76 Li⁺), 64 mA h/g (1.58 Li⁺), 58 mA h/g (1.44 Li⁺), and 55 mA h/g (1.36 Li⁺) were recorded at C rates of C/10, C/5, C/2, 1C, 2C, and 3C, respectively. Charge capacity reverted back to 77 mA h/g (1.9 Li⁺) upon decreasing the C rate to C/5.

FIG. 7.

Rate performance of BaLi₂Ti₆O₁₄ versus lithium at higher current rates up to 3C (726 mA/g). The inset shows voltage profile and their polarization at different C rates.

Overall, the electrochemical performance of BaLi₂Ti₆O₁₄ prepared using rapid solution combustion synthesis was better than solid-state prepared BaLi₂Ti₆O₁₄ without ball milling. It could be due to the formation of smaller particles with porous morphology. In our case, the electrochemical performance can be further improved by ball milling, doping (that improves electronic conductivity), and process optimization (co-fuel/oxidant-fuel ratio) until it can fully utilize the 1D low energy pathways.

Following the efficient (de)lithiation, the presence of open channels inspired us to gauge possible (de)insertion of Na⁺ ions into BaLi₂Ti₆O₁₄. When cycled in a Na half-cell architecture, we observed reversible Na⁺ insertion for the first time (Fig. 8). The first discharge capacity was 195 mA h/g (4.83 Na⁺ insertion) with a plateau ∼0.15 V. In the subsequent cycles, a reversible capacity of 36 mA h/g (0.89 Na⁺) was observed with an average redox potential of 0.6 V. This preliminary result hints at the use of MLi₂Ti₆O₁₄ family of compounds for sodium-ion insertion with proper material optimization.

FIG. 8.

Sodium activity in BaLi₂Ti₆O₁₄ at 20 mA h/g in 1 M NaPF₆ in 1:1 EC:DEC in the voltage range of 0.1 and 2.5 V. The inset shows capacity retention, first cycle loss, and dQ/dV versus V.