From cd58ca342fa989c3109f85cb678bf9144a4d7291 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Tue, 19 May 2020 12:53:10 -0400 Subject: More referee response changes. --- fig-stiffnesses.gplot | 22 ++++---- fig-stiffnesses.pdf | Bin 331750 -> 332053 bytes main.tex | 6 +-- referee_comments.txt | 85 +++++++++++++++++++++---------- referee_response_cutoff_curves.pdf | Bin 81698 -> 80978 bytes referee_response_cutoff_eigenvectors.pdf | Bin 9171 -> 9073 bytes referee_response_cutoff_parameters.pdf | Bin 9444 -> 9194 bytes 7 files changed, 73 insertions(+), 40 deletions(-) diff --git a/fig-stiffnesses.gplot b/fig-stiffnesses.gplot index 4c08763..52762dd 100644 --- a/fig-stiffnesses.gplot +++ b/fig-stiffnesses.gplot @@ -9,11 +9,11 @@ cc5 = "#000000" Tc = 17.26 -a1 = 71.13597475161484 -b1 = 0.010425748328804992 -c1 = 1665.064389018724 -d1 = 6.283065722894796 -e1 = 14.579664883988315 +a1 = 72.9479505025133 +b1 = 0.01227234014115862 +c1 = 2675.3172792640935 +d1 = 8.418493827841177 +e1 = 19.456906054458745 C10(T) = a1 - b1 * T C1(T) = C10(T) / (1 + 1 / (1 / d1 + abs(T - Tc) / c1) / C10(T)) @@ -88,11 +88,11 @@ set title '(c)' offset 12,-2.7 set xlabel '$T / \mathrm K$\\~' offset 0,0.5 set xtics 50, 50, 250 offset 0,0.5 mirror unset y2tics -set ytics 63,1,72 +set ytics 63,1,74 unset y2label set ylabel '$\textcolor{mathc3}{C_{\mathrm{B_{1\mathrm g}}}},\textcolor{mathc4}{C_{\mathrm{B_{1\mathrm g}}}^0} / \mathrm{GPa}$' offset 3 -set yrange [64.5:71.5] +set yrange [64.5:73.5] plot "data/c11mc12.dat" using 1:(100 * $2) with lines lw 6 lc rgb cc3, \ C10(x) dt 3 lw 4 lc rgb cc4 @@ -102,8 +102,8 @@ set y2label '$\Big[C^0_{\mathrm{B_{1\mathrm g}}}(C^0_{\mathrm{B_{1\mathrm g}}}/C set title '(d)' offset 12,-2.7 set format y2 '\scriptsize$%0.2f$' set format y '\scriptsize$%0.2f$' -set yrange [0.14:0.38] -set y2tics 0.15,0.05,0.39 offset -0.7 mirror +set yrange [0.11:0.25] +set y2tics 0.05,0.025,0.25 offset -0.7 mirror plot "data/c11mc12.dat" using 1:(1 / (C10($1)*(C10($1) / (100 * $2) - 1))) with lines lw 6 lc rgb cc3, \ x > Tc ? 1/(C10(x) * (C10(x) / C12(x) - 1)) : 1 / 0 lw 3 lc cc5, \ x <= Tc ? 1/(C10(x) * (C10(x) / C12(x) - 1)) : 1 / 0 dt (2,1) lw 3 lc cc5, \ @@ -127,8 +127,8 @@ plot "data/c11mc12.dat" using 1:(100 * $2) with points pt 2 ps 0.15 lc rgb cc3, unset margins set origin 0.5,0.26 -set yrange [0.154:0.166] -set y2tics 0.155,0.005,0.167 +set yrange [0.1155:0.123] +set y2tics 0.105,0.0025,0.131 set format y2 '\scriptsize$%0.3f$' plot "data/c11mc12.dat" using 1:(1 / (C10($1)*(C10($1) / (100 * $2) - 1))) with points pt 2 ps 0.15 lc rgb cc3, \ diff --git a/fig-stiffnesses.pdf b/fig-stiffnesses.pdf index 4b1d98e..2cc839f 100644 Binary files a/fig-stiffnesses.pdf and b/fig-stiffnesses.pdf differ diff --git a/main.tex b/main.tex index 22b1405..6d133be 100644 --- a/main.tex +++ b/main.tex @@ -476,9 +476,9 @@ corresponding modulus. shows the location of the \ho\ phase. (a) $\Btg$ modulus data and a fit to the standard form.\cite{Varshni_1970} (b) $\Bog$ modulus data and a fit to \eqref{eq:static_modulus} (magenta, dashed) and a fit to \eqref{eq:C0} (black, solid). The fit gives - $C^0_\Bog\simeq\big[71-(0.010\,\K^{-1})T\big]\,\GPa$, $D_\perp - q_*^4/b^2\simeq0.16\,\GPa^{-1}$, and - $a/b^2\simeq6.1\times10^{-4}\,\GPa^{-1}\,\K^{-1}$. Addition of a quadratic + $C^0_\Bog\simeq\big[73-(0.012\,\K^{-1})T\big]\,\GPa$, $D_\perp + q_*^4/b^2\simeq0.12\,\GPa^{-1}$, and + $a/b^2\simeq3.7\times10^{-4}\,\GPa^{-1}\,\K^{-1}$. Addition of a quadratic term in $C^0_\Bog$ was here not needed for the fit.\cite{Varshni_1970} (c) $\Bog$ modulus data and the fit of the \emph{bare} $\Bog$ modulus. (d) $\Bog$ modulus data and the fits transformed by diff --git a/referee_comments.txt b/referee_comments.txt index d4740ff..b10b819 100644 --- a/referee_comments.txt +++ b/referee_comments.txt @@ -17,11 +17,22 @@ > look plausible as the disagreement does not appear to develop slowly > as T decreases, but appears immediately below TN, where eta is small. -The disagreement between the theory at low temperature is largely resolved by -the addition of an additional interaction in the mean-field free energy of the -form ε²η². In a new appendix, we have worked through the mean field modulus -implied with this new interaction and a fit is now shown in Fig. 2 as a thin -black line. +The gross disagreement between the data and theory below the transition is +resolved by the addition of another interaction in the mean-field free energy +of the form ε²η². In a new appendix, we have worked through the mean field +modulus implied with this new interaction and a fit is now shown in Fig. 2 as a +dashed black line. + +Though fine features of the low-temperature behavior are not reproduced, the +addition of further terms in the mean-field free energy supply finer +corrections. For instance, another correction at order η_*² (roughly linear in +ΔT) is produced by a term of the form ε⁴, while myriad terms at yet higher +order provide corrections of order η_*⁴ and up (quadratic and up in ΔT). +Higher order corrections to the mean-field free energy produce arbitrary +analytic corrections to the low-temperature behavior, but do not influence +predictions above the transition. Those can only be further fit by more +complicated dependence of the present free-energy parameters on temperature and +pressure. An inset with a zoom on the critical region has been added to the figure. @@ -35,27 +46,31 @@ temperature window that cuts off at T_max. Our fit's parameters are x_i for i = 1, …, 5 for C₀ = x₁ - x₂ (T / K), x₃ = a/b², x₄ = b²/Dq⁴, and x₅ = b √(-g/u). The variation of these parameters as a function of T_max are shown in referee_response_cutoff_parameters.pdf. The parameter x₁ is fairly stable at -all temperature cutoffs, while the rest vary by at most 20–60% of their 275K -value down to cutoffs of ~90K. +all temperature cutoffs, while the rest vary by 2—5x their 275K value down to +cutoffs of ~90K. More insight into the consistency of the fit comes from examining the linear combinations of parameters that form eigenvectors of the fit covariance matrix, -since these have uncorrelated uncertainties. For the fit including all -temperatures (up to 275K), these are (in order of fit uncertainty): +since—unlike the natural parameters of the mean field theory—these have +uncorrelated uncertainties. For the fit including all temperatures (up to +275K), these are (in order of fit uncertainty): - y₁ = -0.00198126 x₁ + 2.16869 10⁻⁶ x₂ - 0.99998 x₃ - 0.00227756 x₄ - 0.00560291 x₅ - y₂ = -0.0151198 x₁ + 0.0000415145 x₂ - 0.00552438 x₃ - 0.0205384 x₄ + 0.999659 x₅ - y₃ = 0.635138 x₁ - 0.00196902 x₂ - 0.00315925 x₃ + 0.77197 x₄ + 0.0254495 x₅ - y₄ = 0.772222 x₁ - 0.00663886 x₂ - 0.0000753204 x₃ - 0.635317 x₄ - 0.00137316 x₅ - y₅ = 0.00637806 x₁ + 0.999976 x₂ - 4.32279 10⁻⁶ x₃ - 0.00269696 x₄ - 4.93718 10⁻⁷ x₅ + y₁ = -0.0016 x₁ + 1.6 10⁻⁶ x₂ - 1.0 x₃ - 0.0020 x₄ - 0.0044 x₅ + y₂ = -0.017 x₁ + 0.000043 x₂ - 0.0043 x₃ - 0.023 x₄ + 1.0 x₅ + y₃ = -0.62 x₁ + 0.0018 x₂ + 0.0027 x₃ - 0.78 x₄ - 0.029 x₅ + y₄ = -0.78 x₁ + 0.0068 x₂ + 0.000041 x₃ + 0.62 x₄ + 0.0012 x₅ + y₅ = 0.0064 x₁ + 1.0 x₂ - 3.3 10⁻⁶ x₃ - 0.0028 x₄ - 4.4 10⁻⁷ x₅ The variation of these parameter combinations as a function of T_max are shown in referee_response_cutoff_eigenvectors.pdf. The parameter y₁, which is -principally x₃ = a/b², varies the most with the cutoff, at most around 60% of -its value until ~90K. The parameter y₂, which is principally x₅ = b √(-g/u), -varies at most around 15% of its value until ~90K. The other three parameters -are stable at any cutoff, and are likewise mixed combinations of x₁, x₂, and -x₄. +principally x₃ = a/b², varies the most with the cutoff, at most around 5x its +value until ~90K. The parameter y₂, which is principally x₅ = b √(-g/u), varies +at most around 1.5x its value until ~90K. The other three parameters are stable +at any cutoff, and are mixed combinations of x₁, x₂, and x₄. + +Notably, x₁ and x₂ are the only parameters involved in transforming the +experimental data in Figure 2(d), and their stability as a function of the data +window means that transformation is likewise stable. Plots of the fits performed between 90 and 275 K are shown in referee_response_cutoff_curves.pdf. @@ -105,7 +120,9 @@ This is a consistent with hexadecapolar order [hmmmm what can you say to this. i guess we can emphasize the part you did that IS new, at least we think it's new. Also, what's wrong with using existing, working tools for sovling problem?!?!!!!?!?!?! Does everyone have to come up with some black-hole-based-nonsense every time they solve a problem?] -> +In addition, the incorporation of gradient terms into the mean-field free +energy in the context of interpreting ultrasound data appears novel. + > 2) One can generically expect several sources for softening elastic > constants. For example, the authors in ref.25 also succeeded in the > quantitative fits in the framework of a 4f crystal field model for T > @@ -123,13 +140,29 @@ This is a consistent with hexadecapolar order > does not lead uniquely to the authors' arguments on the realization of > the B1g order parameter. -The disagreement between the theory at low temperature is resolved by the -addition of an additional interaction in the mean-field free energy of the form -ε²η², now shown in Fig. 2 as a thin black line. +The gross disagreement between the data and theory below the transition is +resolved by the addition of another interaction in the mean-field free energy +of the form ε²η². In a new appendix, we have worked through the mean field +modulus implied with this new interaction and a fit is now shown in Fig. 2 as a +dashed black line. + +Though fine features of the low-temperature behavior are not reproduced, the +addition of further terms in the mean-field free energy supply finer +corrections. For instance, another correction at order η_*² (roughly linear in +ΔT) is produced by a term of the form ε⁴, while myriad terms at yet higher +order provide corrections of order η_*⁴ and up (quadratic and up in ΔT). +Higher order corrections to the mean-field free energy produce arbitrary +analytic corrections to the low-temperature behavior, but do not influence +predictions above the transition. Those can only be further fit by more +complicated dependence of the present free-energy parameters on temperature and +pressure. -While terms like this provide cusp-like features in the modulus for each strain -symmetry, they cannot explain the 1/ΔT softening seen in the high-temperature -phase, since their contribution to the response function is zero above T_c. +While terms that modify the behavior below the transition provide cusp-like +features in the modulus for each strain symmetry, they cannot explain the 1/ΔT +softening seen in the high-temperature phase, since their contribution to the +response function is zero above T_c. Likewise, mean field theories with a +uniform order parameter cannot explain the finite value of the modulus at the +transition. > 4) The most important point in the HO problem is the microscopic > identification of symmetry breaking and the order parameter. In spite diff --git a/referee_response_cutoff_curves.pdf b/referee_response_cutoff_curves.pdf index f9245c9..8c3e7e5 100644 Binary files a/referee_response_cutoff_curves.pdf and b/referee_response_cutoff_curves.pdf differ diff --git a/referee_response_cutoff_eigenvectors.pdf b/referee_response_cutoff_eigenvectors.pdf index ef69c94..455c060 100644 Binary files a/referee_response_cutoff_eigenvectors.pdf and b/referee_response_cutoff_eigenvectors.pdf differ diff --git a/referee_response_cutoff_parameters.pdf b/referee_response_cutoff_parameters.pdf index 0ed3d76..cee5bc0 100644 Binary files a/referee_response_cutoff_parameters.pdf and b/referee_response_cutoff_parameters.pdf differ -- cgit v1.2.3-70-g09d2