Pith. sign in

REVIEW 4 major objections 5 minor 43 references

A Morphologically Self-Consistent Phase Field Model for the Computational Study of Memristive Thin Film Current-Voltage Hysteresis

T0 review · 4 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read A phase field model predicts memristor filament growth and current-voltage hysteresis without assuming any filament geometry.

desk verdict A promising phase field memristor model undermined by a variational error and hand-tuned parameters; fixable, but the thermodynamic claim is not supported as written. read the letter →

arxiv 2506.17421 v1 pith:XH7VQWRJ submitted 2025-06-20 cond-mat.mtrl-sci physics.app-ph

John F. Sevic , Ambroise Juston , Nobuhiko P. Kobayashi This is my paper
classification cond-mat.mtrl-sciphysics.app-ph
keywords phasefieldmodelmemristivethinfilmconductingfilamentcurrent-voltagehysteresisoxygenvacancyLandauorderparameterCahn-Hilliardequationstochasticlocaldensityofstates
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper proposes a phase field model of memristive thin films in which conducting filaments emerge spontaneously from the stochastic distribution of atomic-scale defects, with no prescribed filament shape. The model couples a Landau order parameter that represents local conductivity to the electrostatic potential, and it reproduces the characteristic pinched current-voltage hysteresis along with a roughly tenfold resistance swing. If correct, it offers a computationally cheap way to predict how fabrication randomness translates into device-to-device variability, uniformity, and endurance across a wafer. The key claim is that the thermodynamic path a filament takes is selected by the random local density of states of the as-fabricated film.

What carries the argument

The order parameter φ(r), with 0 ≤ φ ≤ 1, encodes local electrical conductivity, with two wells in the free energy density corresponding to the insulating and conducting states. The stochastic local density of states enters through Equation 3, which maps the physical disorder to the order parameter domain using two empirical constants, η and φ₀. The Cahn-Hilliard equation (Equation 6) drives φ toward a free-energy minimum while the Poisson equation (Equation 7) enforces charge conservation, and the two are solved self-consistently at each bias step.

What would settle it

Run the same simulation with a perfectly uniform initial conductivity (no random LDOS variations) and observe whether a conducting filament still forms; the model's logic implies no filament should appear in the symmetric case, mirroring its earlier flat-interface result. Alternatively, measure the actual as-fabricated defect distribution of a device, feed it as the initial condition, and compare the predicted filament path to in-situ TEM images.

Watch

Extended reading notes

Core claim

The central discovery is that current-voltage hysteresis and the morphology of conducting filaments can be computed self-consistently from a single order parameter field that tracks the local conductive state, without imposing a filament geometry. Minimizing the Gibbs free energy through a Cahn-Hilliard equation coupled to Poisson's equation yields filaments that grow along paths favored by the random structural and chemical inhomogeneities of the film. This reproduces the expected even-symmetric hysteresis, a maximum high-to-low resistance ratio of about 10 near ±25 mV, and filament shapes consistent with experimental imaging.

Load-bearing premise

The mapping in Equation 3 from the physical local density of states to the order parameter domain relies on two empirical constants, η and φ₀, which the authors set to 1000 and 0.8 to produce the hysteresis they want; if that mapping is not physically grounded, the predicted filament paths and hysteresis are artifacts of parameter choice.

Editorial extensions

If this is right

  • Filament paths and hysteresis become predictions from the as-fabricated defect map rather than inputs, so the same model can be swept over many random realizations to produce wafer-scale statistics.
  • Structurally symmetric films are predicted to show symmetric current-voltage response; any asymmetry in a device must come from symmetry breaking in its initial state.
  • The model is calibrated for tantalum oxide valence-change devices but the parameter set can be re-fit to other transition metal oxides, giving a general computational design tool.
  • Because the filament forms spontaneously, the method can be used to study the statistics of filament nucleation and rupture under repeated cycling, supporting endurance analysis.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The empirical constants η and φ₀, hand-set to 1000 and 0.8, could in principle be tied to physically measurable quantities such as defect density and activation energy; if a calibration relation exists, the model would gain predictive power beyond the fitting set.
  • The same variational machinery could be extended beyond isothermal valence-change switching to electrothermal and phase-change memristors by coupling the order parameter to a temperature field.
  • A direct test would be to seed the model with the measured defect distribution of a real device, run the simulation, and compare the predicted filament position to in-situ TEM images of the same device.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The manuscript proposes a multiphysics phase field model for memristive thin films, combining a Cahn-Hilliard equation for an order parameter φ(r) with a Poisson equation for the electrostatic potential. The charge density is mapped from a stochastic local density of states through n ≈ exp[η(φ−φo)], and the coupled equations are solved with MOOSE for logarithmic voltage sweeps. The reported outputs are current-voltage hysteresis loops and spontaneously forming conducting filaments whose morphology is compared qualitatively with experimental TEM images. The paper claims that the method is morphologically self-consistent in that no a priori filament geometry is prescribed and that filaments evolve on energetically favored thermodynamic paths determined by stochastic atomic-scale variations.

Significance. The central idea—using a variational phase field formulation so that filament morphology and hysteresis emerge from random as-fabricated disorder without a prescribed geometry—is attractive and potentially useful for wafer-scale memristor screening. The model is compact and computationally efficient, and it produces plausible qualitative features: even-symmetric hysteresis for symmetric initial conditions, high/low resistance states, and filamentary structures. However, for the contribution to be significant the dynamics must actually be the gradient flow of the stated free energy, the electrostatic coupling must be correct, and the parameters that control the coupling to the random LDOS must be grounded or demonstrated to be robust. As written, these load-bearing requirements are not met; the paper is a promising demonstration rather than an established predictive method.

major comments (4)
  1. [Section III, Eqs. (5) and (6)] The variational derivation is inconsistent as printed. Equation (5) contains the interface energy density (κ/2)∇²φ. The first variation of ∫(κ/2)∇²φ d³r is, up to boundary terms that vanish under the stated boundary conditions, zero; it is not a bulk gradient-energy density. Consequently, the functional derivative of Eq. (5) cannot produce the −κ∇²φ term in the chemical potential in Eq. (6). That term is the functional derivative of (κ/2)|∇φ|². Equation (6) is therefore the Cahn-Hilliard gradient flow of a different free energy than the one stated in Eq. (5). Additionally, the electrostatic energy in Eq. (5) is φVq, whose variational derivative is Vq, not V as written in Eq. (6). Because the central claim in the abstract and Section IV is that filaments evolve on thermodynamic paths favored by the Gibbs free energy, this inconsistency must be fixed, for example by replacing (κ/2)∇²φ with (κ/2)|∇φ|² and re-deriving Eq. (6).
  2. [Section III, Eq. (7)] Equation (7) is written as ∇²V = n/(εr ε0), but with the standard electrostatic relation E = −∇V used later in Section III, the Poisson equation should be ∇²V = −n/(εr ε0) for positive charge density n. With the printed sign, the electrostatic forcing of the phase field has the wrong polarity relative to the stated potential convention. This affects the predicted filament location and the sign structure of the I-V curves and should be corrected, or the sign convention explicitly redefined.
  3. [Section III, Eq. (3), Table I] The mapping from the LDOS to the order parameter domain uses empirical constants η=1000 and φo=0.8, and the text states these values were chosen as 'producing optimum current-voltage hysteresis.' Because Eq. (3) is the sole physical coupling between φ and the charge density, the observed hysteresis and the selected filament paths may be generated by these hand-set parameters rather than by the stochastic LDOS. The paper should provide an independent calibration procedure for η and φo, or bounds derived from LDOS or transport data, together with a sensitivity analysis over their ranges. Without this, the claim that the model 'correctly predicts' filament evolution is not established.
  4. [Section IV, Figs. 3 and 4] The agreement with experimental imaging is only qualitative: no quantitative metric, no direct side-by-side comparison, and no statistics over realizations of the random LDOS are reported. Since the LDOS is defined as a random process in Section II and the model is stochastic, a single realization does not support the general claims in the abstract about stochastic structural and chemical variations. I recommend reporting an ensemble of simulations with different G(r, φ, t0) realizations and comparing filament diameters, positions, and switching voltages to experimental values.
minor comments (5)
  1. [Section II] The phrase 'the variable to indicates' should read 'the variable t0 indicates'.
  2. [Section III and Fig. 1] The text gives the active region as 'x1 = 5 nm and x1 = 35 nm'; the second coordinate should be x2, and the same correction is needed in the definition of the current integral in Eq. (8).
  3. [Section III, Eq. (4) and footnote [36]] The explicit polynomial forms of a2(φ), a4(φ), and a6(φ) appear only in a footnote; they should be stated in the main text, along with a clear statement of the units of fbulk and of the coefficients.
  4. [Section IV] The sentence describing point A of Fig. 4 contains a typo: 'substantial a increase' should be 'substantial an increase' or 'a substantial increase'.
  5. [Section III, Eq. (8)] The current is evaluated using charge density and electric field only at the top electrode edge; the manuscript should justify why this boundary evaluation is equivalent to the total current, given that no drift-diffusion continuity equation is solved.

Circularity Check

2 steps flagged · score 7.0 of 10

IV hysteresis and filament morphology are hand-tuned outputs: η and φo are set to 'producing optimum current-voltage hysteresis' and κ is set to enforce experimental morphology, so the 'predictions' reduce by construction to fitted parameters.

  1. fitted input called prediction [Section III, Eq. (3), Table I]
    "Since our model is calibrated for transition metal oxides of Type II electric conduction, we therefore anticipate a nonlinear current-voltage response for the positive-going bias and an approximately ohmic response for the negative-going bias. Based on this physical reasoning, emission parameter and mean order parameter η and φo have been set to 1000 and 0.8, respectively, producing optimum current-voltage hysteresis."

    Equation (3) defines the charge density in terms of empirical constants η and φo: n(r,φ) ≈ Σ G(r,φ,t0) exp[η(φ−φo)]. The text then states that η and φo were set to 1000 and 0.8 specifically to produce optimum current-voltage hysteresis. The current-voltage hysteresis is the paper's central claimed result, so the simulated loop is not a prediction from first principles; it is an output manufactured by tuning the mapping constants to the target phenomenon. The positive-going nonlinearity and negative-going ohmic response used as 'physical reasoning' are the same features that the resulting hysteresis exhibits.

  2. fitted input called prediction [Section III, definition of κ; Section IV, morphology claim]
    "The interface energy density, κ, is established to enforce conducting filament morphology consistent with experimental data [38]."

    The parameter κ is not derived from an independent measurement; it is set so that the model's filament morphology matches experimental data. Later the paper claims: 'Conducting filament morphology of our model is consistent with experimental imaging data [31][32][33].' Because the value of κ was chosen to enforce that consistency, the agreement between simulated and observed morphology is built into the model rather than being an independent confirmation. The morphology claim therefore reduces to a restatement of the parameter choice.

full rationale

The most distinctive outputs of the paper—the current-voltage hysteresis loop and the conducting-filament morphology—are not independent predictions. The exponential mapping constants η and φo are explicitly hand-set to 'producing optimum current-voltage hysteresis,' and κ is explicitly chosen to enforce experimentally consistent filament morphology. Thus the headline agreement with target behavior is guaranteed by construction for these features. The paper also states that the model minimizes a Gibbs free energy and then claims it 'correctly predicts conducting filaments evolve on thermodynamic paths that are energetically favored'; to the extent the dynamics are a gradient-flow minimization, this restates the model's design rather than testing it. The self-citations to the authors' own earlier work ([14], [15], [28]) supply the double-well free energy and transport parameters, but by themselves self-citation is not circular; the present computation does solve the coupled Cahn-Hilliard/Poisson system with a random initial LDOS and can produce emergent spatial structure. A separate mathematical inconsistency—Equation (5) writes the interface term as (κ/2)∇²φ, whose bulk variation vanishes, while Equation (6) contains the −κ∇²φ term appropriate to (κ/2)|∇φ|²—further weakens the thermodynamic-path claim, but that is a correctness issue rather than a circularity. Overall, because the central observable is tuned rather than predicted, the circularity score is 7.

Assumptions & free parameters 6 free parameters · 4 assumptions · 0 invented entities

The model depends on several free parameters, most prominently η and φo which are tuned to produce the hysteresis. The free energy coefficients and mobilities are imported from the authors' earlier works. No new physical entities are introduced.

free parameters (6)
  • η (emission parameter) = 1000
    Set by hand to produce 'optimum current-voltage hysteresis' (Section III); no physical derivation.
  • φo (mean order parameter) = 0.8
    Same as η, chosen to produce optimum hysteresis.
  • κ (interface energy density) = 1.0 eV/nm²
    Taken from authors' prior model [14]; chosen to 'enforce conducting filament morphology consistent with experimental data'.
  • M (phase field mobility) = 100 nm/(J·ns)
    Taken from [14]; sets charge carrier relaxation time, not independently measured.
  • μo (electrical mobility) = 100 nm/(V·ns)
    Taken from [15]; assumed uniform.
  • f_bulk coefficients a2, a4, a6 = -0.065, 0.70, 0.25
    Double-well free energy coefficients taken from authors' earlier electrothermal phase field model at 300 K (footnote [36]).
assumptions (4)
  • domain assumption Charge carriers obey a Boltzmann energy distribution (f ≈ exp[-βE]).
    Assumed in Section II; requires weakly interacting charges, reasonable but unverified for filament regime.
  • domain assumption States are sufficiently delocalized that long-range correlation can be ignored.
    Stated in Section II as a 'reasonable approximation for amorphous thin films'; underpins the mean-field treatment.
  • domain assumption Landau mean field theory with a scalar order parameter φ describes the system.
    Adopted in Section II; no derivation of validity for oxygen vacancy transport.
  • ad hoc to paper Electrostatic coupling is gelec(φ,V)=φ V q.
    Introduced in Section III; a simplified linear coupling with no derivation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A Morphologically Self-Consistent Phase Field Model for the Computational Study of Memristive Thin Film Current-Voltage Hysteresis." pith.science (2026). https://pith.science/paper/XH7VQWRJ

@misc{pith2026250617421,
  author       = {Pith},
  title        = {Pith review of: A Morphologically Self-Consistent Phase Field Model for the Computational Study of Memristive Thin Film Current-Voltage Hysteresis},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XH7VQWRJ}},
  note         = {Machine review of arXiv:2506.17421}
}
read the original abstract

A multiphysics phase field model is used for the computational study of memristive thin film morphology and current-voltage hysteresis. In contrast to previous computational methods, no requirements are made on conducting filament geometry. Our method correctly predicts conducting filaments evolve on thermodynamic paths that are energetically favored due to stochastic structural and chemical variations naturally occurring at the atomic-level, due to both latent and intentional fabrication effects. These results have significant implications for the computational design of a broad class of memristive thin films, enabling practical wafer-scale mapping, uniformity, and endurance analysis and optimization.

Figures

Figures reproduced from arXiv: 2506.17421 by the authors.

Figure 2
Figure 2. FIG. 2. Free energy density as a function order parameter [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 1
Figure 1. FIG. 1. Cross section and initial normalized electrical con [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 1
Figure 1. To evaluate the full range of performance of [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figures from the paper (2 more)
Figure 3
Figure 3. Figure 3: FIG. 3. The computational [PITH_FULL_IMAGE:figures/full_fig_p005_3.png]
Figure 3
Figure 3. Figure 3: V. SUMMARY We have presented a multiphysics phase field model inherently compatible with randomness associated with as-fabricated memristive thin films. By treating the mo￾tion of oxygen vacancies and their interactions with ran￾dom states of these amorphous thin films…

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

43 extracted references · 42 canonical work pages

  1. [1]

    Mott and R

    N. Mott and R. Gurney, Electronic Processes and Ionic Crystals, 1st ed. (Oxford University Press, London, Eng- land, 1940)

    work page 1940

  2. [2]

    T. W. Hickmott, Low-frequency negative resistance in thin anodic oxide films, Journal of Applied Physics 33, 2669 (1962)

    work page 1962

  3. [3]

    N. F. Mott, Metal-insulator transition, Rev. Mod. Phys. 40, 677 (1968)

    1968

  4. [4]

    Waser and M

    R. Waser and M. Aono, Nanoionics-based resistive switching memories, Nature materials 6, 833 (2007)

    work page 2007

  5. [5]

    D. B. Strukov, G. S. Snider, D. R. Stewart, and R. S. Williams, The missing memristor found, Nature 453, 80 EP (2008)

    work page 2008

  6. [6]

    Gibson, S

    G. Gibson, S. Musunuru, J. Zhang, K. Vandenberghe, J. Lee, C. Hsieh, W. Jackson, Y. Jeon, D. Henze, Z. Li, and R. Williams, An accurate locally active memristor model for s-type negative differential resistance in nbox, Applied Physics Letters 108, 023505 (2016)

    work page 2016

  7. [7]

    Nardi, S

    F. Nardi, S. Larentis, S. Balatti, D. Gilmer, and D. Ielmini, Resistive switching by voltage-driven ion mi- gration in bipolar rram—part i: Experimental study, IEEE Trans. Electron Devices Lett. 59, 2461 (2012)

    work page 2012

  8. [8]

    Larentis,

    S. Larentis, . F. Nardi, , S. Balatti, D. Gilmer, and D. Ielmini, Resistive switching by voltage-driven ion mi- gration in bipolar rram - part ii: Modeling, IEEE Trans. Electron Devices Lett. 59, 2468 (2012)

    work page 2012

Show all 43 references
  1. [9]

    S. Kim, S. Kim, K. M. Kim, S. R. Lee, M. Chang, E. Cho, Y.-B. Kim, C. J. Kim, U. In Chung, and I. K. Yoo, Physical electrothermal model of resistive switching in bi-layered resistance-change memory, Scientific Reports 3, 1680 EP (2013)

    2013

  2. [10]

    J. F. Sevic and N. P. Kobayashi, Self-consistent continuum-based transient simulation of electrofor- mation of niobium oxide tantalum dioxide selector- memristor structures, Journal of Applied Physics 124, 164501 (2018)

    2018

  3. [11]

    W. Shen, N. Kumari, G. Gibson, Y. Jeon, D. Henze, S. Silverthorn, C. Bash, and S. Kumar, Effect of anneal- ing on structural changes and oxygen diffusion in amor- phous hfo2 using classical molecular dynamics, Journal of Applied Physics 123, 085113 (2018)

    2018

  4. [12]

    Fyrigos, V

    I.-A. Fyrigos, V. Ntinas, G. C. Sirakoulis, P. Dimitrakis, and I. G. Karafyllidis, Quantum mechanical model for filament formation in metal-insulator-metal memristors, IEEE Transactions on Nanotechnology 20, 113 (2021)

    2021

  5. [13]

    Funck and S

    C. Funck and S. Menzel, Comprehensive model of elec- tron conduction in oxide-based memristive devices, ACS Applied Electronic Materials 3, 3674 (2021)

    2021

  6. [14]

    J. F. Sevic and N. P. Kobayashi, A computational phase field study of conducting channel formation in dielectric thin films: A view toward the physical origins of resis- tive switching, Journal of Applied Physics 126, 065305 (2019)

    2019

  7. [15]

    J. F. Sevic and N. P. Kobayashi, Resistive switching con- ducting filament electroformation with an electrothermal phase field method, Applied Physics Letters 123 (2023)

    2023

  8. [16]

    F. Pan, C. Chen, Z. shun Wang, Y. chao Yang, J. Yang, and F. Zeng, Nonvolatile resistive switching memories- characteristics, mechanisms and challenges, Progress in Natural Science: Materials International 20, 1 (2010)

    2010

  9. [17]

    F. Miao, J. P. Strachan, J. J. Yang, M.-X. Zhang, I. Goldfarb, A. C. Torrezan, P. Eschbach, R. D. Kelley, G. Medeiros-Ribeiro, and R. S. Williams, Anatomy of a nanoscale conduction channel reveals the mechanism of a high performance memristor, Advanced Materials 23, 5633 (2011)

    2011

  10. [18]

    J. P. Strachan, M. D. Pickett, J. J. Yang, S. Aloni, A. L. David Kilcoyne, G. Medeiros-Ribeiro, and R. Stan- ley Williams, Direct identification of the conducting channels in a functioning memristive device, Advanced Materials 22, 3573 (2010)

    2010

  11. [19]

    N. Xu, L. Liu, X. Sun, X. Liu, D. Han, Y. Wang, R. Han, J. Kang, and B. Yu, Characteristics and mechanism of conduction and set process in tin-zno-pt resistance switching random-access memories, Applied Physics Let- ters 92, 232112 (2008)

    2008

  12. [20]

    By hopping we mean both discrete jumps in energy and space

  13. [21]

    A. H. Wilson, A note on the theory of rectification, Pro- ceedings of the Royal Society of London 136, 487 (1932). 7

    1932

  14. [22]

    N. F. Mott, Note on the contact between a metal and an insulator or semiconductor, Proceedings of Cambridge Philosophical Society 34, 568 (1938)

    1938

  15. [23]

    S. D. Baranovskii, Theoretical description of charge transport in disordered organic semiconductors, physica status solidi (b) 251, 487 (2014)

    2014

  16. [24]

    Zhang, Evolution of the conductive filament system in hfo2 based memristors observed by direct atomic-scale imaging, Nat Commun 12, 7232 (2021)

    Y. Zhang, Evolution of the conductive filament system in hfo2 based memristors observed by direct atomic-scale imaging, Nat Commun 12, 7232 (2021)

    2021

  17. [25]

    Diaz-Leon, K

    J. Diaz-Leon, K. Norris, J. Sevic, J. Yang, and N. Kobayashi, Integration of a niobium oxide selector on a tantalum oxide memristor by local oxidation using joule heating, SPIE, San Diego (2016)

    2016

  18. [26]

    J. J. Yang, F. Miao, M. D. Pickett, D. A. A. Ohlberg, D. R. Stewart, C. N. Lau, and R. S. Williams, The mechanism of electroforming of metal oxide memristive switches, Nanotechnology 20, 215201 (2009)

    2009

  19. [27]

    K. M. Kim, T. H. Park, and C. S. Hwang, Dual conical conducting filament model in resistance switching tio2 thin films, Scientific Reports 5, 7844 EP (2015)

    2015

  20. [28]

    Diaz-Leon, K

    J. Diaz-Leon, K. Norris, J. Yang, J. Sevic, and N. Kobayashi, A niobium oxide-tantalum oxide selector- memristor self-aligned nanostack, Appl. Phys. Lett. 110, 103102 (2017)

    2017

  21. [29]

    M. Kim, M. A. Rehman, D. Lee, Y. Wang, D.-H. Lim, M. F. Khan, H. Choi, Q. Y. Shao, J. Suh, H.-S. Lee, and H.-H. Park, Filamentary and interface-type memristors based on tantalum oxide for energy-efficient neuromor- phic hardware, ACS Applied Materials & Interfaces 14, 44561 (2022)

    2022

  22. [30]

    Yadav, A

    D. Yadav, A. K. Dwivedi, S. Verma, and D. K. Avasthi, Transition metal oxide based resistive random-access memory: An overview of materials and device perfor- mance enhancement techniques, Journal of Science: Ad- vanced Materials and Devices 9, 100813 (2024)

    2024

  23. [31]

    Y. Yang, P. Gao, S. Gaba, T. Chang, X. Pan, and W. Lu, Observation of conducting filament growth in nanoscale resistive memories, Nature communications 3, 732 (2012)

    2012

  24. [32]

    Ahmed, S

    T. Ahmed, S. Walia, E. L. Mayes, R. Ramanathan, P. Guagliardo, V. Bansal, M. Bhaskaran, J. J. Yang, and S. Sriram, Inducing tunable switching behavior in a sin- gle memristor, Applied Materials Today 11, 280 (2018)

    2018

  25. [33]

    Bejtka, G

    K. Bejtka, G. Milano, C. Ricciardi, C. F. Pirri, and S. Porro, Tem nanostructural investigation of ag- conductive filaments in polycrystalline zno-based resis- tive switching devices, ACS Applied Materials & Inter- faces 12, 29451 (2020)

    2020

  26. [34]

    L. D. Landau and E. M. Lifshitz, Statistical Physics, 3rd ed. (Pergamon Press, Elmsford, New York, 1980)

    1980

  27. [35]

    Provatas and K

    N. Provatas and K. Elder, Phase-Field Methods in Mate- rials Science and Engineering (Wiley-VCH Verlag, Wein- heim, Germany, 2010)

    2010

  28. [36]

    For our model the following appropriately scaled even- order polynomial functions in order parameter φ are used. These functions are based on our earlier electrother- mal phase field model but held constant at 300 K abso- lute temperature for the present isothermal treatment of...

  29. [37]

    J. W. Cahn and J. E. Hilliard, Free energy of a nonuni- form system. i. interfacial free energy, The Journal of Chemical Physics 28, 258 (1958)

    1958

  30. [38]

    S. G. Kim, W. T. Kim, and T. Suzuki, Phase-field model for binary alloys, Phys. Rev. E 60, 7186 (1999)

    1999

  31. [39]

    Gaston, C

    D. Gaston, C. Newman, G. Hansen, and D. Lebrun- Grandi´ e, Moose: A parallel computational framework for coupled systems of nonlinear equations, Nuclear Engi- neering and Design 239, 1768 (2009)

    2009

  32. [40]

    M. R. Tonks, D. Gaston, P. C. Millett, D. Andrs, and P. Talbot, An object-oriented finite element framework for multiphysics phase field simulations, Computational Materials Science 1, 20 (2012)

    2012

  33. [41]

    Balay, S

    S. Balay, S. A. andMark F. Adams, J. Brown, P. Brune, K. Buschelman, L. Dalcin, V. Eijkhout, W. D. Gropp, D. Kaushik, M. G. Knepley, L. C. McInnes, K. Rupp, B. F. Smith, S. Zampini, H. Zhang, and H. Zhang, PETSc Users Manual , Tech. Rep. ANL-95/11 - Revision 3.7 (Ar- gonne Nat...

    2016

  34. [42]

    B. S. Kirk, J. W. Peterson, R. H. Stogner, and G. F. Carey, libMesh: A C++ Library for Parallel Adaptive Mesh Refinement/Coarsening Simulations, Engineering with Computers 22, 237 (2006)

    2006

  35. [43]

    See specifically panels (a) and (b) Figure 4 of [14]

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.