An Efficient Compact Finite Difference Method for the Solution of the Gross-Pitaevskii Equation

<div>The isosurface plots of first excited state in Example <a href="https://static.hindawi.com/articles/acmp/volume-2015/127580/figures/#ex3">3</a> with <svg height="12.7178pt" id="M201" style="vertical-align:-3.42947pt" version="1.1" viewbox="-0.0498162 -9.28833 64.5412 12.7178" width="64.5412pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M162 -163V703H101V-163H162Z" id="g113-9"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="263" vert-adv-y="263"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,3.419,0)"><path d="M542 242C542 369 469 429 379 468C370 472 372 496 376 513L416 681L400 691L344 648L219 51C167 69 113 128 113 217C113 324 166 387 266 428L252 459C124 418 23 336 23 202C23 95 90 13 205 -8L182 -132C165 -222 191 -254 209 -261L215 -258L270 -7C415 -4 542 98 542 242ZM469 232C469 126 403 33 279 39L355 409C399 394 469 331 469 232Z" id="g113-243"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="565" vert-adv-y="565"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,10.439,3.132)"><path d="M389 0V32C297 37 291 46 291 118V635C234 612 176 595 109 583V556L161 554C203 551 207 548 207 497V118C207 46 201 37 110 32V0H389Z" id="g58-47"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,14.871,3.132)"><path d="M245 635C92 635 37 457 37 312C37 148 92 -12 244 -12C395 -12 449 165 449 312C449 469 394 635 245 635ZM243 598C332 598 358 454 358 312C358 172 334 26 245 26C158 26 128 174 128 313S151 598 243 598Z" id="g58-46"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,19.302,3.132)"><path d="M245 635C92 635 37 457 37 312C37 148 92 -12 244 -12C395 -12 449 165 449 312C449 469 394 635 245 635ZM243 598C332 598 358 454 358 312C358 172 334 26 245 26C158 26 128 174 128 313S151 598 243 598Z" id="g58-46"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,24.282,0)"><path d="M162 -163V703H101V-163H162Z" id="g113-9"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="263" vert-adv-y="263"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,31.333,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,42.596,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g121-46"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,48.836,0)"><path d="M113 -12C146 -12 170 11 170 46C170 78 146 103 114 103S58 78 58 46C58 11 82 -12 113 -12Z" id="g113-47"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="228" vert-adv-y="228"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,51.8,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g121-46"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,58.04,0)"><path d="M384 0V27C293 34 287 42 287 114V635C232 613 172 594 109 583V559L157 557C201 555 205 550 205 499V114C205 42 199 34 109 27V0H384Z" id="g121-47"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg>, <svg height="12.7178pt" id="M202" style="vertical-align:-3.42947pt" version="1.1" viewbox="-0.0498162 -9.28833 64.5412 12.7178" width="64.5412pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M162 -163V703H101V-163H162Z" id="g113-9"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="263" vert-adv-y="263"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,3.419,0)"><path d="M542 242C542 369 469 429 379 468C370 472 372 496 376 513L416 681L400 691L344 648L219 51C167 69 113 128 113 217C113 324 166 387 266 428L252 459C124 418 23 336 23 202C23 95 90 13 205 -8L182 -132C165 -222 191 -254 209 -261L215 -258L270 -7C415 -4 542 98 542 242ZM469 232C469 126 403 33 279 39L355 409C399 394 469 331 469 232Z" id="g113-243"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="565" vert-adv-y="565"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,10.439,3.132)"><path d="M389 0V32C297 37 291 46 291 118V635C234 612 176 595 109 583V556L161 554C203 551 207 548 207 497V118C207 46 201 37 110 32V0H389Z" id="g58-47"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,14.871,3.132)"><path d="M389 0V32C297 37 291 46 291 118V635C234 612 176 595 109 583V556L161 554C203 551 207 548 207 497V118C207 46 201 37 110 32V0H389Z" id="g58-47"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,19.302,3.132)"><path d="M245 635C92 635 37 457 37 312C37 148 92 -12 244 -12C395 -12 449 165 449 312C449 469 394 635 245 635ZM243 598C332 598 358 454 358 312C358 172 334 26 245 26C158 26 128 174 128 313S151 598 243 598Z" id="g58-46"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,24.282,0)"><path d="M162 -163V703H101V-163H162Z" id="g113-9"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="263" vert-adv-y="263"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,31.333,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,42.596,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g121-46"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,48.836,0)"><path d="M113 -12C146 -12 170 11 170 46C170 78 146 103 114 103S58 78 58 46C58 11 82 -12 113 -12Z" id="g113-47"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="228" vert-adv-y="228"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,51.8,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g121-46"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,58.04,0)"><path d="M384 0V27C293 34 287 42 287 114V635C232 613 172 594 109 583V559L157 557C201 555 205 550 205 499V114C205 42 199 34 109 27V0H384Z" id="g121-47"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg>, and <svg height="12.7178pt" id="M203" style="vertical-align:-3.42947pt" version="1.1" viewbox="-0.0498162 -9.28833 64.5412 12.7178" width="64.5412pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M162 -163V703H101V-163H162Z" id="g113-9"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="263" vert-adv-y="263"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,3.419,0)"><path d="M542 242C542 369 469 429 379 468C370 472 372 496 376 513L416 681L400 691L344 648L219 51C167 69 113 128 113 217C113 324 166 387 266 428L252 459C124 418 23 336 23 202C23 95 90 13 205 -8L182 -132C165 -222 191 -254 209 -261L215 -258L270 -7C415 -4 542 98 542 242ZM469 232C469 126 403 33 279 39L355 409C399 394 469 331 469 232Z" id="g113-243"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="565" vert-adv-y="565"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,10.439,3.132)"><path d="M389 0V32C297 37 291 46 291 118V635C234 612 176 595 109 583V556L161 554C203 551 207 548 207 497V118C207 46 201 37 110 32V0H389Z" id="g58-47"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,14.871,3.132)"><path d="M389 0V32C297 37 291 46 291 118V635C234 612 176 595 109 583V556L161 554C203 551 207 548 207 497V118C207 46 201 37 110 32V0H389Z" id="g58-47"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,19.302,3.132)"><path d="M389 0V32C297 37 291 46 291 118V635C234 612 176 595 109 583V556L161 554C203 551 207 548 207 497V118C207 46 201 37 110 32V0H389Z" id="g58-47"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,24.282,0)"><path d="M162 -163V703H101V-163H162Z" id="g113-9"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="263" vert-adv-y="263"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,31.333,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,42.596,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g121-46"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,48.836,0)"><path d="M113 -12C146 -12 170 11 170 46C170 78 146 103 114 103S58 78 58 46C58 11 82 -12 113 -12Z" id="g113-47"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="228" vert-adv-y="228"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,51.8,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g121-46"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,58.04,0)"><path d="M384 0V27C293 34 287 42 287 114V635C232 613 172 594 109 583V559L157 557C201 555 205 550 205 499V114C205 42 199 34 109 27V0H384Z" id="g121-47"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg>.</div>

Advances in Condensed Matter Physics

fig5

Figure 5

Figure 5: An Efficient Compact Finite Difference Method for the Solution of the Gross-Pitaevskii Equation 