The Breakdown of Superconductivity Due to Strong Fields for the Ginzburg Landau Model

Size: px
Start display at page:

Download "The Breakdown of Superconductivity Due to Strong Fields for the Ginzburg Landau Model"

Transcription

1 SIAM REVIEW Vol. 44, No., pp c 00 Socety for Industral and Appled Mathematcs The Breakdown of Superconductvty ue to Strong Felds for the Gnzburg Landau Model T. Gorg. Phllps Abstract. We study the behavor of a superconductng materal subjected to a constant appled magnetc feld, H a = he wth e =, usng the Gnzburg Landau theory.we analytcally show the exstence of a crtcal feld h, for whch, when h>h, the normal states are the only solutons to the Gnzburg Landau equatons.we estmate h.as κ 0 we derve h = O(), whle as κ we obtan h = O(κ). Key words. superconductvty, Gnzburg Landau equatons, upper crtcal felds, normal state AMS subject classfcatons. 35J60, 35J65, 35Q40 PII. S Introducton. If a superconductng body s subjected to a suffcently strong appled magnetc feld, ts ablty to act as a superconductor breaks down and only the normally conductng (resstve) state s observed. In ths paper we consder superconductvty as modeled by the Gnzburg Landau theory and establsh ths type of phenomenon. Here, superconductvty s characterzed n terms of a complex valued order parameter ψ (where ψ represents the densty of superconductng electron pars) and a vector feld A the magnetc potental. Consder a superconductng body gven by a bounded doman R n, where n = or 3 and s of class C,α for some 0 <α<. Assume the body has constant permeablty normalzed to and that the exteror conssts of a second materal wth constant permeablty µ e > 0. efne the permeablty densty as follows: µ(x) = forx, = µ e for x R n \. A magnetc feld s appled to all space n the form H a = he, where h s a postve constant and e R 3 s a fxed unt vector. The presence of produces an Publshed electroncally May, 00.Ths paper orgnally appeared n SIAM Journal on Mathematcal Analyss, Volume 30, Number, 999, pages Ths research was sponsored n part by NSF grant MS and also supported by the epartment of Energy grant E-FG0-90ER4547 through the Mdwest Superconductvty Consortum. Mathematcs epartment, Towson Unversty, 8000 York Road, Towson, M (tgorg@towson.edu). 395 epartment of Mathematcs, Purdue Unversty, West Lafayette, IN (phllps@ math.purdue.edu). 37

2 38 T. GIORGI AN. PHILLIPS nduced magnetc feld, µ curl A n R3, and a supercurrent densty j := κ (ψ ψ ψ ψ ) A ψ n. Here κ>0 s the Gnzburg Landau constant determned from the superconductng materal, and the superscrpt denotes complex conjugaton. Accordng to ths theory, the par (ψ, A) s an equlbrum state for the Gbbs free energy ( ) G(ψ, A) := ψ + Aψ κ + ( ψ ) dx (.) + µ curl A he R µ dx + γ ψ ds κ n (see [6], [6]). The constant γ 0 reflects the retardng effect of the materal n the exteror doman on the densty ψ at ; γ s taken to be zero f R n \ s a vacuum and large f the exteror s a magnetc materal. Thus, we consder pars (ψ, A) such that ψ H (; C) H (), A H loc(r n ; R n ), whch are weak solutons to ( ) κ + A ψ ψ + ψ ψ =0 n, ( ) ( ) curl µ curl A + κ (ψ ψ ψ ψ )+A ψ χ =0 nr n, (.) ( ) n κ + A ψ = γψ on, ( ) curl A he L (R n ; R 3 ). µ Here n s the outward normal to at and χ s the characterstc functon for. A prncpal feature of the energy (.) and the solutons to (.) s that they are nvarant under the gauge transformaton where (ψ, A) (ψ, A ), ψ = ψe κη, A = A + η for an arbtrary real valued functon η Hloc (Rn ). Moreover, the ntrnsc quanttes for a soluton are preserved under ths transformaton: ts densty ψ = ψ, magnetc feld µ curl A = µ curl A, current j = j, and the modulus of the dervatve ( κ + A )ψ = ( κ + A)ψ. A soluton s n the normal phase f ψ 0n. Ths s wrtten as (ψ, A) = (0,ha N ), where a N satsfes ( ) curl µ curl a N =0 nr n, (.3) ( ) µ curl a N e L (R n ; R 3 ).

3 ON THE BREAKOWN OF SUPERCONUCTIVITY 39 Such a soluton s called a normal state. It s unquely determned by µ and up to a gauge transformaton; that s, (.3) unquely determnes curl a N. Let κ be fxed. We denote h as the upper crtcal feld for the body, where h := nf{h : normal states are the only solutons to (.) for all h>h }. For the case n whch the body s a bounded doman R 3, we prove the followng statement: Let R 3. Gven κ, µ e, and γ we have h = h(κ, µ e,γ,) < (see Theorem 3.). We show that the normal nducton s contnuous on. In the case that t does not vansh on, we can estmate h. If curl a N 0 n R 3, then there are constants m, φ 0, dependng on µ e and so that ( m ) (.4) h(κ, µ e,γ,) max κ,φκ (see Theorem 3.9). In the classc case where µ e =, t follows that curl a N e, and as such (.4) apples. If µ, then there are constants m and φ such that h max( m κ,φκ) (see Corollary 3.0). We also consder the case of a cylndrcal doman of the form R, where the cross secton s a bounded doman n R wth a C,α boundary and the appled feld H a = he = he 3 s perpendcular to the cross secton. From symmetry the problem reduces to one n two dmensons. We consder ψ(x, y) for (x, y) and A =(A (x, y),a (x, y)) for (x, y) R. The functonal (.) then represents the Gbbs free energy per unt length for the cylnder. We prove the followng theorem. Let R be a cylndrcal body n a parallel appled feld he 3. Gven κ, µ e, and γ there s a fnte upper crtcal feld h, so that f h>h, then the only soluton to (.) wth n =s normal. Moreover, there s a constant φ(µ e, ) so that h(κ, µ e,γ,) max( κ,φκ) (see Theorem.9). Fnally, we consder the case of small κ. We prove the followng result. Let n =wth µ e > 0 or n =3wth µ e =. Then h = O() as κ 0 (see Theorem 4.). It s of nterest to compare these results wth conjectures made by physcsts. For κ fxed, de Gennes and Sant-James have studed the local problem of determnng the smallest value of h for whch all normal states are stable for h h. The nfmum, denoted as h c3, s the value for whch t s possble to have a famly of superconductng solutons bfurcate away from the normal state. In [5] Sant-James and de Gennes dscussed the case of an nfnte slab d <x<d, <y,z< n R 3. The symmetry of the doman reduced the lnear analyss to a one-dmensonal problem. They gave an ansatz for determnng h c3 and predcted lm κ h c3 /κ = c 0 for some constant <c 0 <. Ths can be compared wth our estmates for h from Theorems.9 and 3.9. We have h c3 h = O(κ) asκ. For small κ, physcsts have predcted that h = O() as κ 0 for a slab of fnte thckness d <x<d, <y,z< and that h = O(κ )asκ 0 for the nfntely thck slab <x<0, <y,z< (see [5], [8], and [4]). Our estmate from Theorem 4. gves the result h = O() as κ 0 for our doman. We next comment on past analytc work. In [3] and [4], Bolley and Bolley and Helffer made the ansatz for the slab rgorous and proved asymptotc estmates for h c3.

4 40 T. GIORGI AN. PHILLIPS In [5] they obtaned partal results for estmatng an upper crtcal feld for the slab. They consdered a partcular famly of one-dmensonal functons. For each fxed κ, they showed there s a fnte upper crtcal feld when consderng only solutons n ths famly. In [] Bauman, Phllps, and Tang estmated h c3 for the case of a crcular cylnder, B r R. Ths estmate s relevant here as t plays a central role n our analyss of h for general domans. The central estmates from ths paper were appled by Bauman et al. n [] to nvestgate phase transtons n lqud crystals. A strong analogy exsts between the normal-superconductng phase transton characterzed by (.) and the nematcsmectc phase transton n lqud crystals. The phases for the latter are descrbed n terms of a complex-valued wave functon Ψ and the molecular drector feld n. These are the analogues of the order parameter ψ and the vector potental A, respectvely. Smectc structure s observed through Ψ, and one has Ψ 0 n the nematc phase. The free energy characterzng the nematc-smectc phase transton s the Landau de Gennes energy ntroduced by de Gennes n [9]. The de Gennes model was motvated by the formal analoges between the phase transtons for lqud crystals and superconductvty. There are, however, sgnfcant dfferences between the two theores. For example, the Gnzburg Landau energy (.) s gauge nvarant, whereas the Landau de Gennes energy s not. Furthermore, the drector feld n s requred to satsfy the constrant n =, whereas the correspondng term from superconductvty, A, s not. In secton we consder cylndrcal domans and establsh Theorem.9. In secton 3 we extend these deas to treat bounded domans n R 3. In secton 4 we estmate h for small κ.. Superconductvty wthn an Infntely Long Cylnder n a Parallel Feld. Let (ψ, A) be a weak soluton to (.) wth n = and R. Recall that H a = he 3 s perpendcular to the cross secton. We frst examne the magnetc nducton, curl A, n c. Lemma.. Let (ψ, A) satsfy (.). Then curl A s constant n each component of c. Moreover, curl A = µ e he 3 n the unbounded component. Proof. From (.) we see that curl(curl A) =0 n each component of c. Snce curl(curl A) =( y ( x A y A ), x ( x A y A ), 0), we have that curl A =( x A y A )e 3 s constant n each of these components. The last asserton follows from the fourth equaton n (.). We now determne curl a N. Lemma.. A normal state exsts. Moreover, any normal state (0,ha N ) satsfes curl a N = µe 3. Proof. Consder w =Γ (µ µ e ), where Γ (x) = π ln( x ), x =(x, y). The functon µ µ e has bounded support. As a result, w s well defned wth w Hloc (R ) and w =(µ µ e )nr. Set b N =( w y,w x )+ µe ( y, x). Then curl b N = µe 3 and we see that b N s a weak soluton to (.3), that s, R µ curl b N curl ϕdx = 0 for all ϕ H (R ; R ) such that ϕ has bounded support. Thus, a normal state exsts. Suppose a N s another weak soluton. Takng the dfference of the equatons for b N and a N,weget

5 (.) ON THE BREAKOWN OF SUPERCONUCTIVITY 4 R µ curl(b N a N ) curl ϕdx =0. Let E be the unbounded component of c. From Lemma. we have curl b N = curl a N outsde the bounded set E c. As a result we can take ϕ such that ϕ = b N a N n a neghborhood of E c n (.). Whence, curl b N curl a N. We can now show that a weak soluton has a gauge-equvalent representatve that satsfes a Sobolev estmate. Lemma.3. Let (ζ,b) and (0,ha N ) be weak solutons to (.). Then there s a weak soluton (ψ, A) that s gauge equvalent to (ζ,b) such that (.) A ha N dx C 0 curl(a ha N ) R dx, where C 0 depends only on. Proof. Set curl(b ha N )=fe 3. From Lemmas. and., we have that the support of f s contaned n the bounded set E c and f L (R ). Set w =Γ f; then standard estmates on the Newtonan potental gve w Hloc (R ), w = Γ f, w L (E c ) C 0 () f L (E c ), and w = f (see []). Thus, settng à =( w y,w x ) we have à H loc (R ; R ) and curl à = we 3 = curl(b ha N ). Let A = à + ha N. Then curl(b A) =0. Hence, A = B + η for some η Hloc (R ) and A ha N dx w dx C 0 f dx = C 0 curl(a ha N ) dx. E c E c R We need the followng property for weak solutons. Proposton.4. (see [0]). Let (ψ, A) be a weak soluton to (.); then ψ almost everywhere n. Next we wrte the weak formulaton of (.): (.3) ( ) ( ψ + Aψ ϕ + Aϕ κ κ = γ ψϕ ds κ curl A curl B dx + R µ ) dx + for any ϕ H (), ( ψ )ψϕ dx [ ] κ (ψ ψ ψ ψ )+A ψ B dx =0 for any B H (R ; R ) wth bounded support. Consderng κ ψ + Aψ for (ψ, A) H () Hloc (R ; R ), we have [( ) ] Re ψ + Aψ ψ = κ κ (ψ ψ ψ ψ )+A ψ, [( ) ] Im ψ + Aψ ψ κ = κ ψ. Thus, (.4) {[ ] [ ]} ψ ψ + Aψ = κ κ (ψ ψ ψ ψ )+A ψ ψ + κ ψ ψ

6 4 T. GIORGI AN. PHILLIPS for almost every x such that ψ 0. Moreover, snce ψ = 0 almost everywhere on the set {ψ =0}, t s consstent to defne the term n braces equal to zero on ths set. We conclude that (.4) holds almost everywhere. Lemma.5. Let (ψ, A) and (0,ha N ) be weak solutons satsfyng (.). Then there s a constant C = C (,µ e ) such that (.5) ( + κha N )ψ dx C κ ψ dx. Proof. Let ϕ = ψ n the frst equaton of (.3). Usng (.4), Proposton.4, and γ 0, we obtan ( [ ] κ ψ ) + κ (ψ ψ ψ ψ )+A ψ ψ dx (.6) ( ) = κ + A ψ dx ( ψ ) ψ dx ψ dx. Ths nequalty s also vald for n = 3. Consder the second equaton n (.) for the solutons (ψ, A) and (0,ha N ). Takng the dfference of ther respectve weak equatons we have [ ] R µ curl(a ha N) curl B dx = κ (ψ ψ ψ ψ )+A ψ ψ ψ Bdx. Usng (.6) and the Cauchy Schwarz nequalty, we see R µ curl(a ha N) curl B dx ε ψ dx + ε ψ B dx for any ε>0. Let B be such that B = A ha N n E c, where E s the unbounded component of c. Then snce curl(a ha N )=0 n E and ψ we derve R µ curl(a ha N) dx ε ψ dx + ε A ha N dx. Combnng ths nequalty wth (.), we see that we can take ε suffcently small so that (.7) A ha N dx M ψ dx for some constant M = M(dam,µ e ). Next we wrte ( ) ( ) (.8) κ + A ψ = κ + ha N ψ +(A ha N )ψ. We wll use the elementary nequalty (.9) c b c + b for c, b C.

7 ON THE BREAKOWN OF SUPERCONUCTIVITY 43 Let ( κ + A)ψ = b and ( κ + ha N)ψ = c. Then usng (.6) and (.8) we derve ( ) κ + ha N ψ dx ψ dx + A ha N ψ dx. Snce ψ, we can apply (.7) to obtan ( + κha N )ψ dx ( + M)κ ψ dx. We set C = ( + M) and the lemma s proved. We see that f a superconductng state (.e., a soluton wth ψ 0) exsts, then (.5) mples that the prncpal egenvalue for ( +κha N ) on s bounded by C κ. We wll show that there exsts a constant φ such that f h>max( κ,φκ), then the prncpal egenvalue s greater than C κ. It then follows for such κ and h that there are only normal solutons to (.). The correspondng egenfunctons are expected to take the form of a boundary layer. The followng lemma gves a way of measurng to what extent functons can concentrate near the boundary. For a set O we defne the τ-neghborhood n O of O by O τ = {x O : dst(x, O) <τ}. Lemma.6. Let O be a bounded doman n R n wth a C boundary. Gven λ 0 > 0, there s a constant d(λ 0, O) > 0 such that whenever (.0) f dx λ f dx for some f H (O) wth λ λ 0, then (.) f dx O O O O\O dλ f dx. Proof. Let N k=0f k be an open cover for O such that F 0 O and such that for each k, k N, wehave F k O = {(x,x n ):g k (x ) <x n <g k (x )+δ, x <δ }, where δ and δ are postve constants, (x,x n ) are sutably rotated and translated coordnates, and g k (x )=x n characterzes O F k. We can further assume wthout loss of generalty that g k ( ) s defned for x δ, g k <, and (.) {(x,x n ):g k (x ) <x n <g k (x )+4δ, x < δ } O, {(x,x n ):g k (x ) 4δ <x n <g k (x ), x < δ } R n \ O. Let f H (O), 0 t, v δ, and fx k. We have f ds f ds {x n g k (x )=v} F k {x n g k (x )=t} F k f (x,g k (x )+v) f (x,g k (x )+t) + g k dx { x δ } O F k (f ) x n dx.

8 44 T. GIORGI AN. PHILLIPS Integratng n t from 0 to δ and then dvdng by δ gves f ds ( ) f dx +δ f f dx. {x n g k (x )=v} F k δ O O Next ntegrate v from 0 to d λ, where 0 <d<λ 0δ / s to be determned. We derve (.3) {0<x n g k (x ) d λ } F k f dx d ( ) f dx +δ f f dx. λδ O O From the assumptons on g k, for k N we have { (.4) F k O d F k 0 <x n g k (x ) d } λ λ f d λ s small enough. Indeed, f ths s false for some k we can fnd x =(x,x n ) such that x δ, x n >g k (x )+ d λ, and y =(y,y n ) O such that x y < d λ. If d λ < mn(δ,δ ), we have x y < d λ <δ, mplyng that y < δ and y n x n < d λ <δ. We clam that y n = g k (y ). In fact, y n g k (y ) y n x n + x n g k (x ) + g k (x ) g k (y ), and each term on the rght s bounded by δ. We have shown ths for the frst one. Ths s true for the second snce x F k O. For the last term we use g k < and x y < d λ <δ. Thus, y n g k (y ) < 3δ. From (.) we see that the only possblty for such a y O s y n = g k (y ). As a result, g k (x ) g k (y ) g k (x ) x n g k (y ) x n > d λ d λ = d λ. On the other hand, snce g k < wehave g k (x ) g k (y ) < x y < d λ, and ths s a contradcton. Usng (.3) and (.4) and summng on k for k N, we obtan ( ) f d dx M f dx + λ f dx + λ f dx, λ O O O O dλ where M = M (δ,n). Usng (.0) we have O dλ f dx M d f dx, O where M = M (δ,n,λ 0 ). Settng d = mn( M, λ0δ f dx f dx. O dλ O, λ0δ The asserton (.) follows from ths nequalty. We wll use the followng result from [] for B r (0) R. ), we conclude that

9 ON THE BREAKOWN OF SUPERCONUCTIVITY 45 Proposton.7. There s a contnuous functon σ( ) :t [0, ) R wth σ(t) > 0 for t>0 for whch lm σ(t) exsts wth 0 < lm σ(t) <, and such that t t (.5) B r(0) ) ( + ω ( y, x) ζ for all ζ H (B r (0)) and ω 0. Indeed, n [, secton ], t s shown that nf ζ L = B r ζ W, (B r;c) ( + ω dx ω σ(ωr) ζ dx B r(0) ) ( y, x) ζ dx ω σ, where σ = σ(ωr). Furthermore, σ s characterzed by σ(t) = nf n z σ(t, n), where for each n, σ(t, n) s analytc and postve on 0 <t<. Moreover, lm t 0 σ(t, 0)=0. In [, secton 6], t s also shown that σ(t) = mn σ(t, n) for 0 t n 0. 0 n n 0 As a result, t follows that σ(t) s postve and contnuous. The lm t σ(t) s analyzed n [, secton 6], as well. Remark. Ifb s another vector feld such that b H (B r (0); R ) wth curl b = e 3, then (.5) s also vald wth ( y, x) replaced by b. Indeed, we can defne a functon q H (B r (0)) such that q = b ( y, x). Wth ths we can defne a local gauge transformaton ζ = ζe ωq, b = ( y, x)+ q, for whch ζ H (B r (0)) provded ζ H (B r (0)). Moreover, ) ( + ω b)ζ = ( + ω ( y, x) ζ and ζ = ζ so that (.6) B r(0) ( + ω b)ζ dx ω σ(ωr) ζ dx B r(0) for all ζ H (B r (0)), b H (B r (0); R ) such that curl b = e 3. We now derve an estmate smlar to (.6) for provded ω s bounded away from zero. Lemma.8. Gven m>0, there s a constant C = C (m, ), 0 <C, such that f ω m, then (.7) C ω ζ dx ( + ω b)ζ dx for all ζ H () and b H (; R ) for whch curl b = e 3. Proof. Let ζ H () such that ζ dx > 0 and ( + ω b)ζ dx ω ζ dx

10 46 T. GIORGI AN. PHILLIPS for some ω, ω m. If no such ζ exsts, then (.7) s vald wth C = and we are fnshed. From (.4) we see ζ ( + ω b)ζ. Thus, ζ dx ω ζ dx. As a result, we can apply Lemma.6 to conclude that (.8) ζ dx ζ dx. Next we choose a cover for \ d ω \ dω consstng of a fnte collecton of dsks {B d (x k ), ω k =,...,N(ω)}, each contaned n n such a way that N(ω) k= χ B dω (x k ) K, where K s ndependent of ω. We see ( + ω b)ζ dx K N k= B d ω (x k ) N k= Usng Proposton.7, the last term bounds K ω σ(d) N k= B dω (x k ) ζ dx K (d)ω ( + ω b)ζ dx B dω (x k ) ( + ω b)ζ dx. \ dω ζ dx K ω ζ dx, where the fnal nequalty follows from (.8). Set C = K /. Ths chan of nequaltes establshes the lemma. We now establsh the prncpal result n ths secton. Here we prove the exstence of an upper crtcal feld h and obtan a bound for t as κ and κ 0. Theorem.9. There s a constant φ = φ(µ e, ) so that f h>max( κ,φκ), then any weak soluton for (.) wth n =s normal. Proof. Let (0, a N ) be a normal state for (.3) and (ψ, A) be a weak soluton for (.). A state s normal f and only f ts entre gauge equvalence class s normal. Therefore, we can assume wthout loss of generalty that (ψ, A) and (0,ha N ) satsfy (.). Set ω = hκ. Then ω by hypothess. We apply (.7) wth m =to derve (.9) C hκ ψ dx ( + hκa N )ψ dx, and by Lemma.5, the rght-hand sde of (.9) s bounded by C κ ψ dx. Let φ = C /C. We have h ψ dx φκ ψ dx. By assumpton, h>φκ. Hence t must hold that ψ dx =0.

11 ON THE BREAKOWN OF SUPERCONUCTIVITY Three-mensonal Bodes. In ths secton, we consder a superconductng body gven by a bounded doman R 3 subjected to a unform appled feld H a = he. We wll assume wthout loss of generalty that e = e 3 throughout ths secton. enote by Ȟ (R 3 ) the completon of C0 (R 3 ; R 3 ) wth respect to the norm ) B Ȟ (R 3 ) (R = B dx. 3 One can show that elements B Ȟ (R 3 ) satsfy the followng relatonshps: (3.) B L 6 (R 3 ;R 3 ) θ B Ȟ (R 3 ), where θ s ndependent of B, and (3.) B Ȟ (R 3 ) = R 3 ( dv B + curl B )dx (see []). In order to represent magnetc felds we need the followng lemma. Lemma 3.. Let g L (R 3 ; R 3 ) such that dv g =0n (R 3 ). Then there s a unque u Ȟ (R 3 ) such that curl u = g and dv u =0. Proof. Consder k Γ g, where Γ(x) =Γ 3 (x) = 3ω 3 x s the Newtonan potental for R 3 and k 3. We clam that k Γ g Ȟ (R 3 ). To see ths, let us frst assume that g has bounded support. Then k Γ g exsts as a weakly sngular ntegral and and k Γ g = O( x ) ( k Γ g) = O( x 3 )as x. Let ϕ R (x) be a standard C cutoff functon such that ϕ R = for x R and ϕ R = 0 for x R +. It follows drectly that {ϕ R(n) ( k Γ g)} s a Cauchy sequence n Ȟ (R 3 ) that converges to k Γ g pontwse for any sequence R(n). Thus, k Γ g Ȟ (R 3 ), assumng g has bounded support. Fnally, by standard L -sngular ntegral theory, (3.3) k Γ g Ȟ (R 3 ) g L (R 3 ;R 3 ), and as a consequence, k Γ g Ȟ (R 3 ) for all g L (R 3 ; R 3 ). efne u : g L (R 3 ; R 3 ) Ȟ (R 3 )by u(g) = ( Γ g 3 3 Γ g, 3 Γ g Γ g 3, Γ g Γ g ). Let g ε be a mollfcaton of g; then g ε g n L as ε 0, and dv g ε = 0 for each ε>0. We defne g ε,r = ϕ R g ε. Note that dv(g ε,r )= ϕ R g ε. If we choose sequences ε(n) 0 and R(n) as n, then g n = ϕ R(n) g ε(n) C 0 (R 3 ; R 3 ), g n g n L, and dv g n 0nL as n. Set w n =Γ g n. These are well defned snce the g n have bounded support. Usng (3.3), we see that curl w n = u(g n ) u(g) nȟ as n.

12 48 T. GIORGI AN. PHILLIPS Consder (3.4) curl u(g n )= curl curl w n = w n (dv w n )=g n (dv w n ). We know that (dv w n )= Γ (dv g n ) 0 n Ȟ (R 3 )asn snce dv g n 0n L. Thus, usng (3.) we conclude that (dv w n ) 0 n L 6 as n. Furthermore, we have curl u(g n ) curl u(g) and g n g n L as n. As a consequence, curl u(g) =g n R 3. Snce u(g n )= curl w n we have dv u(g n ) = 0, whch mples dv u(g) =0. Fnally, usng (3.) we see that u s unque n Ȟ (R 3 ). We can apply the precedng lemma to characterze weak solutons. Lemma 3.. Let (ζ,b) be a weak soluton to (.). Then there s a gaugeequvalent soluton (ψ, A) such that dv A =0and (A µeh ( y, x, 0)) Ȟ (R 3 ). Moreover, f ( ψ, Ã) s another such soluton, then ψ = aψ for some a C, a =, and à = A. Proof. Set g = (curl B µ e he 3 ) L (R 3 ; R 3 ). From the prevous lemma there s a unque element u Ȟ (R 3 ) such that curl u = g and dv u = 0. Therefore, we fnd A = u + µeh ( y, x, 0). We now characterze the normal state n three dmensons. Lemma 3.3. There s a unque normal state satsfyng (.3) such that (a N µ e ( y, x, 0)) Ȟ (R 3 ) and dv a N =0. Proof. Consder the strctly convex functonal E(b) =G(0, b)+ (dv b) dx = R 3 dx + R 3 ( µ µ curl b e 3 + (dv b) ) dx for the class S = {b :(b µe ( y, x, 0)) Ȟ (R 3 )}. A unque equlbrum exsts that also mnmzes E( ). Let b be ths equlbrum. If dv b 0, then by Lemma 3. we can fnd another vector feld b S such that curl b = curl b and dv b = 0. Ths would mply E( b) <E( b), whch s mpossble. Thus, dv b = 0. It follows that b satsfes (.3), and as a result a normal state (0, a N ) exsts. Conversely, a normal state satsfyng the hypothess s an equlbrum for E( ), and so a N s unque. Recall that the nducton curl a N, not a N, s the physcally relevant quantty. Below we show that t s unquely determned. Lemma 3.4. Let (0, a N ) be a normal state. Then curl a N s unquely determned, curl a N s harmonc n R 3 \, and Moreover, f dv a N = 0, then curl a N C,α () C,α ( c ). a N C,α () C,α ( c ). Proof. Usng Lemma 3. we see that any normal state s gauge equvalent to the normal state descrbed n Lemma 3.3. Snce a gauge transformaton leaves the curl of a vector feld nvarant, we conclude that curl a N s unquely determned for solutons to (.3). We can use the frst equaton n (.3) to prove that there exsts a functon p H loc (R3 ) such that curl a N = µ p. Snce (3.5) dv(µ p) =0nR 3

13 ON THE BREAKOWN OF SUPERCONUCTIVITY 49 and µ s constant on the components of R 3 \, the functon p (and thus curl a N ) s harmonc n each component. We apply the results from [3, Chap. 5, secton 4] to the soluton p for (3.5), to derve that p C,α () C,α ( c ), and as a consequence, curl a N C,α () C,α ( c ). Assume dv a N = 0. Let U be an open neghborhood of and consder w H (U; R 3 ) such that w = curl a N n U. From [3, Chap. 5], we have w C,α () C,α (U\). The dentty curl(curl w)+ w = (dv w) nu yelds from whch we obtan curl(curl w + a N )= (dv w), curl curl(curl w + a N )=0 n (U). By hypothess, dv(curl w + a N ) = 0. Whence, from the dentty above, (curl w + a N ) = curl curl(curl w + a N )=0 n (U). Ths mples (curl w + a N ) C (U), and we conclude that a N C,α () C,α ( c ). Consder the case µ. Gven e, we can fnd a lnear functon a(x) such that curl a e. Clearly a satsfes (.3). It follows from the prevous lemma, then, that curl a N e when µ. We now derve a Sobolev estmate analogous to Lemma.3. Lemma 3.5. Let (ζ,b) be a weak soluton to (.). Let (ψ, A) be the gaugeequvalent soluton found n Lemma 3. and (0,ha N ) be the normal state found n Lemma 3.3. Then there s a constant C 0 dependng only on such that A ha N dx C 0 curl(a ha N ) R dx. 3 Proof. Usng (3.) and (3.) we see A ha N L 6 (R 3 ;R 3 ) θ (A ha N ) L (R 3 ;R 3 ) = θ curl(a ha N ) L (R 3 ;R 3 ). Snce s bounded we have A ha N L (;R 3 ) M() A ha N L 6 (;R 3 ) and the lemma follows. We proceed to derve the three-dmensonal counterpart to Lemma.5. Lemma 3.6. Let (ψ, A) and (0,ha N ) be as n Lemma 3.5. Then there s a constant C = C (,µ e ) so that ( + hκa N )ψ dx C κ ψ dx. Proof. We proceed just as n Lemma.5 to obtan R µ [curl(a ha N) curl B]dx ε ψ dx + ε 3 ψ B dx

14 50 T. GIORGI AN. PHILLIPS for any ε>0 and B H (R 3 ; R 3 ) wth bounded support. However, snce A ha N Ȟ (R 3 ), we can take B = B j A ha N n Ȟ (R 3 )asj. As a result, we have R µ curl(a ha N) dx ε ψ dx + ε ψ A ha N dx. 3 The remander of the proof s just as before. We next gve a three-dmensonal analogue for the egenvalue estmate from []. Let v R 3 \{0} such that v = and x 0 R 3. Let T (x 0,r,v) be a cylnder wth central axs parallel to v, heght r, and whose mddle cross secton s the dsk of radus r wth center x 0. Lemma 3.7. Let b H (T (x 0,r,v); R 3 ) such that curl b = v. Then (3.6) ( + ω b)ζ dx ω σ(ωr) ζ dx T for all ζ H (T ), where σ( ) s as n Proposton.7. Proof. We frst transfer the problem to T (0,r,e 3 )=B r (0) ( r, r). Let Q SO(3) such that e 3 = Qv, and set y(x) =Q(x x 0 ). Then y : x T (x 0,r,v) T (0,r,e 3 ). Gven ζ H (T (x 0,r,v)) we defne ξ(y) =ζ(x(y)). Then ( y + ω bq t )ξ(y) =( x + ω b)ζ(x)q t. By changng varables we see that (3.6) s equvalent to showng the followng nequalty: ( + ω bq t )ξ dy ω σ(ωr) ξ dy. T (0,r,e 3) T (0,r,e 3) For any w R 3,wehave w t Q t w t Q t w t Qw curl y (bq t ) = det y = det x Q t = det x = w curl x b. bq t bq t b Therefore, curl y (bq t )=Q(curl x b)=qv = e 3. By changng the gauge f necessary, we can assume bq t = ( y, x, 0). Then ( + ω bq t )ξ dy (( x, y, 0) + ω bq t )ξ dy T (0,r,e 3) T (0,r,e 3) r ) = (( x, y, 0) + ω ( y, x, 0) ξ(x, y, z) r B r dx dy dz r ω σ(ωr) ξ(x, y, z) dx dy dz r B r = ω σ(ωr) ξ dy, T (0,r,e 3) where we have appled Proposton.7 for each r z r. T

15 ON THE BREAKOWN OF SUPERCONUCTIVITY 5 We go on to prove the three-dmensonal counterpart to the egenvalue estmate n Lemma.8. Lemma 3.8. Let (0, a N ) be the normal state from Lemma 3.3. Assume curl a N 0 n. Then there exst constants m and 0 <C so that f ω m, t holds that (3.7) C ω ζ dx ( + ω a N )ζ dx for all ζ H (). Proof. We argue as n Lemma.8. There exsts a constant d>0 so that, gven ξ H () and ω, ether (3.7) s true wth C = and ξ = ζ or (3.8) ξ dx ξ dx. \ dω Assume the latter. In ths case we cover \d by a famly of cylnders {T k : k = ω,...,n(ω)} such that T k = T (x k, d ω, curl a N(x k )/ curl a N (x k ) ) wth T k for each k and N k= χ T k K, where K s ndependent of ω for ω. As a consequence, (3.9) ( + ω a N )ξ dx K N k= T k ( + ω a N )ξ dx. In each T k we wrte a N (x) =l k (x)+q k (x), where l k (x) =a N (x k )+ a N (x k ) (x x k ). Note that curl l k (x) = curl a N (x k ). Usng (.9), for each k we obtan ( + ω a N )ξ dx ( + ω l k )ξ dx ω 4 q k ξ dx. T k T k T k From Lemma 3.7 we have T k ( + ω l k )ξ dx ω curl a N (x k ) σ ω M 0 T k ξ dx, ( ) curl a N (x k ) d ξ dx T k where M 0 > 0 depends on nf curl a N > 0 and the structure of σ( ) (see Proposton.7). Snce a N C,α () wehave q k M (dam T k ) +α M ω α, where M s ndependent of k. As a result, we see for each k that (3.0) ( ) ( + ω a N )ξ M0 dx T k ω M ω α ξ dx M 0 T k 4 ω ξ dx, T k provded ω m = m(,µ e ) suffcently large.

16 5 T. GIORGI AN. PHILLIPS From (3.9) and (3.0), then, ( + ω a N )ξ dx M 3 ω N k= T k ξ dx for some M 3 > 0 ndependent of the cover. Usng \d N ω k= T k and (3.8), we derve ξ dx ξ dx ξ dx. \ dω N k= T k Settng C = M3, we have our lemma. The followng theorem s proved n the same manner as Theorem.9. We establsh the exstence of h and derve an upper bound for t provded curl a N does not vansh on. Theorem 3.9. Assume that curl a N 0 n. There are constants m and φ, dependng on and µ e, so that f h>max( m κ,φκ), then any weak soluton to (.) wth n =3s normal. For the case µ, we have curl a N, and we can recover the followng result. Corollary 3.0. If µ e =, then there exst constants m and φ dependng on so that f h>max( m κ,φκ), then any weak soluton to (.) wth n =3s normal. In general, one does not know whether or not curl a N vanshes somewhere n. Nevertheless, snce curl a N s harmonc n t can vansh only on a small set. Lemma 3.. Let (0, a N ) be a normal state. Then L 3 ({x : curl a N (x) = 0}) =0. Proof. Snce curl a N s harmonc n, ether L 3 ({x : curl a N (x) =0}) =0 or curl a N 0 n. From Lemma 3.4, we know that there s a functon p C,α () C,α ( c ) C(R 3 ) such that Hence, p satsfes curl a N = µ p n R 3. (3.) { dv(µ p) =0 nr 3, ( p e 3 ) L (R 3 ). Assume that curl a N 0 n. Then p = constant = p 0 n. It follows from the frst equaton n (3.) that p solves p =0 n c, p = p 0, p n =0 on c, where n s the exteror normal to. The Cauchy problem has the unque soluton p = p 0. Ths would contradct the second equaton n (3.). To conclude, we show that there s a fnte upper crtcal feld for each κ. Theorem 3.. Let κ, µ e, and γ be fxed. There s a constant h = h(κ, µ e,γ,) so that f h>h, then any weak soluton to (.) wth n =3s normal. Proof. Let (0, a N ) be as n Lemma 3.3. Assume that there exsts a sequence {(ψ j, A j )}, where for each j the par s as n Lemma 3.5, solvng (.) wth h = h j for whch lm j h j = and ψ j dx > 0.

17 ON THE BREAKOWN OF SUPERCONUCTIVITY 53 Set ϕ j (x) = ψ j (x) / ψ j L (). From (.6) we have ϕ j dx κ, and we can fnd a subsequence ϕ j ϕ 0 n L () asj wth ϕ 0 L () =. From Lemma 3. the set Q = {x : curl a N (x) =0} s a closed set of measure zero. It follows that there exsts a ball B r \Q such that B r ϕ 0 dx δ >0 and nf Br curl a N > 0. Note that for j suffcently large we have ψ j dx δ ψ j dx. B r For each j we cover B r by a fnte famly of cylnders, {T (x k,ω j, curl a N (x k )/ curl a N (x k ) )} = {T kj, k N(j)}, such that ωj = h jκ, x k B r, and such that the famly has overlap of at most K unformly n x ndependent of j. Forj suffcently large, each T kj B r and, as n Lemma 3.8, we derve ( + ωj a N )ψ j dx ωj M 0 ψ j dx, T kj T kj where M 0 > 0 depends on nf Br curl a N and a N C,α () but not on j. Whence we can wrte N(j) ( + h j κa N )ψ j dx K ( + h j κa N )ψ j dx T kj k= M h j κ ψ j dx M h j κ N(j) k= T kj ψ j dx M h j κ B r ψ j dx, where M and M are postve and depend on nf Br curl a N, a N, and δ. Lemma 3.6 we know ( + h j κa N )ψ j dx C κ ψ j dx, From whch leads to M h j κ ψ j dx C κ ψ j dx for all j suffcently large. Snce h j as j, we must have ψ j dx = 0 for j large, and ths s a contradcton. 4. Estmates for Small κ. In ths secton we consder h for small κ n cases where curl a N e n. Theorem 4.. Let n =wth µ e > 0 or n =3wth µ e =. Then h = O() as κ 0. Proof. We consder the case n = 3. The argument for n = s dentcal.

18 54 T. GIORGI AN. PHILLIPS Let κ and assume that (ψ, A) solves (.) wth ψ 0. From (.6) we have ψ dx κ ψ dx ψ dx. If we apply Lemma.6 wth λ = λ 0 =, we fnd that there s a constant d>0so that ψ dx ψ dx. \ d We take r, 0<r<d, to be determned and cover \ d by a famly of cylnders {T (x k,r,e)} such that the cylnders have fnte overlap ndependent of r and each s contaned n. Then, argung just as n Lemma.8, there exsts a constant C > 0 for whch C hκ σ((hκ) r) ψ dx ( + hκa N )ψ dx. Applyng Lemma 3.6 to the rght-hand sde and recallng that ψ dx > 0, the followng nequalty holds: (4.) C hκ σ((hκ) r) C κ. From Corollary 3.0 we have hκ M for some M < for all κ. We now choose r so that M r. The reason for ths choce s that f 0 <t, then σ(t) s the prncpal egenvalue for the Sturm Louvlle problem g (s)+ g (s) s g(s) = σ(t)g(s) for 0 <s<t, s 4 g (t) =0andg s bounded (see [, secton 6]). As such, from [7, Prop. 3.4], we have σ(t) = t 8 + o(t )ast 0. It follows that there s a constant M > 0 so that M t σ(t) for 0 t. Combnng ths estmate wth (4.), we derve h C C M r for 0 <κ. Snce ths holds for all superconductng solutons, we conclude that h C C M r for 0 <κ. The length scale for varatons n superconductng solutons s κ. Because of ths, t s of nterest to consder domans wth dmensons comparable to κ. To ths end, let R n and defne the dlated doman (κ) ={x R n : κx }.

19 ON THE BREAKOWN OF SUPERCONUCTIVITY 55 Then dam((κ)) = κ dam(). Let µ so that curl a N e n (κ). Consder (ψ, A) satsfyng (.) on (κ). Let and ψ(x) =ψ(κ x) for x Ã(x) =A(κ x) for x R 3. Then ( ψ, Ã) satsfes ( + Ã) ψ ψ + ψ ψ =0 n, ( curl à = κ (4.) ( ψ ψ ψ ψ ) )+à ψ χ n R 3, n ( + Ã) ψ = γ ψ on, (curl à hκ e) L (R 3 ; R 3 ). Assume that 0 <κ. Argung as n Lemma.5, takng nto account the multple κ of the rght-hand sde of the second equaton n (4.), we obtan the analogue of (.7), à hκ a N dx κ 4 M ψ dx, where M = M(dam ). Usng κ we then obtan the analogue of (.5), ( + hκ a N ) ψ dx C κ 4 ψ dx. It follows as n Theorem.9, then, that there s a constant φ so that f hκ φκ 4 then ψ 0. Returnng to (ψ, A), we conclude that h(κ,,γ,(κ)) = O(κ 3 ) as κ 0. As we noted n the ntroducton, there are superconductng solutons for the case of the slab <x<0, <y,z< for whch h dverges as κ 0 (see [5]). However, ths occurs at the slower rate h = O(κ ). REFERENCES [] P. Bauman, M.C. Calderer, C. Lu, and. Phllps, The Phase Transton between Chral Nematc and Smectc A Lqud Crystals, preprnt. [] P. Bauman,. Phllps, and Q. Tang, Stable nucleaton for the Gnzburg Landau system wth an appled magnetc feld, Arch.Ratonal Mech.Anal., 4 (998), pp. 43. [3] C. Bolley, Modélzaton de champ de retard à la condensaton d un supraconducteur par un problème de bfurcaton, Math.Modél.Anal.Numér., 6 (99), pp [4] C. Bolley and B. Helffer, An applcaton of sem-classcal analyss to the asymptotc study of the supercoolngfeld of a superconductngmateral, Ann.Inst.H.Poncaré, 58 (993), pp [5] C. Bolley and B. Helffer, On the asymptotcs of the crtcal felds for the Gnzburg Landau equaton, n Progress n Partal fferental Equatons: The Metz Surveys, Ptman Res. Notes Math.Ser.34, Longman Scentfc and Techncal, Harlow, UK, 994, pp.8 3.

20 56 T. GIORGI AN. PHILLIPS [6] J. Chapman, Q. u, and M. Gunzburger, A Gnzburg Landau type model for superconductng/normal junctons ncludng Josephson junctons, European J.Appl.Math., 6 (995), pp [7] M. auge and B. Helffer, Egenvalues varaton, I, Neumann problem for Sturm Louvlle operators, J.fferental Equatons, 04 (993), pp [8] P. G. de Gennes, Superconductvty: Selected topcs n sold state physcs and theoretcal physcs, n Proceedngs of the 8th Latn Amercan School of Physcs, Caracas, 966. [9] P. G. de Gennes, An analogy between superconductvty and semectcs A, Sold State Commun., 0 (97), pp [0] Q. u, M.. Gunzburger, and J. S. Peterson, Analyss and approxmaton of the Gnzburg Landau model of superconductvty, SIAM Rev., 34 (99), pp []. Glbarg and N. S. Trudnger, Ellptc Partal fferental Equatons of Second Order, Sprnger-Verlag, Berln, Hedelberg, 977. [] O. A. Ladyzhenskaya, The Mathematcal Theory of Vscous Incompressble Flow, Gordon and Breach, New York, 969. [3] O. A. Ladyzhenskaya, The Boundary Value Problems of Mathematcal Physcs, Sprnger- Verlag, Berln, New York, 985. [4] The Orsay Group, Strongfeld effects at the surface of a superconductor, n Quantum Fluds,.F.Bewer, ed., North-Holland, Amsterdam, 966, p.6. [5]. Sant-James and P. G. de Gennes, Onset of superconductvty n decreasngfelds, Phys. Lett., 7 (963), pp [6] M. Tnkham, Introducton to Superconductvty, nd ed., McGraw-Hll, New York, 996.

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity

Uniqueness of Weak Solutions to the 3D Ginzburg- Landau Model for Superconductivity Int. Journal of Math. Analyss, Vol. 6, 212, no. 22, 195-114 Unqueness of Weak Solutons to the 3D Gnzburg- Landau Model for Superconductvty Jshan Fan Department of Appled Mathematcs Nanjng Forestry Unversty

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Math 217 Fall 2013 Homework 2 Solutions

Math 217 Fall 2013 Homework 2 Solutions Math 17 Fall 013 Homework Solutons Due Thursday Sept. 6, 013 5pm Ths homework conssts of 6 problems of 5 ponts each. The total s 30. You need to fully justfy your answer prove that your functon ndeed has

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation

Asymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation Nonl. Analyss and Dfferental Equatons, ol., 4, no., 5 - HIKARI Ltd, www.m-har.com http://dx.do.org/.988/nade.4.456 Asymptotcs of the Soluton of a Boundary alue Problem for One-Characterstc Dfferental Equaton

More information

REAL ANALYSIS I HOMEWORK 1

REAL ANALYSIS I HOMEWORK 1 REAL ANALYSIS I HOMEWORK CİHAN BAHRAN The questons are from Tao s text. Exercse 0.0.. If (x α ) α A s a collecton of numbers x α [0, + ] such that x α

More information

Affine transformations and convexity

Affine transformations and convexity Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens THE CHINESE REMAINDER THEOREM KEITH CONRAD We should thank the Chnese for ther wonderful remander theorem. Glenn Stevens 1. Introducton The Chnese remander theorem says we can unquely solve any par of

More information

Exercise Solutions to Real Analysis

Exercise Solutions to Real Analysis xercse Solutons to Real Analyss Note: References refer to H. L. Royden, Real Analyss xersze 1. Gven any set A any ɛ > 0, there s an open set O such that A O m O m A + ɛ. Soluton 1. If m A =, then there

More information

Canonical transformations

Canonical transformations Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Perron Vectors of an Irreducible Nonnegative Interval Matrix Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

Group Analysis of Ordinary Differential Equations of the Order n>2

Group Analysis of Ordinary Differential Equations of the Order n>2 Symmetry n Nonlnear Mathematcal Physcs 997, V., 64 7. Group Analyss of Ordnary Dfferental Equatons of the Order n> L.M. BERKOVICH and S.Y. POPOV Samara State Unversty, 4430, Samara, Russa E-mal: berk@nfo.ssu.samara.ru

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

3 Basic boundary value problems for analytic function in the upper half plane

3 Basic boundary value problems for analytic function in the upper half plane 3 Basc boundary value problems for analytc functon n the upper half plane 3. Posson representaton formulas for the half plane Let f be an analytc functon of z throughout the half plane Imz > 0, contnuous

More information

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS These are nformal notes whch cover some of the materal whch s not n the course book. The man purpose s to gve a number of nontrval examples

More information

First day August 1, Problems and Solutions

First day August 1, Problems and Solutions FOURTH INTERNATIONAL COMPETITION FOR UNIVERSITY STUDENTS IN MATHEMATICS July 30 August 4, 997, Plovdv, BULGARIA Frst day August, 997 Problems and Solutons Problem. Let {ε n } n= be a sequence of postve

More information

STEINHAUS PROPERTY IN BANACH LATTICES

STEINHAUS PROPERTY IN BANACH LATTICES DEPARTMENT OF MATHEMATICS TECHNICAL REPORT STEINHAUS PROPERTY IN BANACH LATTICES DAMIAN KUBIAK AND DAVID TIDWELL SPRING 2015 No. 2015-1 TENNESSEE TECHNOLOGICAL UNIVERSITY Cookevlle, TN 38505 STEINHAUS

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

The Feynman path integral

The Feynman path integral The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space

More information

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP C O L L O Q U I U M M A T H E M A T I C U M VOL. 80 1999 NO. 1 FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP BY FLORIAN K A I N R A T H (GRAZ) Abstract. Let H be a Krull monod wth nfnte class

More information

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY. Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Appendix B. Criterion of Riemann-Stieltjes Integrability

Appendix B. Criterion of Riemann-Stieltjes Integrability Appendx B. Crteron of Remann-Steltes Integrablty Ths note s complementary to [R, Ch. 6] and [T, Sec. 3.5]. The man result of ths note s Theorem B.3, whch provdes the necessary and suffcent condtons for

More information

Perfect Competition and the Nash Bargaining Solution

Perfect Competition and the Nash Bargaining Solution Perfect Competton and the Nash Barganng Soluton Renhard John Department of Economcs Unversty of Bonn Adenauerallee 24-42 53113 Bonn, Germany emal: rohn@un-bonn.de May 2005 Abstract For a lnear exchange

More information

DIFFERENTIAL FORMS BRIAN OSSERMAN

DIFFERENTIAL FORMS BRIAN OSSERMAN DIFFERENTIAL FORMS BRIAN OSSERMAN Dfferentals are an mportant topc n algebrac geometry, allowng the use of some classcal geometrc arguments n the context of varetes over any feld. We wll use them to defne

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Maximizing the number of nonnegative subsets

Maximizing the number of nonnegative subsets Maxmzng the number of nonnegatve subsets Noga Alon Hao Huang December 1, 213 Abstract Gven a set of n real numbers, f the sum of elements of every subset of sze larger than k s negatve, what s the maxmum

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014)

10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014) 0-80: Advanced Optmzaton and Randomzed Methods Lecture : Convex functons (Jan 5, 04) Lecturer: Suvrt Sra Addr: Carnege Mellon Unversty, Sprng 04 Scrbes: Avnava Dubey, Ahmed Hefny Dsclamer: These notes

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS BOUNDEDNESS OF THE IESZ TANSFOM WITH MATIX A WEIGHTS Introducton Let L = L ( n, be the functon space wth norm (ˆ f L = f(x C dx d < For a d d matrx valued functon W : wth W (x postve sem-defnte for all

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

Geometry of Müntz Spaces

Geometry of Müntz Spaces WDS'12 Proceedngs of Contrbuted Papers, Part I, 31 35, 212. ISBN 978-8-7378-224-5 MATFYZPRESS Geometry of Müntz Spaces P. Petráček Charles Unversty, Faculty of Mathematcs and Physcs, Prague, Czech Republc.

More information

Randić Energy and Randić Estrada Index of a Graph

Randić Energy and Randić Estrada Index of a Graph EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 5, No., 202, 88-96 ISSN 307-5543 www.ejpam.com SPECIAL ISSUE FOR THE INTERNATIONAL CONFERENCE ON APPLIED ANALYSIS AND ALGEBRA 29 JUNE -02JULY 20, ISTANBUL

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

The internal structure of natural numbers and one method for the definition of large prime numbers

The internal structure of natural numbers and one method for the definition of large prime numbers The nternal structure of natural numbers and one method for the defnton of large prme numbers Emmanul Manousos APM Insttute for the Advancement of Physcs and Mathematcs 3 Poulou str. 53 Athens Greece Abstract

More information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

THE WEIGHTED WEAK TYPE INEQUALITY FOR THE STRONG MAXIMAL FUNCTION

THE WEIGHTED WEAK TYPE INEQUALITY FOR THE STRONG MAXIMAL FUNCTION THE WEIGHTED WEAK TYPE INEQUALITY FO THE STONG MAXIMAL FUNCTION THEMIS MITSIS Abstract. We prove the natural Fefferman-Sten weak type nequalty for the strong maxmal functon n the plane, under the assumpton

More information

General viscosity iterative method for a sequence of quasi-nonexpansive mappings

General viscosity iterative method for a sequence of quasi-nonexpansive mappings Avalable onlne at www.tjnsa.com J. Nonlnear Sc. Appl. 9 (2016), 5672 5682 Research Artcle General vscosty teratve method for a sequence of quas-nonexpansve mappngs Cuje Zhang, Ynan Wang College of Scence,

More information

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

More information

Complete subgraphs in multipartite graphs

Complete subgraphs in multipartite graphs Complete subgraphs n multpartte graphs FLORIAN PFENDER Unverstät Rostock, Insttut für Mathematk D-18057 Rostock, Germany Floran.Pfender@un-rostock.de Abstract Turán s Theorem states that every graph G

More information

Estimation: Part 2. Chapter GREG estimation

Estimation: Part 2. Chapter GREG estimation Chapter 9 Estmaton: Part 2 9. GREG estmaton In Chapter 8, we have seen that the regresson estmator s an effcent estmator when there s a lnear relatonshp between y and x. In ths chapter, we generalzed the

More information

Online Appendix. t=1 (p t w)q t. Then the first order condition shows that

Online Appendix. t=1 (p t w)q t. Then the first order condition shows that Artcle forthcomng to ; manuscrpt no (Please, provde the manuscrpt number!) 1 Onlne Appendx Appendx E: Proofs Proof of Proposton 1 Frst we derve the equlbrum when the manufacturer does not vertcally ntegrate

More information

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2 Physcs 607 Exam 1 Please be well-organzed, and show all sgnfcant steps clearly n all problems. You are graded on your wor, so please do not just wrte down answers wth no explanaton! Do all your wor on

More information

Georgia Tech PHYS 6124 Mathematical Methods of Physics I

Georgia Tech PHYS 6124 Mathematical Methods of Physics I Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends

More information

CSCE 790S Background Results

CSCE 790S Background Results CSCE 790S Background Results Stephen A. Fenner September 8, 011 Abstract These results are background to the course CSCE 790S/CSCE 790B, Quantum Computaton and Informaton (Sprng 007 and Fall 011). Each

More information

ON THE EXTENDED HAAGERUP TENSOR PRODUCT IN OPERATOR SPACES. 1. Introduction

ON THE EXTENDED HAAGERUP TENSOR PRODUCT IN OPERATOR SPACES. 1. Introduction ON THE EXTENDED HAAGERUP TENSOR PRODUCT IN OPERATOR SPACES TAKASHI ITOH AND MASARU NAGISA Abstract We descrbe the Haagerup tensor product l h l and the extended Haagerup tensor product l eh l n terms of

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

FINITELY-GENERATED MODULES OVER A PRINCIPAL IDEAL DOMAIN

FINITELY-GENERATED MODULES OVER A PRINCIPAL IDEAL DOMAIN FINITELY-GENERTED MODULES OVER PRINCIPL IDEL DOMIN EMMNUEL KOWLSKI Throughout ths note, s a prncpal deal doman. We recall the classfcaton theorem: Theorem 1. Let M be a fntely-generated -module. (1) There

More information

SL n (F ) Equals its Own Derived Group

SL n (F ) Equals its Own Derived Group Internatonal Journal of Algebra, Vol. 2, 2008, no. 12, 585-594 SL n (F ) Equals ts Own Derved Group Jorge Macel BMCC-The Cty Unversty of New York, CUNY 199 Chambers street, New York, NY 10007, USA macel@cms.nyu.edu

More information

Dirichlet s Theorem In Arithmetic Progressions

Dirichlet s Theorem In Arithmetic Progressions Drchlet s Theorem In Arthmetc Progressons Parsa Kavkan Hang Wang The Unversty of Adelade February 26, 205 Abstract The am of ths paper s to ntroduce and prove Drchlet s theorem n arthmetc progressons,

More information

Lecture 3. Ax x i a i. i i

Lecture 3. Ax x i a i. i i 18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest

More information

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011 Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected

More information

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,

More information

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1 P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

Finding Primitive Roots Pseudo-Deterministically

Finding Primitive Roots Pseudo-Deterministically Electronc Colloquum on Computatonal Complexty, Report No 207 (205) Fndng Prmtve Roots Pseudo-Determnstcally Ofer Grossman December 22, 205 Abstract Pseudo-determnstc algorthms are randomzed search algorthms

More information

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

Professor Terje Haukaas University of British Columbia, Vancouver  The Q4 Element Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

Modelli Clamfim Equazioni differenziali 7 ottobre 2013

Modelli Clamfim Equazioni differenziali 7 ottobre 2013 CLAMFIM Bologna Modell 1 @ Clamfm Equazon dfferenzal 7 ottobre 2013 professor Danele Rtell danele.rtell@unbo.t 1/18? Ordnary Dfferental Equatons A dfferental equaton s an equaton that defnes a relatonshp

More information

Google PageRank with Stochastic Matrix

Google PageRank with Stochastic Matrix Google PageRank wth Stochastc Matrx Md. Sharq, Puranjt Sanyal, Samk Mtra (M.Sc. Applcatons of Mathematcs) Dscrete Tme Markov Chan Let S be a countable set (usually S s a subset of Z or Z d or R or R d

More information

Solutions HW #2. minimize. Ax = b. Give the dual problem, and make the implicit equality constraints explicit. Solution.

Solutions HW #2. minimize. Ax = b. Give the dual problem, and make the implicit equality constraints explicit. Solution. Solutons HW #2 Dual of general LP. Fnd the dual functon of the LP mnmze subject to c T x Gx h Ax = b. Gve the dual problem, and make the mplct equalty constrants explct. Soluton. 1. The Lagrangan s L(x,

More information

Solution Thermodynamics

Solution Thermodynamics Soluton hermodynamcs usng Wagner Notaton by Stanley. Howard Department of aterals and etallurgcal Engneerng South Dakota School of nes and echnology Rapd Cty, SD 57701 January 7, 001 Soluton hermodynamcs

More information

Lecture Note 3. Eshelby s Inclusion II

Lecture Note 3. Eshelby s Inclusion II ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte

More information

A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS

A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS Journal of Mathematcal Scences: Advances and Applcatons Volume 25, 2014, Pages 1-12 A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS JIA JI, WEN ZHANG and XIAOFEI QI Department of Mathematcs

More information

Lecture 5 Decoding Binary BCH Codes

Lecture 5 Decoding Binary BCH Codes Lecture 5 Decodng Bnary BCH Codes In ths class, we wll ntroduce dfferent methods for decodng BCH codes 51 Decodng the [15, 7, 5] 2 -BCH Code Consder the [15, 7, 5] 2 -code C we ntroduced n the last lecture

More information

Remarks on the Properties of a Quasi-Fibonacci-like Polynomial Sequence

Remarks on the Properties of a Quasi-Fibonacci-like Polynomial Sequence Remarks on the Propertes of a Quas-Fbonacc-lke Polynomal Sequence Brce Merwne LIU Brooklyn Ilan Wenschelbaum Wesleyan Unversty Abstract Consder the Quas-Fbonacc-lke Polynomal Sequence gven by F 0 = 1,

More information

In this section is given an overview of the common elasticity models.

In this section is given an overview of the common elasticity models. Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp

More information

On the symmetric character of the thermal conductivity tensor

On the symmetric character of the thermal conductivity tensor On the symmetrc character of the thermal conductvty tensor Al R. Hadjesfandar Department of Mechancal and Aerospace Engneerng Unversty at Buffalo, State Unversty of New York Buffalo, NY 146 USA ah@buffalo.edu

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity

Module 1 : The equation of continuity. Lecture 1: Equation of Continuity 1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum

More information

Eigenvalues of Random Graphs

Eigenvalues of Random Graphs Spectral Graph Theory Lecture 2 Egenvalues of Random Graphs Danel A. Spelman November 4, 202 2. Introducton In ths lecture, we consder a random graph on n vertces n whch each edge s chosen to be n the

More information

Example: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41,

Example: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41, The greatest common dvsor of two ntegers a and b (not both zero) s the largest nteger whch s a common factor of both a and b. We denote ths number by gcd(a, b), or smply (a, b) when there s no confuson

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 151 Lecture 3 Lagrange s Equatons (Goldsten Chapter 1) Hamlton s Prncple (Chapter 2) What We Dd Last Tme! Dscussed mult-partcle systems! Internal and external forces! Laws of acton and

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

CHAPTER 14 GENERAL PERTURBATION THEORY

CHAPTER 14 GENERAL PERTURBATION THEORY CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves

More information

On Finite Rank Perturbation of Diagonalizable Operators

On Finite Rank Perturbation of Diagonalizable Operators Functonal Analyss, Approxmaton and Computaton 6 (1) (2014), 49 53 Publshed by Faculty of Scences and Mathematcs, Unversty of Nš, Serba Avalable at: http://wwwpmfnacrs/faac On Fnte Rank Perturbaton of Dagonalzable

More information

Norm Bounds for a Transformed Activity Level. Vector in Sraffian Systems: A Dual Exercise

Norm Bounds for a Transformed Activity Level. Vector in Sraffian Systems: A Dual Exercise ppled Mathematcal Scences, Vol. 4, 200, no. 60, 2955-296 Norm Bounds for a ransformed ctvty Level Vector n Sraffan Systems: Dual Exercse Nkolaos Rodousaks Department of Publc dmnstraton, Panteon Unversty

More information

Another converse of Jensen s inequality

Another converse of Jensen s inequality Another converse of Jensen s nequalty Slavko Smc Abstract. We gve the best possble global bounds for a form of dscrete Jensen s nequalty. By some examples ts frutfulness s shown. 1. Introducton Throughout

More information

Fixed points of IA-endomorphisms of a free metabelian Lie algebra

Fixed points of IA-endomorphisms of a free metabelian Lie algebra Proc. Indan Acad. Sc. (Math. Sc.) Vol. 121, No. 4, November 2011, pp. 405 416. c Indan Academy of Scences Fxed ponts of IA-endomorphsms of a free metabelan Le algebra NAIME EKICI 1 and DEMET PARLAK SÖNMEZ

More information

Integrals and Invariants of Euler-Lagrange Equations

Integrals and Invariants of Euler-Lagrange Equations Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,

More information

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space.

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space. Lnear, affne, and convex sets and hulls In the sequel, unless otherwse specfed, X wll denote a real vector space. Lnes and segments. Gven two ponts x, y X, we defne xy = {x + t(y x) : t R} = {(1 t)x +

More information

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 21: Numerical methods for pricing American type derivatives Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

LECTURE V. 1. More on the Chinese Remainder Theorem We begin by recalling this theorem, proven in the preceeding lecture.

LECTURE V. 1. More on the Chinese Remainder Theorem We begin by recalling this theorem, proven in the preceeding lecture. LECTURE V EDWIN SPARK 1. More on the Chnese Remander Theorem We begn by recallng ths theorem, proven n the preceedng lecture. Theorem 1.1 (Chnese Remander Theorem). Let R be a rng wth deals I 1, I 2,...,

More information