Steady-State Entanglement in the Nuclear Spin Dynamics of a Double Quantum Dot
Rok akademicki: 2021
Academic year: 2021
(2) PHYSICAL REVIEW LETTERS. PRL 111, 246802 (2013) (a). (b). (c). FIG. 1 (color online). (a) Schematic illustration of nuclear entanglement generation via electron transport. Whenever the Pauli blockade is lifted via the HF interaction with the nuclear spins, a nuclear flip can occur in either of the two dots. The local nature of the HF interaction is masked by the nonlocal character of the electronic level j2 i. (b) Spectrum of Hel for ¼ 40 eV and t ¼ 30 eV. The three eigenstates jk i are displayed in red. The triplets jT i are degenerate for !0 ¼ 0. In this setting, lifting of the spin blockade due to HF interaction is predominantly mediated by the nonlocal jump operators required for two-mode squeezing, namely, L2 and L2 . The ellipse refers to a potential operational area of our scheme. (c) The resulting effective three-level system fjT i; j2 ig including coherent HF coupling and the relevant dissipative processes: j2 i decays according to its overlap with jS02 i with an effective decay rate 2 ¼ jh2 jS02 ij2 . Within this three-level subspace, purely electronic Pauli-blockade lifting mechanisms like cotunneling or spin-orbital effects result in effective dephasing and dissipative mixing rates, labeled as deph and , respectively.. X z z aHF X þ þ ðSi Ai þ S Si Ai : (2) i Ai Þ þ aHF 2 i¼1;2 i¼1;2 P Here, Si and Ai ¼ j ai;j i;j for ¼ , z denote electron and collective nuclear spin operators, and ai;j defines the unitless HF coupling constant between the electron spin P i in dot i and the jth nucleus: N j¼1 ai;j ¼ N, where N ¼ 6 ðN1 þ N2 Þ=2 10 refers to the average number of nuclei per dot. The individual nuclear spin operators i;j are assumed to be spin- 12 and we neglect the nuclear Zeeman and dipole-dipole terms . The second term in Eq. (2) can be split into an effective nuclear magnetic P field and residual quantum fluctuations, Hzz ¼ aHF i¼1;2 Szi Azi , where Azi ¼ Azi hAzi it . The (time-dependent) semiclassical OH field exhibits a homogeneous ! OH ¼ ðaHF =2Þ ðhAz1 it þ hAz2 it Þ and inhomogeneous component OH ¼ ðaHF =2ÞðhAz2 it hAz1 it Þ, which can be absorbed into the definitions of !0 and in Eq. (1) as !0 ¼ ! OH þ !ext and ¼ OH þ ext , respectively. For now, we assume HHF ¼. week ending 13 DECEMBER 2013. the symmetric situation of vanishing external fields !ext ¼ ext ¼ 0 . Thus, !0 and are dynamic variables depending on the nuclear polarizations. The flip-flop dynamics, given by the first term in Eq. (2), and the OH fluctuations described by Hzz can be treated perturbatively with respect to the effective electronic Hamiltonian Hel . Its eigenstates within the Sztot ¼Sz1 þSz2 ¼0 subspace can be expressed as jk i ¼ k j*+i þ k j+*i þ k jS02 i (k ¼ 1, 2, 3) with corresponding pﬃﬃﬃﬃ eigenenergies k . For t !0 , gHF , where gHF ¼ N aHF , j1;3 i are far detuned, and the electronic subsystem can be simplified to an effective three-level system comprising the levels fjT i; j2 ig. Effects arising due to the presence of j1;3 i will be discussed below. Within this reduced scheme, Hff reads a Hff ¼ HF ½L2 j2 ihTþ j þ L2 j2 ihT j þ H:c:; (3) 2 þ where the nonlocal nuclear operators L2 ¼ 2 Aþ 1 þ 2 A2 and L2 ¼ 2 A1 þ 2 A2 are associated with lifting the Pauli blockade from jTþ i and jT i via j2 i, respectively. They can be controlled via the external parameters t and defining the amplitudes 2 and 2 . The dynamical evolution of the system is described in terms of a Markovian master equation for the reduced density matrix of the DQD system describing the relevant electronic and nuclear degrees of freedom . Besides the HF dynamics described above, it accounts for other purely electronic mechanisms like, e.g., cotunneling. These effects and their implications for the nuclear dynamics are described in  and lead to effective decay and dephasing processes in the T subspace with rates , deph ; see Fig. 1(c). For fast electronic dynamics ( , deph gHF ) and a sufficiently high gradient * 3 eV (see ), the hybridized electronic level j2 i exhibits a significant overlap with the localized singlet jS02 i and the electronic subsystem settles in the desired quasisteady state, el ss ¼ ðjTþ ihTþ j þ jT ihT jÞ=2, on a time scale much shorter than the nuclear dynamics. One can then adiabatically eliminate all electronic coordinates yielding a coarsegrained equation of motion for the nuclear density matrix ¼ Trel ½ , where Trel ½. . . denotes the trace over the electronic degrees of freedom: _ ¼ Lid ½ þ Lnid ½. Here, the first dominant term describes the desired nuclear squeezing dynamics. Lid ½ ¼ ½D½L2 þ D½L2 þ i ð½Ly2 L2 ; 2 2 y (4) þ ½L2 L2 ; Þ; where D½c ¼ c cy 12 fcy c; g. It arises from coupling to the level j2 i, while Lnid ½ results from coupling to the far detuned levels j1;3 i and OH fluctuations described by Hzz . Here, and refer to a HF-mediated decay rate and Stark shift, respectively, . Pure stationary solutions j ss i associated with the dynamics generated by Eq. (4) can be obtained from the. 246802-2.
(3) PRL 111, 246802 (2013). week ending 13 DECEMBER 2013. PHYSICAL REVIEW LETTERS. dark-state condition L2 j ss i ¼ L2 j ss i ¼ 0. First, we consider the limit of equal dot sizes (N1 ¼ N2 ) and uniform HF coupling (ai;j ¼ N=Ni ) and generalize our results later. The nuclear system can be described via Dicke states jJi ; ki i, where ki ¼ 0; 1; . . . ; 2Ji and Ji refer to the spin-z^ projection and total spin quantum numbers, respectively. For J1 ¼ J2 ¼ J, one readily checks that the dark-state condition is satisfied by the (unnormalized) pure state j ss i ¼ P2J k k¼0 jJ; ki jJ; 2J ki, representing an entangled state closely similar to the two-mode squeezed state . The parameter ¼ 2 =2 quantifies the entanglement and polarization of the nuclear system. j j < 1 (j j > 1) corresponds to states of large positive (negative) OH gradients, respectively. The system is invariant under the symmetry transformation (2 $ 2 , Az1;2 ! Az1;2 ) which gives rise to a bistability in the steady state, as for every solution with positive OH gradient ( > 0), we find another one with < 0. For a given j j 1, the individual nuclear polarizations in the state j ss i approach one as we increase the system size J, and we can describe the system dynamics in the vicinity of the respective steady state in the framework of a HolsteinPrimakoff (HP) transformation . This allows for a detailed analysis of the nuclear dynamics including perturbative effects from the processes P described by Lnid . The collective nuclear spins Ii ¼ j i;j are mapped to bosonic operators  and the (unique) ideal steady state is well known to be a two-mode squeezed state [5,26] which represents j ss i within the HP picture. Since in the bosonic case the modulus of is confined to j j < 1, the HP analysis refers to one of the two symmetric steady-state solutions mentioned above. Within the HP approximation the dynamics generated by _ ¼ Lid ½ þ Lnid ½ are quadratic in the new bosonic creation and annihilation operators. Therefore, the nuclear dynamics are purely Gaussian and exactly solvable. The generation of entanglement can be certified via the EPR entanglement condition [5,30], EPR < 1, where EPR ¼ ½varðI1x þ I2x Þ þ varðI1y þ I2y Þ=ðjhI1z ij þ jhI2z ijÞ. While EPR 1 for separable states, the ideal dynamics Lid drive the nuclear spins into an EPR state with id EPR ¼ ð1 j jÞ=ð1 þ j jÞ < 1. As illustrated in Fig. 2, we numerically find that the generation of steady-state entanglement persists even for asymmetric dot sizes of 20%, classical uncertainty in the total spins Ji , and the undesired terms Lnid . When tuning t from 10 eV to 35 eV, the squeezing parameter j j increases from 0:2 to 0:6, respectively. For j j 0:2, we obtain a relatively high fidelity F with the ideal two-mode squeezed state, close to 80%. For stronger squeezing, the target state becomes more susceptible to the undesired noise terms, first leading to a reduction of F and eventually to a breakdown of the HP approximation. The associated critical behavior can be understood in terms of a dissipative phase transition [28,32]. We now turn to the experimental realization of our scheme : In the analysis above, we discussed the. (a). (b). FIG. 2 (color online). Steady-state entanglement between the two nuclear spin ensembles quantified via (a) the EPR uncertainty EPR and (b) fidelity F of the nuclear steady state with the two-mode squeezed target state. The black solid curve refers to the idealized setting where the undesired HF coupling to j1;3 i has been ignored and where J1 ¼ J2 ¼ pJmax , p ¼ 0:8, and N1 ¼ N2 ¼ 2Jmax ¼ 106 , corresponding to OH ¼ 40 eV. The blue-dashed line then also takes into account coupling to j1;3 i while the red-dashed curve, in addition, accounts for an asymmetric dot size: N2 ¼ 0:8N1 ¼ 8 105 . The amount of entanglement decreases for a smaller nuclear polarization: p ¼ 0:7 (green dashed curve). Classical uncertainty (symbols) in the total spin Ji quantum numbers leads to less entanglement, but does not destroy it.pﬃﬃﬃﬃﬃ Here, we have set the range of the ﬃ distribution to Ji ¼ 50 Ni . Other numerical parameters: !0 ¼ 0, ¼ 25 eV, ¼ 30 eV, and þ deph =2 ¼ 1 eV.. idealized case of uniform HF coupling. However, our scheme also works for nonuniform coupling, provided that the two dots are sufficiently similar: If the coupling is completely inhomogeneous, that is ai;j ai;k for all j k, but the two QDs are identical (a1;j ¼ a2;j 8 j ¼ 1; . . . ; N1 N2 ), Eq. (4) supports a unique pure entangled stationary state. Up to normalization, it reads j ss i ¼ N j¼1 j ij , where j ij ¼ j#j ; "j i þ j"j ; #j i is an entangled state of two nuclear spins belonging to different nuclear ensembles . j ss i features a (large) polarization gradient Iz ¼ hI2z i hI1z i ¼ Nð1 2 Þ=ð1 þ 2 Þ. The buildup of a large OH gradient is corroborated within a semiclassical calculation which neglects correlations among the nuclear spins . This is valid on time scales that are long compared to nuclear dephasing mechanisms [14,34,36]. Assuming equal dot sizes, N1 ¼ N2 ¼ N, we use a semiclassical factorization scheme  resulting in decoupled equations of motion for the two nuclear polarization variables hI1z it and hI2z it . In particular, Iz evolves as. 246802-3. d
(4) Iz ¼ eff Iz N ; dt. eff. (5).
(5) PRL 111, 246802 (2013) (a). PHYSICAL REVIEW LETTERS (b). (c). FIG. 3 (color online). Semiclassical solution to the nuclear polarization dynamics. (a) Instantaneous nuclear polarization rate _ Iz as a function of the gradient for t ¼ 20 eV (blue dashed), t ¼ 30 eV (red solid), and t ¼ 50 eV (black dashdotted). FPs are found at _ Iz ¼ 0. The ovals mark stable highgradient steady-state solutions. The background coloring refers to the sign of _ Iz (for t ¼ 30 eV) which determines the stable FP the nuclear system is attracted to (see arrows). (b) Zoom-in of (a) into the low-gradient regime: The unpolarized FP lies at ¼ 0, whereas critical, instable points crt OH (marked by stars) can be identified with _ Iz ¼ 0 and d_ Iz =d > 0. (c) Stable high-polarization FPs ss OH (see ovals) as a function of t; for t 10 eV we obtain a nuclear polarization of 90%. Other numerical parameters: ¼ 25 eV, ¼ 30 eV, ¼ 0:3 eV, and deph ¼ 0:5 eV.. where the HF-mediated depolarization eff and pumping rate
(6) (see  for their connection to microscopic parameters) depend on the gradient defined in Eq. (1), in particular on the OH gradient OH / Iz . The electronnuclear feedback loop can then be closed self-consistently by identifying steady-state solutions of Eq. (5) in which the parameter is provided by the nuclear OH gradient only. The instantaneous polarization rate _ Iz , given in Eq. (5), is displayed in Fig. 3 as a function of , with the electronic subsystem in its respective steady state, yielding a nonlinear equation for the nuclear equilibrium polarizations. Stable fixed points (FPs) are determined by _ Iz ¼ 0 and d_ Iz =d < 0 as opposed to instable ones where d_ Iz =d > 0 [13,15,38]. We can identify parameter regimes in which the nuclear system features three FPs which are interspersed by two instable points. Two of the stable FPs are high-polarization solutions of opposite sign, supporting a macroscopic OH gradient, while one is the trivial, zero polarization solution. If the initial gradient lies outside the instable points, the system turns self-polarizing and the OH gradient approaches a highly polarized FP. For typical parameter values we estimate that the OH gradient at the instable points is ð1–2Þ eV; compare Fig. 3(b). This comparatively moderate initial gradient could be. week ending 13 DECEMBER 2013. achieved via, e.g., a nanomagnet [39,40] or alternative dynamic nuclear polarization schemes [14,20,35,41]. Next, we address the effects of weak nuclear interactions: First, we have neglected nuclear dipole-dipole interactions. However, we estimate the time scale for the entanglement creation as t ¼ @=N & 10 s which is fast compared to typical nuclear decoherence times, recently measured to be 1 ms in vertical DQDs . Thus, it should be possible to create entanglement between the two nuclear spin ensembles faster than it gets disrupted by dipole-dipole interactions among the nuclei. Second, we have disregarded nuclear Zeeman terms since our scheme requires no external homogeneous magnetic field for sufficiently strong tunneling t . Finally, entanglement could be detected by measuring the OH shift in each dot separately ; in combination with NMR techniques to rotate the nuclear spins  we can obtain all spin components and their variances which are sufficient to verify the presence of entanglement (similar to the proposal ). To conclude, we have presented a scheme for the dissipative entanglement generation among the two nuclear spin ensembles in a DQD. This may provide a long-lived, solid-state entanglement resource and a new route for nuclear-spin-based information storage and manipulation. We acknowledge support by the DFG within SFB 631, the Cluster of Excellence NIM and the project MALICIA within the 7th Framework Programme for Research of the European Commission, under FET-Open Grant No. 265522. E. M. K. acknowledges support by the Harvard Quantum Optics Center and the Institute for Theoretical Atomic and Molecular Physics. L. V. acknowledges support by the Dutch Foundation for Fundamental Research on Matter (FOM)..  B. Kraus and J. I. Cirac, Phys. Rev. Lett. 92, 013602 (2004).  F. Verstraete, M. M. Wolf, and J. I. Cirac, Nat. Phys. 5, 633 (2009).  S. Diehl, A. Micheli, A. Kantian, B. Kraus, H. P. Bu¨chler, and P. Zoller, Nat. Phys. 4, 878 (2008).  R. Sanchez and G. Platero, Phys. Rev. B 87, 081305(R) (2013).  C. A. Muschik, E. S. Polzik, and J. I. Cirac, Phys. Rev. A 83, 052312 (2011); H. Krauter, C. A. Muschik, K. Jensen, W. Wasilewski, J. M. Petersen, J. I. Cirac, and E. S. Polzik, Phys. Rev. Lett. 107, 080503 (2011).  R. Hanson, J. R. Petta, S. Tarucha, and L. M. K. Vandersypen, Rev. Mod. Phys. 79, 1217 (2007).  D. D. Awschalom, N. Smarth, and D. Loss, Semiconductor Spintronics and Quantum Computation (Springer, New York, 2002).  E. A. Chekhovich, M. N. Makhonin, A. I. Tartakovskii, A. Yacoby, H. Bluhm, K. C. Nowack, and L. M. K. Vandersypen, Nat. Mater. 12, 494 (2013).. 246802-4.
(7) PRL 111, 246802 (2013). PHYSICAL REVIEW LETTERS.  M. S. Rudner, L. M. K. Vandersypen, V. Vuletic, and L. S. Levitov, Phys. Rev. Lett. 107, 206806 (2011).  M. S. Rudner, F. H. L. Koppens, J. A. Folk, L. M. K. Vandersypen, and L. S. Levitov, Phys. Rev. B 84, 075339 (2011).  M. J. A. Schuetz, E. M. Kessler, J. I. Cirac, and G. Giedke, Phys. Rev. B 86, 085322 (2012).  K. Ono and S. Tarucha, Phys. Rev. Lett. 92, 256803 (2004).  I. T. Vink, K. C. Nowack, F. H. L. Koppens, J. Danon, Y. V. Nazarov, and L. M. K. Vandersypen, Nat. Phys. 5, 764 (2009).  M. Gullans, J. J. Krich, J. M. Taylor, H. Bluhm, B. I. Halperin, C. M. Marcus, M. Stopa, A. Yacoby, and M. D. Lukin, Phys. Rev. Lett. 104, 226807 (2010).  J. Danon, I. Vink, F. Koppens, K. Nowack, L. Vandersypen, and Yu. Nazarov, Phys. Rev. Lett. 103, 046601 (2009).  A. C. Johnson, J. R. Petta, J. M. Taylor, A. Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, and A. C. Gossard, Nature (London) 435, 925 (2005).  F. H. L. Koppens et al., Science 309, 1346 (2005).  A. V. Khaetskii, D. Loss, and L. Glazman, Phys. Rev. Lett. 88, 186802 (2002).  H. Bluhm, S. Foletti, I. Neder, M. Rudner, D. Mahalu, V. Umansky, and A. Yacoby, Nat. Phys. 7, 109 (2010).  S. Foletti, H. Bluhm, D. Mahalu, V. Umansky, and A. Yacoby, Nat. Phys. 5, 903 (2009).  J. M. Taylor, C. M. Marcus, and M. D. Lukin, Phys. Rev. Lett. 90, 206803 (2003).  W. M. Witzel and S. Das Sarma, Phys. Rev. B 76, 045218 (2007).  H. Ribeiro, J. R. Petta, and G. Burkard, Phys. Rev. B 82, 115445 (2010).  K. Ono, D. G. Austing, Y. Tokura, and S. Tarucha, Science 297, 1313 (2002).  J. Schliemann, A. Khaetskii, and Daniel Loss, J. Phys. Condens. Matter 15, R1809 (2003).  See Supplemental Material at http://link.aps.org/ supplemental/10.1103/PhysRevLett.111.246802 for further details regarding the Master equation, the ideal nuclear steady state, the HP transformation, the semiclassical analysis and a summary of the experimental requirements. 2þ ~  Microscopically, and are given by ¼ a2HF =2½ 2 2 ~ ~ ~ and ¼ ð2 =2Þ , respectively. Here, ¼ 2 þ. =2 þ deph =4.. week ending 13 DECEMBER 2013.  E. M. Kessler, G. Giedke, A. Imamoglu, S. F. Yelin, M. D. Lukin, and J. I. Cirac, Phys. Rev. A 86, 012116 (2012).  Here, we consider the subspace with large collective spin quantum numbers, Ji OðN=2Þ. The zeroth-order HP mapping can be justified self-consistently, provided that the occupations in the bosonic modes are small compared to 2Ji .  M. G. Raymer, A. C. Funk, B. C. Sanders, and H. deGuise, Phys. Rev. A 67, 052104 (2003).  We average over an uniform distribution of fJ1 ; J2 g subpﬃﬃﬃﬃﬃﬃ spaces with a range of Ji ¼ 50 Ni . The center of the distribution Ji has been taken as Ji ¼ pNi =2, where the polarization p is set by the OH gradient via p ¼ max OH =max OH ; here, OH ¼ AHF =2 50 eV.  M. J. A. Schuetz et al. (unpublished).  Numerical evidence (for small systems) indicates that small deviations from perfect symmetry between the QDs still yield an entangled (mixed) steady state close to j ss i .  We estimate 1 eff 1 s; this is compatible with the semiclassical approximation and in agreement with typical polarization time scales [14,35].  R. Takahashi, K. Kono, S. Tarucha, and K. Ono, Phys. Rev. Lett. 107, 026602 (2011).  H. Christ, J. I. Cirac, and G. Giedke, Phys. Rev. B 75, 155324 (2007).  The results obtained within this approximative factorization scheme have been confirmed by numerical simulations for small sets of nuclei .  H. Bluhm, S. Foletti, D. Mahalu, V. Umansky, and A. Yacoby, Phys. Rev. Lett. 105, 216803 (2010).  M. Pioro-Ladrie`re, T. Obata, Y. Tokura, Y.-S. Shin, T. Kubo, K. Yoshida, T. Taniyama, and S. Tarucha, Nat. Phys. 4, 776 (2008).  G. Petersen, E. A. Hoffmann, D. Schuh, W. Wegscheider, G. Giedke, and S. Ludwig, Phys. Rev. Lett. 110, 177602 (2013).  J. Petta, J. Taylor, A. Johnson, A. Yacoby, M. Lukin, C. Marcus, M. Hanson, and A. Gossard, Phys. Rev. Lett. 100, 067601 (2008).  Note that any initial OH splitting ! OH is damped to zero in the steady state .  For fast recharging of the DQD, ¼ R =2, where R is the sequential tunneling rate to the right lead .. 246802-5.
Warto podkreślić, że czas odgrywa istotną rolę w procesie leczenia skutków udaru mózgu i ma bardzo duże znaczenie na etapie od wystąpienia jego objawów do pod- jęcia
Nauczanie Piusa XII na temat podmiotów realizujących wycho wawcze posłannictwo Kościoła jest oparte na teologii przedsoboro- wej uznającej, że w Kościele jest jeden
- (ustęp 2) „W zakresie, który nie podlega jej wyłącznej kompetencji, Wspólnota po dejmuje działania, zgodnie z zasadą subsydiamości, tylko wówczas i tylko w
congruity between being and form, ontology and aesthetics, and how the subject, being in dialogue with itself and an other, may conceive of itself as form; my interpretation will
Jasno rysuje nam się teraz odpowiedź na postawione na początku szkicu pytanie o to, jakie tłumione wartości duchowe społeczeństwa polskiego u progu
The KM value for free formaldehyde 21 lM with ferredoxin as an electron acceptor is approximately 3 times lower than the value measured when benzyl viologen is used as an acceptor..
Here, we demonstrate the heralded initialisation, control and entanglement of individual dark spins associated to multiple P1 centers, which are part of a spin bath surrounding
Our analysis is based on the typical separation of time scales between (fast) electron spin evolution and (slow) nuclear spin dynamics, yielding a coarse-grained quantum master
Poczynając od maja 1917 г., do powołanej przy Uniwersytecie w Rostowie Komisji Likwi- dacyjnej Uniwersytetu Warszawskiego zaczęły wpływać pisma od po- szczególnych
Kultura i Polityka : zeszyty naukowe Wyższej Szkoły Europejskiej im... 1997 .A ll President’s Words: The Bully Pulpit and the Creation
In this paper the case of the rotating ring under a stationary constant load and that of a stationary ring subjected to a moving load are compared in terms of their resonance speeds,
A signiﬁcant variable inﬂuencing trust between coopera!ng ﬁrms that focus on collabora!ve learning is convergence over their strategies (Valkokari and Helander, 2007).
Therefore, the present study was performed in intention to analyze the reliability of maximal lactate production rate ( VLa · max ) by using a sprint test, maximum oxygen
Z połączenia dostępnych materiałów londyńskich i nowojorskich powstaje zbiór 237 listów wymienianych przez Lechonia i Wierzyńskiego w latach 1941–1956, przy czym Lechoń
We show that this results in tuning as well as detuning of the resonance by DNSP, depending on the mismatch be- tween the driving frequency and g B jB 0 j= @: At small mismatch,
produktie transport GROEVE transport laden HAVEN transport lossen HAVEN transport •transport ^^6 it verwerken WERK. figuur 4 Transport
tej. Nadanie przez monarchę indygenatu obcemu szlachcicowi było równoznaczne z uzyskaniem szlacheckiego obywatelstwa. Przysługiwał osobie, która go uzyskała, ale też jej
Forced vibrations are presented in terms of maximal amplitudes of the steady-state torsional stresses in the thrust shaft, intermediate shaft and propeller shaft of the actual
substrate binding site, as was manifested by the higher apparent K m value at 1 M
Doyon, Entanglement negativity and entropy in non-equilibrium conformal field theory, Nucl.. Schalm, Energy flow in quantum critical systems far from equilibrium,
In section 3 and 4 we explain the two complementary numerical methods, shooting and matching, in detail and present the consistent results on the time evolution of the
Recent mathematical, theoretical and numerical studies have introduced new concepts on peritoneal transport and yielded a better description of particular processes such as pure